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HUMAN GENOME -DERIVED SINGLE EXON NUCLEIC ACID PROBES USEFUL 
FOR ANALYSIS OF GENE EXPRESSION IN HUMAN FETAL LIVER 



CROSS REFERENCE TO RELATED APPLICATIONS 

5 

The present application is a continuation-in-part of U.S. 
patent application serial nos , 09/632,366, filed August 3, 
2000 and 09/608,408, filed June 30, 2000; claims the 
benefit under 35 U.S.C. s 119(e) of U.S . provisional patent 

10 application serial nos. 60/236;, 3 59 , filed September 27, 
2000, 60/234,687, filed September 21, 2000, 60/207,456, 
filed May 26, 2000, and 60/180,312, filed February 4, 2000; 
and further claims the benefit under 35 U.S.C. s 119(a) of 
UK patent application no. 0024263.6, filed October 4, 2000, 

15 the disclosures of which are incorporated herein by 
reference in their entireties. 

REFERENCE TO SEQUENCE LISTING AND INCORPORATION BY 
REFERENCE THEREOF 

20 

The present application includes a Sequence Listing in 
electronic format, filed pursuant to PCT Administrative 
Instructions 801 - 806 on a single CD-R disc, in 
triplicate, containing a file named pto_FETAL_LIVER , txt , 
25 created 24 January 2001, having 25,630,231 bytes. The 
Sequence Listing contained in said file on said disc is 
incorporated herein by reference in its entirety. 

Field of the Invention 

30 

The present invention relates to genome -derived 
single exon microarrays useful for verifying the expression 
of regions of genomic DNA predicted to encode protein. In 
particular, the present invention relates to" unique genome- 
35 derived single exon nucleic acid probes expressed in human 

1 
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Fetal liver and single exon nucleic acid microarrays that 
include such probes. 



Background of the Invention 
5 For almost two decades following the invention of 

general techniques for nucleic acid sequencing, Sanger et 
al., Proc.Natl. Acad. Sci . USA 70 (4 ): 1209 - 13 (1973); 
Gilbert et al . , Proc , Natl. Acad. Sci. USA 70(12) :3581-4 
(1973) , these techniques were used principally as tools to 
10 further the understanding of proteins — known or 

suspected - about which a basic foundation of biological 
knowledge had already been built. In many cases, the 
cloning effort that preceded sequence identification had_ 
been both informed and directed by that antecedent 
15 biological understanding. 

For example, the cloning of the T cell receptor 
for antigen was predicated upon its known or suspected cell 
type-specific expression, by its suspected membrane 
association, and by the predicted assembly of its gene via 
20 T cell- specif ic somatic recombination. Subsequent 
sequencing efforts at once confirmed and extended 
understanding of this family of proteins. Hedrick et al . , 
Na.ture 308 (5955 ): 153 -8 (1984). 

More recently, however, the development of high 
25 throughput sequencing methods and devices, in concert with 
large public and private undertakings to sequence the human 
and other genomes, has altered this investigational 
paradigm: today, sequence information often precedes 
understanding of the basic biology of the encoded protein 
30 product . 

One of the approaches to large-scale sequencing 
is predicated upon the proposition that expressed \ 
sequences - that is, those accessible through isolation of- 
mRNA - are of greatest initial interest. This "expressed , 
35 sequence tag" ("EST") approach has already yielded vast 

2 
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amounts of sequence data (see for example Adams et al . , 
Science 252:1651 (1991); Williamson, Drug Discov. Today 
4:115 (1999)) . For nucleic acids sequenced by this 
approach, often the only biological information that is 
5 known a priori with any certainty is the likelihood of 
biologic expression itself. By virtue of the species and 
tissue from which the mRNA had originally been obtained, 
most such sequences are also annotated with the identity of 
the species and at least one tissue in which expression 
10 appears likely. 

More recently, the pace of genomic sequencing has 
accelerated dramatically.' When genomic DNA serves as the 
initial substrate for sequencing efforts, expression cannot 
be presumed; often the only a priori biological information 
15 about the sequence includes the species and chromosome (and 
perhaps chromosomal map location) of origin. 

With the ever-accelerating pace of sequence 
accumulation by directed, EST, and genomic sequencing 
approaches - and in particular, with the accumulation of 
20 sequence information from multiple genera, from multiple 

species within genera, and from multiple individuals within 
a species - there is an increasing need for methods that 
rapidly and effectively permit the functions of nucleic 
sequences to be elucidated. And as such functional 
25 information accumulates, there is a further need for 
methods of storing such functional information in 
meaningful and useful relationship to the sequence itself; 
that is, there is an increasing need for means and 
apparatus for annotating raw sequence data with known or 
30 predicted functional information. 

Although the increase in the pace of genomic 
sequencing is due in large part to technological changes in 
sequencing strategies and instrumentation. Service, Science 
280:995 (1998); Pennisi, Science 283: 1822-1823 (1999), 
35 there is an important functional motivation as well. 

3 
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While it was understood that the EST approach 
would rarely be able to yield sequence information about 
the noncoding portions of the genome, it now also appears 
the EST approach is capable of capturing only a fraction of 

5 a genome's actual expression complexity. 

For example, when the C. elegrans genome was fully 
sequenced, gene prediction algorithms identified over 
19,000 potential genes, of which only 7,000 had been found 
by EST sequencing. C. elegans Sequencing Consortium, 

10 Science 282:2012 (1998). Analogously, the recently 

completed sequence of chromosome 2 of Arabidopsis predicts 
over 4000 genes, Lin et al . , Nature, 402:761 (1999), of 
which only about 6% had- previously been identif ied via EST 
sequencing efforts. Although the human genome has the 

15 greatest depth of EST coverage, it is still woefully short 
of surrendering all of its genes. One recent estimate 
suggests that the human genome contains more than 146,000 
genes, which would at this point leave greater than half of 
the genes undiscovered. It is now predicted that many 

20 genes, perhaps 20 to 50%, will only be found by genomic 
sequencing. 

There is, therefore, a need for methods that 
permit the functional regions of genomic sequence - and 
most importantly, but not exclusively, regions that 

25 function to encode genes - to be identified. 

Much of the coding sequence of the human genome 
is not homologous to known genes, making detection of open 
reading frames ("ORFs") and predictions of gene function 
difficult. Computational methods exist for predicting 

30 coding regions in eukaryotic genomes. Gene prediction 
programs such as GRAIL and GRAIL II, Uberbacher et al . , 
Proc. Natl, Acad, Sci . USA 88 (24) ; 11261-5 (1991); Xu et 
al., Genet. Eng. 16:241-53 (1994); Uberbacher et ai . , 
Methods Enzymol. 266:259-81 (1996); GENEFINDER, Solovyev et 

35 ai., Nucl, Acids. Res. 22:5156-63 (1994); Solovyev et al . , 
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Ismb 5:294-302 (1997); and GENESCAN, Burge et al . , J. Mol . 
Biol. 268:78-94 (1997), predict many putative genes without 
known homology or function. Such programs are known, 
however, to give high false positive rates. Burset et al., 
5 Genomics 34:353-367 (1996). Using a consensus obtained by 
a plurality of such programs is known to increase the 
reliability of calling exons from genomic sequence. 
Ansari-Lari et al , , Genome Res. 8(1): 29-40 (1998) 

Identification of functional genes from genomic 
10 data remains, however, an imperfect art. For example, in 
reporting the full sequence of human chromosome 21, the 
Chromosome 21 Mapping and Sequencing Consortium reports 
that prior bioinf ormatic estimates of human gene number may 
need to be revised substantially downwards. Nature 
15 405:311-199 (2000); Reeves, Nature 405:283-284 (2000). 

Thus, there is a need for methods and apparatus 
that permit the functions of the regions identified 
bioinformatically - and specifically, that permit the 
expression of regions predicted to encode protein - readily 
20 to be confirmed experimentally. 

Recently, the development of nucleic acid 
microarrays has made possible the automated and highly 
parallel measurement of gene expression. Reviewed in 
Schena (ed.), DNA Microarrays : A Practical Approach 
25 (Practical Approach Series ) , Oxford University Press (1999) 
(ISBN: 0199637768); Nature Genet. 21 (1) (suppl ) : 1 - 60 
(1999); Schena (ed.), Microarray Biochip: Tools and 
Technology , Eaton Publishing Company/BioTechniques Books 
Division (2000) (ISBN: 1881299376) . 
30 It is common for microarrays to be derived from 

cDNA/EST libraries, either from those previously described 
in the literature, such as those from the I.M.A.G.E. 
consortium, Lennon et al . , Genomics 33(1): 151-2 (1996), or 
from the construction of "problem specific" libraries 
35 targeted at a particular biological question, R.S. Thomas 
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et al., Cancer Res. (in press) . Such microarrays by 
definition can measure expression only of those genes found 
in EST libraries, and thus have not been useful as probes 
for genes discovered solely by genomic sequencing. 

The utility of using whole genome nucleic acid 
microarrays to answer certain biological questions has been 
demonstrated for the yeast Saccharomyces cerevisiae, De 
Risi et al.. Science 278:680 (1997). The vast majority of 
yeast nuclear genes, approximately 95% however, are single 
exon genes, i.e., lack introns, Lopez et al . , RNA 5:1135- 
1137 (1999); Goffeau etal., Science 274 : 563 -67 (1996), 
permitting coding regions more readily to be identified. 
Whole. genome nucleic, acid .microarrays have not generally 
been used to probe gene expression from more complex 
eukaryotic genomes, and in particular from those averaging 
more than one intron per gene. 

Diseases of the liver are a significant cause of 
human morbidity and mortality. Increasingly, genetic 
factors are being found that contribute to predisposition, 
onset, and/or aggressiveness of most, if not all, of these 
diseases; although causative mutations in single genes have 
been identified for some, these disorders are believed for 
the most part to have polygenic etiologies. There is a need 
for methods and apparatus that permit prediction of 
diseases of the liver, particularly those diseases with 
polygenic etiology, from diagnosis of fetal liver. 



Summary of the Invention 

30 The present invention solves these and other 

. problems in the art by providing methods and apparatus for 
predicting, confirming, and displaying functional \ 
information derived from genomic sequence. The present 
invention also provides apparatus for verifying the 

35 expression of putative genes identified within genomic 

6 
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sequence . 

In particular, the invention provides novel 
genome-derived single exon nucleic acid microarrays useful 
for verifying the expression of putative genes identified 
5 within genomic sequence. 

The present invention also provides compositions 
and kits for the ready production of nucleic acids 
identical in sequence to, or substantially identical in 
sequence to, probes on the genome-derived single exon 
10 microarrays of the present invention. 

Accordingly, in a first aspect of the invention, 
there is provided a spatially-addressable set of single 
exon nucleic acid probes for measuring gene expression in a 
sample derived from human Fetal liver, comprising a 
15 plurality of single exon nucleic acid probes according to 
any one of the nucleotide sequences set out in SEQ ID NOs : 
1 - 12,673 or a complementary sequence, or a portion of 
such a sequence. 

By plurality is meant at least two, suitably at 
20 least 20, most suitably at least 100, preferably at least 
1000 and, most preferably, upto 5000, 

In one embodiment of the first aspect, each of ^ 
said plurality of probes is separately and addressably 
amplifiable. 

25 In an alternative embodiment, each of said 

plurality of probes is separately and addressably 
isolatable from said plurality. 

In a preferred embodiment, each of said plurality 
of probes is amplifiable using at least one common primer. 
Preferably, each of said plurality of probes is amplifiable 
using a first and a second common primer. 

In yet another embodiment, said set of single 
exon nucleic acid probes comprises between 50 - 20,000 
probes, for example, 50 - 5000. 

Suitably, said set of single exon nucleic acid 
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probes comprises at least 50 - 1000 discrete single exon 
nucleic acid probes having a sequence as set out in any of 
SEQ ID NOS.: 1 - 25,129 or a complimentary sequence, or a 
portion of such a sequence. 

Preferably, the average length of the single exon 
nucleic acid probes is between 200 and 500 bp. It is 
preferred that the average length should be at least 200bp, 
suitably at least 250bp, most suitably at least 300bp, 
preferably at least 400bp and, most preferably, 500 bp. 

In another embodiment, the single exon nucleic 
acid probes lack prokaryotic and bacteriophage vector 
sequence. It is preferred that at least 50%, suitably at 
least 6,0%,, most suitably. at least, 70%, preferably at least 
75%, more preferably at least 80, 85, 90, 95 or 99% of said 
single exon nucleic acid probes lack prokaryotic and 
bacteriophage vector sequence. 

In another preferred embodiment, said single exon 
nucleic acid lack homopolymeric stretches of A or T. It is 
preferred that at least 50%, suitably at least 60%, most 
suitably at least 70%, preferably at least 75%, more 
preferably at least 80, 85, 90, 95 or 99% of said single 
exon nucleic acid probes lack homopolymeric stretches of A 
or T. 

Preferably, a spat ially- addressable set of single 
25 exon nucleic .acid probes in accordance with the first 

aspect of the invention is is addressably disposed upon a 
substrate . 

Suitable substrates include a filter membrane 
which may, preferably, be nitrocellulose or nylon. The 
nylon may preferably, be positively-charged. Other suitable 
substrates include glass, amorphous silicon, crystalline 
silicon, and plastic. Further suitable materials include 
polymethylacrylic, polyethylene, polypropylene, 
polyacrylate, polymethylmethacrylate, polyvinylchloride , 
polytetrafluoroethylene, polystyrene, polycarbonate, 



20 



30 



35 



8 



20 



25 



WO 01/57277 ; PCT/U SO 1/00669 

polyacetal, polysulfone, celluloseacetate , 
cellulosenitrate, nitrocellulose, and mixtures thereof. 

in a second aspect of the invention, there is 
provided a microarray comprising a spatially addressable 
set of single exon nucleic acid probes in accordance with 
the first aspect of the invention. 

in one embodiment, a genome -derived single -exon 
microarray is packaged together with such an ordered set of 
amplifiable probes corresponding to the probes, or one or 
more subsets of probes, thereon. In alternative 
embodiments, the ordered set of amplifiable probes is 
packaged separately from the genome -derived single exon 
microarray. 

In another aspect, the invention provides genome- 
; derived single exon nucleic acid probes useful for gene 
expression analysis, and particularly for gene expression 
analysis by microarray. In particular embodiments of this 
aspect, the present invention provides human single -exon 
probes that include specif ically-hybridizable fragments of 
SEQ ID Nos, 12,674 - 25,129, wherein the fragment 
hybridizes at high stringency to an expressed human gene, 
in particular embodiments, the invention provides single 
exon probes comprising SEQ ID Nds. 1 - 12,673. 

Accordingly, in a third aspect of the invention, 
there is provided a single exon nucleic acid probe for 
measuring human gene expression in a sample derived from 
human Fetal liver which is a nucleic acid molecule 
comprising a nucleotide sequence as set out in any of SEQ 
ID NOs.: 1 - 12,673 or a complementary sequence or a 
fragment thereof wherein said probe hybridizes at high 
stringency to a nucleic acid expressed in the human Fetal 
liver. 

In one embodiment, a single exon nucleic acid 
probe in accordance with the third aspect comprises a 
35 nucleotide sequence as set out in any of SEQ ID NOs. : 
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12,674 - 25,129 or a complementary sequence or a fragment 
thereof . 

In a fourth aspect of the invention, there is 
provided a single exon nucleic acid probe for measuring 

5 human gene expression in a sample derived from human Fetal 
liver which is a nucleic acid molecule having a sequence 
encoding a peptide comprising a peptide sequence as set out 
in any of SEQ ID NOs . : 25,130 - 37,156 or a complementary 
sequence or a fragment thereof wherein said probe 

iO hybridizes at high stringency to a nucleic acid expressed 
in the human Fetal liver. 

Preferably, a single exon nucleic acid probe in 
accordance with the third or fourth aspects of .the 
invention comprises between at least 15 and 50 contiguous 

15 nucleotides of said SEQ ID NO:. It is preferred that the 
single exon nucleic acid probe comprises at least 15, 
suitably at least 20, more suitably at least 25 or 
preferably at least 50 contiguous nucleotides of said SEQ 
ID NO : . 

20 In another preferred embodiment, a single exon 

nucleic acid probe in accordance with the third or fourth 
aspects of the invention is between 3kb and 25kb in length. 
It is preferred that said probe is no more than 3kb, 
suitably no more than 5kb, more suitably no more than lOkb, 

25 preferably 15kb, more preferably 20kb or, most preferably, 
no more than 20kb in length. 

Preferably, a single exon nucleic acid probe in 
accordance with either the fifth or sixth aspect of the 
invention is DNA, preferably single -stranded DNA, RNA or 

30 PNA. 

In another embodiment of either the third or 
fourth aspect of the invention, a single exon nucleic acid 
probe is detectably labeled. Suitable detectable labels 
include a radionuclide, a fluorescent label or a first 
35 member of a specific binding pair. Suitable fluorescent 
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labels include dyes such as cyanine dyes, preferably Cy3 
and Cy5 although other suitable dyes will be known to those 
skilled in the art . 

In a particularly preferred embodiment, a single 
5 exon nucleic acid probe in accordance with either the third 
or fourth aspect of the invention lacks prokaryotic and 
bacteriophage vector sequence. In yet another embodiment, a 
single exon nucleic acid probe in accordance with either 
the third or fourth aspect of the invention lacks 
10 homopolymeric stretches of A or T. 

In a fifth aspect of the invention, there is 
provided an amplifiable nucleic acid composition, 
comprising : 

the single exon nucleic acid probe in accordance 
15 with either of the third or fourth aspects of the 
invention; and at least one nucleic acid primer; 

wherein said at least one primer is sufficient to 
prime enzymatic amplification of said probe. 

In an sixth aspect of the invention, there is 
20 provided a method of measuring gene expression in a sample 
derived from human Fetal liver, comprising: 

contacting the single exon microarray in 
accordance with the second aspect of the invention, with a 
first collection of detectably labeled nucleic acids, said 
25 first collection of nucleic acids derived from mRNA of 
human Fetal liver; and then 

measuring the label detectably bound to each 
probe of said microarray. 

In a seventh aspect of the invention, there is 
30 provided a method of identifying exons in a eukaryotic 
genome, comprising: 

algorithmically predicting at least one exon from 
genomic sequence of said eukaryote; and then 

detecting specific hybridization of detectably 
35 labeled nucleic acids to a single exon probe, 

11 
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wherein said detectably labeled nucleic acids are 
derived from mRNA from- the Fetal liver of said eukaryote, 
said probe is a single exon probe having a fragment 
identical in sequence to, or complementary in sequence to, 

5 said predicted exon, said probe is included within a single 
exon microarray in accordance with the first aspect of the 
invention, and said fragment is selectively hybridizable at 
high stringency. 

In a eighth aspect of the invention, there is 

10 provided a method of assigning exons to a single gene, 
comprising : 

identifying a plurality of exons from genomic 

sequence in accordance with the seventh aspect of the 

invention; and then 
15 measuring the expression of each of said exons in 

a plurality of tissues and/or cell types using 

hybridization to single exon microarrays having a probe 

with said exon, 

wherein a common pattern of expression of said 
20 exons in said plurality of tissues and/or cell types 

indicates that the exons should be assigned to a single 

gene . 

In an ninth aspect of the invention, there is 
provided a nucleic acid sequence as set out in any of SEQ 
25 ID NOs: 1 - 25,129 wherein said sequence encodes a peptide. 

In a tenth aspect of the invention, there is 
provided a peptide encoded by a sequence comprising a 
sequence as set out in any of SEQ ID NOs: 12,674 - 25,129, 
or a complementary sequence or coding portion thereof. 
30 In a preferred embodiment, a peptide may be 

encoded by a sequence comprising a sequence set out in any 
of SEQ ID NOS.: 1 -12,673. \ 

In a further aspect, the invention provides 
peptides comprising an amino acid sequence translated from 
35 the DNA fragments, said amino acid sequences comprising SEQ 
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ID NOS. : 25,130 - 37,156. 

Accordingly in a eleventh aspect of the invention 
there is provided a peptide comprising a sequence as set 
out in any of SEQ ID NOs: 25,130 - 37,156, or fragment 
5 thereof. 

In another aspect, the invention provides means 
for displaying annotated sequence, and in particular, for 
displaying sequence annotated according to the methods and 
apparatus of the present invention. Further, such display 
10 can be used as a preferred graphical user interface for 
electronic search, query, and analysis of such annotated 
sequence . 

15 Detailed Description of the Invention 
Definitions 

As used herein, the term "microarray" and phrase 
"nucleic acid microarray" refer to a substrate-bound 

20 collection of plural nucleic acids, hybridization to each 
of the plurality of bound nucleic acids being separately 
detectable. The substrate can be solid or porous, planar 
or non-planar, unitary or distributed. 

As so defined, the term "microarray" and phrase 

25 "nucleic acid microarray" include all the devices so called 
in Schena (ed.), DNA Microarrays: A Practical Approach 
(Practical Approach Series ) , Oxford University Press (1999) 
(ISBN: 0199637768); Nature Genet. 21 (1) (suppl) : 1 - 60 
(1999); and Schena (ed.), Microarray Biochip: Tools and 

30 Technology , Eaton Publishing Company/BioTechniques Books 
Division (2000) (ISBN: 1881299376) . As so defined, the 
term "microarray" and phrase "nucleic acid microarray" 
further include substrate -bound collections of plural 
nucleic acids in which the nucleic acids are distributably 

35 disposed on a plurality of beads, rather than on a unitary 

13 
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planar substrate, as is described, inter alia, in Brenner 
et al., Proc, Natl, Acad. Sci. USA 97 (4) : 166501670 (2000); 
in such case, the term "microarray" and phrase "nucleic 
acid microarray" refer to the plurality of beads in 

5 aggregate. 

As used herein with respect to a nucleic acid 
microarray, the term "probe" refers to the nucleic acid 
that is, or is intended to be, bound to the substrate; in 
such context, the term "target" thus refers to nucleic acid 

10 intended to be bound thereto by Watson-Crick 

complementarity. As used herein with respect to solution 
phase hybridization, the term "probe" refers to the nucleic 
acid of known sequence that is detectably.Jabeled . 

As used herein, the expression "probe comprising 

15 SEQ ID NO,", and variants thereof, intends a nucleic acid 
probe, at least a portion of which probe has either (i) the 
sequence directly as given in the referenced SEQ ID NO., or 
(ii) a sequence complementary to the sequence as given in 
the referenced SEQ ID NO., the choice as between sequence 

20 directly as given and complement thereof dictated by the 
requirement that the probe hybridize to mRNA. 

As used herein, the term "open reading frame" and 
the equivalent acronym "ORF" refer to that portion of an 
exon that can be translated in its entirety into a sequence 

25 of contiguous amino acids i.e. a nucleic acid sequence 

that, in at least one reading frame, does not possess stop 
codons; the term does not require that the ORF encode the 
entirety of a natural protein. 

As used herein, the term ''amplicon" refers to a 

30 PGR product amplified from human genomic DNA, containing 
the predicted exon. 

As used herein the' term '"exon" refers to the 
consensus prediction of the various exon and gene 
predicting algorithms i.e. a nucleic acid sequence 

35 bioinf ormatically predicted to encode a portion of a 

14 
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natural protein. 

As used herein, the term "peptide" refers to a 
sequence of amino acids. The sequences referred to as 
PEPTIDE SEQ ID NOS . : are thq predicted peptide sequences 

5 that would be translated from one of the exons, or a 

portion thereof set out in exon SEQ ID NOS,:. The codons 
encoding the peptide are wholly contained within the exon. 

As used herein, a "portions" of a defined 
nucleotide sequence or sequences can be and, preferably, 

10 are fragments unique to that sequence or to one or a 
combination of those sequences. A fragment unique to a 
nucleic acid'molecule is one that is a signature for the 
larger nucleic acid molecule. 

As used herein, the phrase "expression of a 

15 probe" and its linguistic variants means that the ORE 
present within the probe,' or its complement, is present 
within a target mRNA. 

As used herein, ''stringent conditions" refers to 
parameters well known to those skilled in the art. When a 

20. nucleic acid molecule is said to be hybridisable to another 
of a given sequence under ''stringent conditions" it is 
meant that it is homologous to the given sequence. 

As used herein, the phrase "specific binding 
pair" intends a pair of molecules that bind to one another 

25 with high specificity. Binding pairs are said to exhibit 
specific binding when they exhibit avidity of at least 10^, 
preferably at least 10^, more preferably at least 10^ 
liters/mole. Nonlimiting examples of specific binding 
pairs are: antibody and antigen; biotin and avidin,- and 

30 biotin and streptavidin . 

As used herein with respect to the visual display 
of annotated genomic sequence, the term "rectangle" means 
any geometric shape that has at least a first and a second 
border, wherein the first and second borders each are 

35 capable of mapping uniquely to a point of another visual 

15 
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object of the display. 

As used herein, a "Mondrian" means a visual 
display in which a single genomic sequence is annotated 
with predicted and experimentally confirmed functional 
5 information. 



Brief Description of the Drawings 

10 The present invention is further illustrated with 

reference to the following non-limiting figures and 
examples in which: 

FIG. 1 illustrates a process for predicting _ ^ 
functional regions from genomic sequence, confirming the 

15 functional activity of such regions experimentally, and 
associating and displaying the data so obtained in 
meaningful and useful relationship to the original sequence 
data; 

FIG. 2 further elaborates that portion of the 
20 process schematized in FIG. 1 for predicting functional 
regions from genomic sequence; 

FIG. 3 illustrates a Mondrian visual display; 

FIG. 4 presents a Mondrian showing a hypothetical 
annotated genomic sequence; 
25 FIG. 5 is a histogram showing the distribution of 

ORF length and PCR products as obtained, with ORF length 
shown in black and PCR product length shown in dotted 
lines; 

FIG. 6 is a histogram showing the distribution, 
30 among exons predicted according to the methods described, 
of expression as measured using simultaneous two color 
hybridization to a genome -derived single exon microarray. 
The graph shows the number of sequence -verified products 
that were either not expressed {"0"), expressed in one or 
35 more but not all tested tissues ("1" - "9"), or expressed 
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in all tissues tested ("10"); 

FIG. 7 is a pictorial representation of the 
expression of verified sequences that showed expression 
with signal intensity greater than 3 in at least one 
5 tissue, with: FIG. 7A showing the expression as measured by 
microarray hybridization in each of the 10 measured 
tissues, and the expression as measured "bioinf ormatically" 
by query of EST, NR and SwissProt databases; with FIG. 7B 
showing the legend for display of physical expression 
10 (ratio) in FIG. 7A; and with FIG. 7C showing the legend for 
scoring EST hits as depicted in FIG. 7A; 

FIG. 8 shows a comparison of normalized CY3 
signal intensity for arrayed sequences that were identical 
to sequences in existing EST, NR and SwissProt databases or 
15 that were dissimilar (unknown) , where black denotes the 
signal intensity for all sequence-verified products with a 
BLAST Expect ("E") value of greater than le-30 (1 x 10'^°) 
("unknown") and a dotted line denotes sequence-verified 
spots with a BLAST expect ("E") value of less than le-30 (1 
20 X 10'^°) ("known") ; 

FIG. 9 presents a Mondrian of BAG AC008172 (bases 
25,000 to 130,000), containing the carbamyl phosphate 
synthetase gene (AF154830 . 1) ; and 

FIG. 10 is a Mondrian of BAG A04 983 9. 

25 

Methods and Apparatus for Predicting, Confirming, 
Annotating, and Displaying Functional Regions From Genomic 
Sequence Data 

30 

FIG. 1 is a flow chart illustrating in broad 
outline a process for predicting functional regions from 
genomic sequence, confirming and characterizing the 
functional activity of such regions experimentally, and 
35 then associating and displaying the information so obtained 

17 
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in meaningful and useful relationship to the original 
sequence data. 

The initial input into process 10 of the present 
invention is drawn from one or more databases 100 
5 containing genomic sequence data. Because genomic sequence 
is usually obtained from subgenomic fragments, the sequence 
data typically will be stored in a series of records 
corresponding to these subgenomic sequenced fragments. 
Some fragments will have been catenated to form larger 

10 contiguous sequences ("contigs"); others will not. A 
finite percentage of sequence data in the database will 
typically be erroneous, consisting inter alia of vector 
sequence, sequence created from aberrant cloning. events, . 
sequence of artificial polylinkers, and sequence that was 

15 erroneously read. 

Each sequence record in database 100 will 
minimally contain as annotation a unique sequence 
identifier {accession number) , and will typically be 
annotated further to identify the date of accession, 

20 species of origin, and depositor. Because database 100 can 
contain nongenomic sequence, each sequence will typically 
be annotated further to permit query for genomic sequence . 
Chromosomal origin, optionally with map location, can also 
be present. Data can be, and over time increasingly will 

25 be, further annotated with additional information, in part 
through use of the present invention, as described below. 
Annotation can be present within the data records, in 
information external to database 100 and linked to the 
records thereto, or through a combination of the two. 

30 Databases useful as genomic sequence database 100 

in the present invention include GenBank, and particularly 
include several divisions thereof, including the \ 
htgs{draft), NT {nucleotide, command line), and NR 
{nonredundant ) divisions. GenBank is produced by the 

35 National Institutes of Health and is maintained by the 
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National Center for Biotechnology Information (NCBI) . 
Databases of genomic sequence from species other than 
human, such as mouse, rat, Arabidopsis, C. elegans, C. 
hrigsii, Drosophila, zebra fish, and other higher 
5 eukaryotic organisms will also prove useful as genomic 
sequence database 100 . 

Genomic sequence obtained by query of genomic 
sequence database 100 is then input into one or more 
processes 200 for identification of regions therein that 

10 are predicted to have a biological function as specified by 
the user. Such functions include, but are not limited to, 
encoding protein, regulating transcription, regulating 
message transport after transcription into mRNA, regulating 
message splicing after transcription into mRNA, of 

15 regulating message degradation after transcription into 
mRNA, and the like. Other functions include directing 
somatic recombination events, contributing to chromosomal 
stability or movement, contributing to allelic exclusion or 
X chromosome inactivation , and the like. 

20 The particular genomic sequence to be input into 

process 200 will depend, upon the function for which 
relevant sequence is to be identified as well as upon the_ 
approach chosen for such identification. Process step 200 
can be iterated to identify different functions within a 

25 given genomic region. In such case, the input often will 
be different for the several iterations. 

Sequences predicted to have the requisite 
function by process 200 are then input into process 300, 
where a subset of the input sequences suitable for 

30 experimental confirmation is identified. Experimental 
confirmation can involve physical and/or bioinf ormatic 
assay. Where the subsequent experimental assay is 
bioinf ormatic, rather than physical, there are fewer 
constraints on the sequences that can be tested, and in 

35 this latter case therefore process 300 can output the 

19 
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entirety of the input sequence. 

The subset of sequences output from process 3 00 
is then used in process 400 for experimental verification 
and characterization of the function predicted in 
5 process 200, which experimental verification can, and often 
will, include both physical and bioinf ormat ic assay. 

Process 500 annotates the sequence data with the 
functional information obtained in the physical and/or 
bioinf ormatic assays of process 400. Such annotation can 
10 be done using any technique that usefully relates the 

functional information to the sequence, as, for example, by 
incorporating the functional data into the sequence data 
record itself, by linking records in a hierarchical or 
relational database, by linking to external databases, by a 
15 combination thereof, or by other means well known within 
the database arts. The data can even be submitted for 
incorporation into databases maintained by others, such as 
GenBank, which is maintained by NCBI . 

As further noted in FIG. 1, additional annotation 
20 can be input into process 500 from external sources 600. 

The annotated data is then displayed in process 
800, either before^ concomitantly with, or after optional 
storage 700 on nontransient media, such as magnetic disk, 
optical disc, magnetoopt ical disk, flash memory, or the 
25 like. 

FIG. 1 shows that the experimental data output 
from process 400 can be used in each preceding step of 
process 10: e.g., facilitating identification of functional 
sequences in process 200, facilitating identification of an 

30 experimentally suitable subset thereof in process 300, and 
facilitating creation of physical and/or informational 
substrates for, and performance of subsequent assay, of 
functional sequences in process 400. 

Information from each step can be passed directly 

35 to the succeeding process, or stored in permanent or 
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interim form prior to passage to the succeeding process. 
Often, data will be stored after each, or at least a 
plurality, of such process steps. Any or all process steps 
can be automated. i 

5 FIG. 2 further elaborates the prediction of 

functional sequence within genomic sequence according to 
process 200. 

Genomic sequence database 100 is first queried 20 
for genomic sequence. 
10 The sequence required to be returned by query 20 

will depend, in the first instance, upon the function to be 
identified. 

For example, genomic sequences that function to 
encode protein can be identified inter alia using gene 

15 prediction approaches, comparative sequence analysis 

approaches, or combinations of the two. In gene prediction 
analysis, sequence from one genome is input into process 
200 where at least one, preferably a plurality, of 
algorithmic methods are applied to identify putative coding 

20 regions. In comparative sequence analysis, by contrast, 

corresponding, e.g., syntenic, sequence from a plurality of 
sources, typically a plurality of species, is input into 
process 200, where at least one, possibly a plurality, of 
algorithmic methods are applied to compare the sequences 

25 and identify regions of least variability. 

The exact content of query 20 will also depend 
upon the database queried. For example, if the database 
contains both genomic and nongenomic sequence, perhaps 
derived from multiple species, and the function to be 

30 determined is protein coding regions in human genomic 
sequence, the query will accordingly require that the 
sequence returned be genomic and derived from humans. 

Query 20 can also incorporate criteria that 
compel return of sequence that meets operative requirements 

35 of the subsequent analytical method. Alternatively, or in 
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addition, such operative criteria can be enforced in 
subsequent preprocess step 24 . 

For example, if the function sought to be 
identified is protein coding, query 20 can incorporate 

5 criteria that return from genomic- sequence database 100 
only those sequences present within contigs sufficiently 
long as to have obviated substantial fragmentation of any 
given exon among a plurality of separate sequence 
fragments . 

10 Such criteria can, for example, consist of a 

required minimal individual genomic sequence fragment 
length, such as 10 kb, more typically 20 kb, 30 kb, 40kb, 
and preferably 50 kb or more, as well as an optional, 
further or alternative requirement that sequence from any 

15 given clone, such as a bacterial artificial chromosome 
("BAC"), be presented in no more than a finite maximal 
number of fragments, such as no more than 2 0 separate 
pieces, more typically no more than 15 fragments, even more 
typically no more than about 10 - 12 fragments. 

20 Results using the present invention have shown 

that genomic sequence from bacterial artificial chromosomes 
(BACs) is sufficient for gene prediction analysis according 
to the present invention if the sequence is at least 50 kb 
in length, and if additionally the sequence from any given 

25 BAC is presented in fewer than 15, and preferably fewer 

than 10, fragments.' Accordingly, query 20 can incorporate 
a requirement that data accessioned from BAC sequencing be 
in fewer than 15, preferably fewer than 10, fragments. 

An additional criterion that can be incorporated 

30 into the query can be the date, or range of dates, of 
sequence accession. Although the process has been 
described above as if genomic sequence database 100 were 
static, it is of course understood that the genomic 
sequence databases need not be static, and indeed are 

35 typically updated on a frequent, even hourly, basis. Thus, 
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as further described in Examples 1 and 2, infra, it is 
possible to query the database for newly added sequence, 
either newly added after an absolute date, or newly added 
relative to a prior analysis performed using the methods 
5 and apparatus of the present invention. In this way, the 
process herein described can incorporate a dynamic, 
temporal component . 

One utility of such temporal limitation is to 
identify, from newly accessioned genomic sequence, the 

10 presence of novel genes, particularly those not previously 
identified by EST sequencing {or other sequencing efforts 
that are similarly based upon gene expression) . As further 
described in Example 1, such an approach has shown that 
newly accessioned human genomic sequence, when analyzed for 

15 sequences that function to encode protein, readily 

identifies genes that are novel over those in existing EST 
and other expression databases. This makes the methods of 
the present invention extremely powerful gene discovery 
tools. And as would be appreciated, such gene discovery 

20 can be performed using genomic sequence from species other 
than human. 

If query 20 incorporates multiple criteria, such 
as above-described, the multiple criteria, can be performed 
as a series of separate queries or as a single query, 

25 depending in part upon the query language, the complexity 
of the query, and other considerations well known in the 
database arts. 

If query 20 returns no genomic sequence meeting 
the query criteria, the negative result can be reported by 

30 process 22, and process 200 (and indeed, entire process 10) 
ended 23, as shown. Alternatively, or in addition to 
report and termination of the initial inquiry, a new query 
20 can be generated that takes into account the initial 
negative result. 

35 When query 20 returns sequence meeting the query 
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criteria, the returned sequence is then passed to optional 
preprocessing 24, suitable and specific for the desired 
analytical approach and the particular analytical methods 
thereof to be used in process 25. 
5 Preprocessing 24 can include processes suitable 

for many approaches and methods thereof, as well as 
processes specifically suited for the intended subsequent 
analysis . 

Preprocessing 24 suitable for most approaches and 

10 methods will include elimination of sequence irrelevant to, 
or that would interfere with, the subsequent analysis. 
Such sequence includes repetitive sequence, such as Alu 
repeats and LINE elements,, vector sequence, .artificial 
sequence, such as artificial polylinkers, and the like. 

15 Such removal can readily be performed by identification and 
subsequent masking of the undesired sequence. 

Identification can be effected by comparing the 
genomic sequence returned by query 20 with public or 
private databases containing known repetitive sequence, 

20 vector sequence, artificial sequence, and other artif actual 
sequence. Such comparison can readily be done using 
programs well known in the art, such as CROSS_MATCH, or by 
proprietary sequence comparison programs the engineering of 
which is well within the skill in the art. 

25 Alternatively, or in addition, undesirable, 

including artif actual, sequence can be identified 
algorithmically without comparison to external databases 
and thereafter removed. For example, synthetic polyl inker 
sequence can be identified by an algorithm that identifies 

30 a significantly higher than average density of known 

restriction sites. As another example, vector sequence can 
be identified by algorithms that identify nucleotide or 
codon usage at variance with that of the bulk of the 
genomic sequence. 

35 Once identified, undesired sequence can be 
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removed. Removal can usefully be done by masking the 
undesired sequence as, for example, by converting the 
specific nucleotide references to one that is unrecognized 
by the subsequent bioinf ormatic algorithms, such as "X" . 

5 Alternatively, but at present less preferred, the undesired 
sequence can be excised from the returned genomic sequence, 
leaving gaps. 

Preprocessing 24 can further include selection 
from among duplicative sequences of that one sequence of 

10 highest quality. Higher quality can be measured as a lower 
percentage of, fewest number of, or least densely clustered 
occurrence of ambiguous nucleotides, defined as those 
nucleotides that are identified in the genomic sequence 
using symbols indicating ambiguity. Higher quality can 

15 also or alternatively be valued by presence in the longest 
contig. 

Preprocessing 24 can, and often will, also 
include formatting of the data as specifically appropriate 
for passage to the analytical algorithms of process 25. 

20 Such formatting can and typically will include, inter alia,, 
addition of a unique sequence identifier, either derived 
from the original accession number in genomic sequence 
database 100, or newly applied, and can further include 
additional annotation. Formatting can include conversion 

25 from one to another sequence listing standard, such as 

conversion to or from FASTA or the like, depending upon the 
input expected by the subsequent process. 

Preprocessing, which can be optional depending 
upon the function desired to be identified and the 

30 informational requirements of the methods for effecting 

such identification, is followed by sequence processing 25, 
where sequences with the desired function are identified 
within the genomic sequence. 

As mentioned above, such functions can include, ' 

35 but are not limited to, encoding protein, regulating 
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transcription, regulating message transport after 
transcription into mRNA, regulating message splicing after 
transcription, of regulating message degradation, and the 
like. Other functions include directing somatic 
5 recombination events, contributing to chromosomal stability 
or movement, contributing to allelic exclusion or X 
chromosome inactivation, or the like. 

The methods of the present invention are 
particularly useful for gene discovery, that is, for 

10 identifying, from genomic sequence, regions that function 
to encode genes, and in a particularly useful embodiment, 
for identifying regions that function to encode genes not 
hitherto identified by expression-based or directed cloning 
and sequencing. In conjunction with verification using the 

15 novel single exon microarrays of the present invention, as 
further described below, the methods herein described 
become powerful gene discovery tools . 

Accordingly, in a preferred embodiment of the 
present invention, process 25 is used to identify putative 

20 coding regions. Two preferred approaches in process 25 for 
identifying sequence that encodes putative genes are gene 
prediction and comparative sequence analysis. 

Gene prediction can be performed using any of a 
number of algorithmic methods, embodied in one or more 

25 software programs, that identify open reading frames (ORFs) 
using a variety of heuristics, such as GRAIL, DICTION, and 
GENEFINDER. Comparative sequence analysis similarly can be 
performed using any of a .variety of known programs that 
identify regions with lower sequence variability. 

30 As further described in Example 1, below, gene 

finding software programs yield a range of results. For 
the newly accessioned human genomic sequence input in 
Example 1, for example, GRAIL identified the greatest 
percentage of genomic sequence as putative coding region, 

35 2% of the data analyzed; GENEFINDER was second, calling 1%; 
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and DICTION yielded the least putative coding region, with 
0.8% of genomic sequence called as coding region. 

Increased reliability can be obtained when 
consensus is required among several such methods. Although 
5 discussed herein particularly with respect to exon calling, 
consensus among methods will in general increase 
reliability of predicting other functions as well. 

Thus, as indicated by query 26, sequence 
processing 25, optionally with preprocessing 24, can be 

10 repeated with a different method, with consensus among such 
iterations determined and reported in process 27. 

Process 27 compares the several outputs for a 
given input genomic sequence and identifies consensus among 
the separately reported results. The consensus itself, as 

15 well as the sequence meeting that consensus, is then stored 
in process 29a, displayed in process 2 9b, and/or output to 
process 300 for subsequent identification of a subset 
thereof suitable for assay. 

Multiple levels of consensus can be calculated 

20 and reported by process 27. For example, as further 
described in Example 1, infra, process 27 can report 
consensus as between all specific pairs of methods of gene 
prediction, as consensus among any one or more of the pairs 
of methods of gene prediction, or as among all of the gene 

25 prediction algorithms used. Thus, in Example 1., process 27 
reported that GRAIL and GENEFINDER programs agreed on 0.7% 
of genomic sequence, that GRAIL and DICTION agreed on 0.5% 
of genomic sequence, and that the three programs together 
agreed on 0.25% of the data analyzed. Put another way, 

30 0.25% of the genomic sequence was identified by all three 
of the programs as containing putative coding region. 

Furthermore, consensus can be required among 
different approaches to identifying a chosen function. 

For example, if the function desired to be 

35 identified is coding of protein sequence, and a first used 
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approach to exon calling is gene prediction, the process 
can be repeated on the same input sequence, or subset 
thereof, with another approach, such as comparative 
sequence analysis. In such a case, where comparative 

5 sequence analysis follows gene prediction, the comparison 
can be performed not only on genomic nucleic acid sequence, 
but additionally or alternatively can be performed on the 
predicted amino acid sequence translated from the ORFs 
prior identified by the gene prediction approach. 
10 Although shown as an iterative process, the 

multiple analyses required to achieve consensus can be done 
in series, in parallel, or some combination thereof. 

Predicted functional sequence, optionally 
representing a consensus among a plurality of methods and 
15 approaches for determination thereof, is passed to process 
300 for identification of a subset thereof for functional 
assay . 

In the preferred embodiment of the methods of the 
present invention, wherein the function sought to be 

20 identified is protein coding, process 300 is used to 
identify a subset thereof suitable for experimental 
verification by physical and/or bioinformatic approaches. 

For example, putative ORFs identified in process 
200 can be classified, or binned, bioinf ormatically into 

25 putative genes. This binning can be based inter alia upon 
consideration of the average number of exons/gene in the 
species chosen for analysis, upon density of exons that 
have been called on the genomic sequence, and other 
empirical rules. Thereafter, one or more among the gene- 

30 specific ORFs can be chosen for subsequent use in gene 
expression assay. 

Where such subsequent gene expression assay uses 
amplified nucleic acid, considerations such as desired 
amplicon length, primer synthesis requirements, putative 

35 exon length, sequence GC content, existence of possible 

28 
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secondary structure, and the like can be used to identify 
and select those ORFs that appear most likely successfully 
to amplify. Where subsequent gene expression assay relies 
upon nucleic acid hybridization, whether or not using 
5 amplified product, further considerations involving 

hybridization stringency can be applied to identify that 
subset of sequences that will most readily permit sequence- 
specific discrimination at a chosen hybridization and wash 
stringency. One particular such consideration is avoidance 

10 of putative exons that span repetitive sequence; such 

sequence can hybridize spuriously to nonspecific message, 
reducing specific signal in the hybridization. 

For bioinf ormatic assay, there are fewer 
constraints on the sequences that can be tested 

15 experimentally, and in this latter case therefore process 
3 00 can output the entirety of the input sequence. 

The subset of sequences identified by process 300 
as suitable for use in assay is then used in process 400 to 
create the physical and/or informational substrate for 

20 experimental verification of the predictions made in 
process 200, and thereafter to assay those substrates. 

As mentioned, the methods of the present 
invention are particularly useful for identifying potential 
coding regions within genomic sequence. In a preferred 

25 embodiment of process 400, therefore, the expression of the 
sequences predicted to encode protein is verified. The 
combination of the predictive and experimental methods 
provides a powerful gene discovery engine. 

Thus, in another aspect, the present invention 

30 provides methods and apparatus for verifying the expression 
of putative genes identified within genomic sequence. In 
particular, the invention provides a novel method of 
verifying gene expression in which expression of predicted 
ORFs is measured and confirmed using a novel type of 

35 nucleic acid microarray, the genome-derived single exon 
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nucleic acid microarrays of the present invention. 

Putative ORFs as predicted by a consensus of gene 
calling, particularly gene prediction, algorithms in 
process 200, and as further identified as suitable by 
5 process 300, are amplified from genomic DNA using the 
polymerase chain reaction (PCR) . Although PGR is 
conveniently used, other amplification approaches can also 
be used. 

Amplification schemes can be designed to capture 
10 the entirety of each predicted ORF in an amplicon with 
minimal additional (that is, intronic or intergenic) 
sequence. Because ORFs predicted from human genomic 
sequence using the methods of the present invention differ 
in length, such an approach results in amplicons of varying 
15 length. 

However, most predicted ORFs are shorter than 500 
bp in length, and although amplicons of at least about 100 
or 200 base pairs can be immobilized as probes on nucleic 
acid microarrays, early experimental results using the 

20 methods of the present invention have suggested that longer 
amplicons, at least about 400 or 500 base pairs, are more 
effective. Furthermore, certain advantages derive from 
application to the microarray of amplicons of defined size. 
Therefore, amplification schemes can 

25 alternatively, and preferably, be designed to amplify 

regions of defined size, preferably at least about 300, 400 
or 500 bp, centered about each predicted ORF. Such an 
approach results in a population of amplicons of limited 
size diversity, but that typically contain intronic and/or 

30 intergenic nucleic acid in addition to putative ORF. 

Conversely, somewhat fewer than 10% of ORFs 
predicted from human genomic sequence according to the 
methods of the present invention exceed 500 bp in length. 
Portions of such extended ORFs, preferably at least about 

35 300,400 or 500 bp in length, can be amplified. However, it 
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has been discovered that the percentage success at 
amplifying pieces of such ORFs is low, and that such 
putative exons are more effectively amplified when larger 
fragments, at least about 1000 or 1500 bp, and even as 
5 large as 2000 bp are amplified. 

The putative ORFs selected in process 300 are 
thus input into one or more primer design programs, such as 
PRIMER3 {available online for use at 

http://www-genome.wi.mit.edu/cgi-bin/primer/ ), with a goal 

10 of amplifying at least about 500 base pairs of genomic 

sequence centered within or about ORFs predicted to be no 
more than about 500 bp, or at least about 1000 - 1500 bp of 
genomic sequence for ORFs predicted to exceed 500 bp in 
length, and the primers synthesized by standard techniques. 

15 Primers with the requisite sequences can be purchased 
commercially or synthesized by standard techniques. 

Conveniently, a first predetermined sequence can 
be added commonly to the ORF- specific 5' primer and a 
second, typically different, predetermined sequence 

20 commonly added to each 3" ORF-unique primer. This serves 
to immortalize the amplicon, that is, serves to permit 
further amplification of any amplicon using a single set of 
primers complementary respectively to the common 5' and 
common 3' sequence elements. The presence of these 

25 "universal" priming sequences further facilitates later 
sequence verification, providing a sequence common to all 
amplicons at which to prime sequencing reactions. The 
common 5 ' and 3 ' sequences further serve to add a cloning 
site should any of the ORFs warrant further study. 

30 Such predetermined sequence is usefully at least 

about 10, 12 or 15 nt in length, and usually does not 
exceed about 25 nt in length. The "universal" priming \^ 
sequences used in the examples presented infra were each .16 
nt long. 

35 The genomic DNA to be used as substrate for 
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amplification will come from the eukaryotic species from 
which the genomic sequence data had originally been 
obtained, or a closely related species, and can 
conveniently be prepared by well known techniques from 
5 somatic or germline tissue or cultured cells of the 

organism. See, e.g., Short Protocols in Molecular Biology 
: A Compendium of Methods from Current Protocols in 
Molecular Biology , Ausubel et al . (eds.), 4*^^ edition 
(April 1999) , John Wiley & Sons {ISBN: 047132938X) and 

10 Maniatis et al . , Molecular Cloning : A Laboratory Manual , 
2""^ edition (December 1989), Cold Spring Harbor Laboratory 
Press (ISBN: 0879693096) . Many such prepared genomic DNAs 
are available commercially, with th^ human genomic DNAs 
additionally having certification of donor informed 

15 consent. 

Although the intronic and intergenic material 
flanking putative coding regions in the ampl icons could 
potentially interfere with hybridizations during microarray 
experiments, we have found, surprisingly, that differential 

20 expression ratios are not significantly affected. Rather, 
the predominant effect of exon size is to alter the 
absolute signal intensity, rather than its ratio. Equally 
surprising, the art had suggested that single exon probes 
would not provide sufficient signal intensity for high 

25 stringency hybridization analyses; we find that such probes 
not only provide adequate signal, but have substantial 
advantages, as herein described. 

After partial purification, as by size exclusion 
spin column, with or without confirmation as to amplicon 

30 quality as by gel electrophoresis, each amplicon (single 
exon probe) is disposed in an array upon a support 
substrate . 

Methods for creating microarrays by deposition 
and fixation of nucleic acids onto support substrates are 
35 well known in the art (Reviewed by Schena et al., see 
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above) . 

Typically, the support substrate will be glass, 
although other materials, such as amorphous or crystalline 
silicon or plastics. Such plastics include 
5 polymethylacrylic, polyethylene, polypropylene, 

polyacrylate, polymethylmethacrylate, polyvinylchloride , 
polytetraf luoroethylene, polystyrene, polycarbonate, 
polyacetal, polysulfone, celluloseacetate , 

cellulosenitrate, nitrocellulose, or mixtures thereof, can 
10 also be used. Typically, the support will be rectangular, 
although other shapes, particularly circular disks and even 
spheres, present certain advantages. Particularly 
advantageous alternatives to glass slides as support 
substrates for array of nucleic acids are optical discs, as 
15 described in WO 98/12559. 

The amplified nucleic acids can be attached 
covalently to a surface of the support substrate or, more 
typically, applied to a derivatized surface in a chaotropic 
agent that facilitates denaturation and adherence by 
20 presumed noncovalent interactions, or some combination 
thereof . 

Robotic spotting devices useful for arraying 
nucleic acids on support substrates can be constructed 
using public domain specifications {The MGuide, version 

25 2.0, http://cmgm.stanford.edu/pbrown/mguide/index.html), or 
can conveniently be purchased from commercial sources 
(MicroArray Genii Spotter and MicroArray Genlll Spotter, 
Molecular Dynamics, Inc., Sunnyvale, CA) . Spotting can 
also be effected by printing .methods , including those using 

30 ink jet technology. 

As is well known. in the art, microarrays 
typically also contain immobilized control nucleic acids. 
For controls useful in providing measurements of background 
signal for the genome -derived single exon microarrays of 

35 the present invention, a plurality of E. coli genes can 
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readily be used. As further described in Example 1, 16 or 
32 E. coli genes suffice to provide a robust measure of 
background noise in such microarrays. 

As is well known in the art, the amplified 
5 product disposed in arrays on a support substrate to create 
a nucleic acid microarray can consist entirely of natural 
nucleotides linked by phosphodiester bonds, or 
alternatively can include either nonnative nucleotides, 
alternative internucleotide linkages, or both, so long as 

10 complementary binding can be obtained in the hybridization. 
If enzymatic amplification is used to produce the 
immobilized probes, the amplifying enzyme will impose 
certain further constraints upon the types of nucleic acid 
analogs that can be generated. 

15 Although particularly described herein as using 

high density microarrays constructed on planar substrates, 
the methods of the present invention for confirming the 
expression of ORFs predicted from genomic sequence can use 
any of the known types of microarrays, as herein defined, 

20 including lower density planar arrays, and microarrays on 
nonplanar, nonunitary, distributed substrates. 

For example, gene expression can be confirmed 
using hybridization to lower density arrays, such as those 
constructed on membranes, such as nitrocellulose, nylon, 

25 and positively-charged derivatized nylon membranes. 
Further, gene expression can also be confirmed using 
nonplanar, bead-based microarrays such as are described in 
Brenner et al , , Proc, Natl. Acad. Sci, USA 97 {4 ): 166501670 
(2000); U.S. Patent No. 6,057,107; and U.S. Patent No. 

30 5,736,330. In theory, a packed collection of such beads 
provides in aggregate a higher density of nucleic acid 
probe than can be achieved with spotting or lithography 
techniques on a single planar substrate. 

Planar microarrays on solid substrates, however, 

35 provide certain useful advantages, including high 
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throughput and compatibility with existing readers. For 
example, each standard microscope slide can include at 
least 1000, typically at least 2000, preferably 5000 and 
upto 10,000 - 50,000 or more nucleic acid probes of 
5 discrete sequence. The number of sequences deposited will 
depend on their required application. 

Each putative gene can be represented in the 
array by a single predicted ORF. Alternatively, genes can 
be represented by more than one predicted ORF. For 

10 purposes of measuring differential splicing, more than one 
predicted ORF will be provided for a putative gene. And as 
is well known in the art, each probe of defined sequence, 
representing a single predicted ORF, can be deposited in a 
plurality of locations on a single microarray to provide 

15 redundancy of signal. 

The genome-derived single exon microarrays 
described above differ in several fundamental and 
advantageous ways from microarrays presently used in the 
gene expression art, including (1) those created by 

20 deposition of mRNA-derived nucleic acids, (2). those created 
by in situ synthesis of oligonucleotide probes, and (3) 
those constructed from yeast genomic DNA. 

Most nucleic acid microarrays that are in use for 
study of eukaryotic gene expression have as immobilized 

25 probes nucleic acids that are derived - either directly or 
indirectly — from expressed message. As discussed above, 
it is common, for example, for such microarrays to be 
derived from cDNA/EST libraries, either from those 
previously described in the literature, see Lennon et al . , 

30 or from the de novo construction of "problem specific" 
libraries targeted at a particular biological question, 
R.S. Thomas et al . , Cancer Res. (in press) . Such 
microarrays are herein collectively denominated "EST 
microarrays" . 

35 Such EST microarrays by definition can measure 
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expression only of those genes found in EST libraries, 
shown herein to represent only a fraction of expressed 
genes. Furthermore, such libraries - and thus microarrays 
based thereupon - are biased by the tissue or cell type of 
5 message origin, by the expression levels of the respective 
genes within the tissues, and by the ability of the message 
successfully to have been reverse- transcribed and cloned. 

Thus, as further discussed in Example 1, the 
methods of the present invention enable sequences that do 

10 not appear in EST or other expression databases to be 
determined - subsequently arrayed for expression 
measurements could not, therefore, have been represented as 
probes on an EST microarray. And as further demonstrated 
in the examples, infra, the remaining population of genes 

15 identified from genomic sequence by the methods of the 
present invention - that is, the one third of sequences 
that had previously been accessioned in EST or other 
expression databases - are biased toward genes with higher 
expression levels. 

20 Representation of a message in an EST and/or cDNA 

library depends upon the successful reverse transcription, 
optionally but typically with subsequent successful 
cloning, of the message. This introduces substantial bias 
into the population of probes available for arraying in EST 

25 microarrays. 

^ In contrast, neither reverse transcription nor 
cloning is required to produce the probes arrayed on the 
genome-derived single exon microarrays of the present 
invention. And although the ultimate deposition of a probe 

30 on the genome -derived single exon microarray of the present 
invention depends upon a successful amplification from 
genomic material, a priori knowledge of the sequence of the 
desired amplicon affords greater opportunity to recover any 
given probe sequence recalcitrant to amplification than is 

35 afforded by the requirement for successful reverse 
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transcription and cloning of unknown message in EST 
approaches . 

Thus, the genome -derived single exon microarrays 
of the present invention present a far greater diversity of 

5 probes for measuring gene expression, with far less bias, 
than do EST microarrays presently used in the art. 

As a further consequence of their ultimate origin 
from expressed message, the probes in EST microarrays often 
contain poly-A {or complementary poly-T) stretches derived 

10 from the poly-A tail of mature mRNA. These homopolymeric 
stretches contribute to cross -hybridization, that is, to a 
spurious signal occasioned by hybridization to the 
homopolymeric tail of a labeled cDNA that lacks sequence 
homology to the gene-specific portion of the probe. 

15 In contrast, the probes arrayed in the genome - 

derived single exon microarrays of the present invention 
lack homopolymeric stretches derived from message 
polyadenylation, and thus can provide more specific signal. 
Typically, at least about 50, 60 or 75% of the probes on 

20 the genome-derived single exon microarrays of the present 
invention lack homopolymeric regions consisting of A or T, 
where a homopolymeric region is defined for purposes herein 
as stretches of 25 or more, typically 30 or more, identical 
nucleotides . 

25 A further distinction, which also affects the 

specificity of hybridization, is occasioned by the typical 
derivation of EST microarray probes from cloned material. 
Because much of the probe material disposed as probes on 
EST microarrays is excised or amplified from plasmid, 

30 phage, or phagemid vectors, EST microarrays typically 

include a fair amount of vector sequence, more so when the 
probes are amplified, rather than excised, from the vector. 

In contrast, the vast majority of probes in the 
genome-derived single exon microarrays of the present 

35 invention contain no prokaryotic or bacteriophage vector 
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sequence, having been amplified directly or indirectly from 
genomic DNA. Typically, therefore, at least about 50, 60, 
70 or 80% or more of individual exon- including probes 
disposed on a genome -derived single exon microarray of the 
5 present invention lack vector sequence, and particularly 
lack sequences drawn from plasmids and bacteriophage. 
Preferably, at least about 85, 90 or more than 90% of exon- 
including probes in the genome -derived single exon 
microarray of the present invention lack vector sequence. 

10 With attention to removal of vector sequences through 

preprocessing 24, percentages of vector- free exon-including 
probes can be as high as 95 - 99%. The substantial absence 
of vector sequence from the genome-derived single exon 
microarrays of the present invention results in greater 

15 specificity during hybridization, since spurious cross- 
hybridization to a probe vector sequence is reduced. 

As a further consequence of excision or 
amplification of probes from vectors in construction of EST 
microarrays, the probes arrayed thereon often contain 

20 artificial sequence, derived from vector polylinker 

multiple cloning sites, at both 5' and 3' ends. The probes 
disposed upon the genome-derived single exon microarrays 
need have no such artificial sequence appended thereto. 

As mentioned above, however, the ORF-specific 

25 primers used to amplify putative ORFs can include 

artificial sequences, typically 5* to the ORF-specific 
primer sequence, useful for "universal" (that is, 
independent of ORF sequence) priming of subsequent 
amplification or sequencing reactions, when such 

30 "universal" 5' and/or 3' priming sequences are appended to 
the amplification primers, the probes disposed upon the 
genome -derived single exon microarray will include 
artificial sequence similar to that found in EST 
microarrays. However, the genome -derived single exon 

35 microarray of the present invention can be made without 
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such sequences, and if so constructed, presents an even 
smaller amount of nonspecific sequence that would 
contribute to nonspecific hybridization. 

Yet another consequence of typical use of cloned 
5 material as probes in EST microarrays is that such 
microarrays contain probes that result from cloning 
artifacts, such as chimeric molecules containing coding 
region of two separate genes . Derived from genomic 
material, typically not thereafter cloned, the probes of 

10 the genome -derived single exon microarrays of the present 
invention lack such cloning artifacts, and thus provide 
greater specificity of signal in gene expression 
measurements . 

A further consequence of the cloned origin of 

15 probes on many EST microarrays is that the individual 
probes often have disparate sizes, which can cause the 
optimal hybridization stringency to vary among probes on a 
single microarray. In contrast, as discussed above, the 
probes arrayed on the genome -derived single exon 

20 microarrays of the present invention can readily be 

designed to have a narrow distribution in sizes, with the 
range of probe sizes no greater than about 10% of the 
average size, typically no greater than about 5% of the 
average probe size . 

25 Because of their^ origin from fully- or partially- 

spliced message, probes disposed upon EST arrays will often 
include multiple exons. The percentage of such exon- 
spanning probes in an EST microarray can be calculated, on 
average, based upon the predicted number of exons/gene for 

30 the given species and the average length of the immobilized 
probes. For human genes, the near-complete sequence of 
human chromosome 22, Dunham ,et ai . , Nature 402 { 6761) : 489-95 
(1999), predicts that human genes average 5.5 exons/gene. 
Even with probes of 200 - 500 bp, the vast majority of 

35 human EST microarray probes include more than one exon. 
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In contrast, by virtue of their origin from 
algorithmically identified ORFs in genomic sequence, the 
probes in the genome -derived single exon microarrays of the 
present invention can consist of individual exons . Thus, 
5 in contrast to EST microarrays, at least about 50, 60, 70, 
75, 80, 85, 95 or 99% of probes deposited in the genome- 
derived microarray of the present invention consist of, or 
include, no more than one predicted ORF. 

This provides the ability, not readily achieved 

10 using EST microarrays, to use the genome -derived single 
exon microarrays of the present invention to measure 
tissue -specif ic expression of individual exons, which in 
turn allows differential splicing events to be detected and 
characterized, and in particular, allows the correlation of 

15 differential splicing to tissue- specif ic expression 
patterns. 

Furthermore, the exons that are represented in 
EST microarrays are often biased toward the 3' or 5 ' end of 
their respective genes, since sequencing strategies used 

20 for EST identification are so biased. In contrast, no such 
3' or 5 ' bias necessarily inheres in the selection of exons 
for disposition on the genome -derived single exon 
microarrays of the present invention. 

Conversely, the probes provided on the genome - 

25 derived single exon microarrays of the present invention 
typically, but need not necessarily, include intronic 
and/or intergenic sequence that is absent from EST 
microarrays, which are derived from mature mRNA. 
Typically, at least about 50, 60, 70, 80 or 90% of the 

30 exon- including probes on the genome-derived single exon 

microarrays of the present invention include sequence drawn 
from noncoding regions. As discussed above, the additional 
presence of noncoding region does not significantly 
interfere with measurement of gene expression, and provides 

35 the additional opportunity to assay prespliced RNA, and 
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thus measure such phenomena such as nuclear export control. 

The genome-derived single exon microarrays of the 
present invention are also quite different from in situ 
synthesis microarrays, where probe size is severely 
5 constrained by inadequacies in the photolithographic 
synthesis process. 

Typically, probes arrayed on in situ synthesis 
microarrays are limited to a maximum of about 25 bp. As a 
well known consequence, hybridization to such chips must be 

10 performed at low stringency. In order, therefore, to 
achieve unambiguous sequence - specif ic hybridization 
results, the in situ synthesis microarray requires 
substantial redundancy, with concomitant programmed 
arraying for each probe of probe analogues with altered 

15 {i.e., mismatched) sequence. 

In contrast, the longer probe length of the 
genome-derived single exon microarrays of the present 
invention allows much higher stringency hybridization and 
wash. Typically, therefore, exon- including probes on the 

20 genome -derived single exon microarrays of the present 
invention average at least about 100, 200, 300, 400 or 
500 bp in length. By obviating the need for substantial 
probe redundancy, this approach permits a higher density of 
probes for discrete exons or genes to be arrayed on the 

25 microarrays of the present invention than can be achieved 
for in situ synthesis microarrays. 

A further distinction is that the probes in in 
situ synthesis microarrays typically are covalently linked 
to the substrate surface. In contrast, the probes disposed 

30 on the genome -derived microarray of the present invention 
typically are-, but need not necessarily be, bound 
noncovalently to the substrate. 

Furthermore, the short probe size on in situ 
microarrays causes large percentage differences in the 

35 melting temperature of probes hybridized to their 
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complementary target sequence, and thus causes large 
percentage differences in the theoretically optimum 
stringency across the array as a whole. 

In contrast, the larger probe size in the 
5 microarrays of the present invention create lower 

percentage differences in melting temperature across the 
range of arrayed probes . 

A further significant advantage of the 
microarrays of the present invention over in situ 
10 synthesized arrays is that the quality of each individual 
probe can be confirmed before deposition. In contrast, the 
quality of probes cannot be assessed on a probe -by-probe 
basis for the in situ synthesized microarrays presently 
being used. 

15 The genome -derived single exon microarrays of the 

present invention are also distinguished over, and present 
substantial benefits over, the genome -derived microarrays 
from lower eukaryotes such as yeast. Lashkari et al . , 
Proc. Natl. Acad. Sci . USA 94:13057-13062 (1997). 

20 Only about 220 - 250 of the 6100 or so nuclear 

genes in Saccharowyces cerevisiae - that is, only about 4 
- 5% — have standard, spliceosomal, introns, Lopez et al . , 
Nucl . Acids Res, 28:85-86 {2 000); Spingola et al . , RNA 
5(2):221-34 {1999). Furthermore, the entire yeast genome 

25 has already been sequenced. These two facts permit the 

ready amplification and disposition of single-ORF amplicons 
on such microarray without the requirement for antecedent 
use of gene prediction and/or comparative sequence 
analyses . 

30 Thus, a significant aspect of the present 

invention is the ability to identify and to confirm 
expression of predicted coding regions in genomic sequence 
drawn from eukaryotic organisms that have a higher 
percentage of genes having introns than do yeast such as 

35 Saccharoiuyces cerevisiae, particularly in genomic sequence 
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drawn from eukaryotes in which at least about 10, 20 or 50% 
of protein-encoding genes have introns. In preferred 
embodiments, the methods and apparatus of the present 
invention are used to identify and confirm expression of 

5 novel genes from genomic sequence of eukaryotes in which 
the average number of introns per gene is at least about 
one, two or three or more. 

After the physical substrate is prepared, 
experimental verification of predicted function is 
10 performed. 

In a preferred embodiment of the present 
invention, where the function sought to be identified in 
genomic sequence is protein coding, experimental 
verification is performed by measuring expression of the 

15 putative ORFs, typically through nucleic acid hybridization 
experiments, and in particularly preferred embodiments, 
through hybridization to genome -derived single exon 
microarrays prepared as above- described. 

Expression is conveniently measured and expressed 

20 for each probe in the microarray as a ratio of the 

expression measured concurrently in a plurality of mRNA 
sources, according to techniques well known in the 
microarray art, Reviewed in Schena et al . , and as further 
described in Example 2, below. The mRNA source for the 

25 reference against which specific expression is measured can 
be drawn from a homogeneous mRNA source, such as a single 
cultured cell-type, or alternatively can be heterogeneous, 
as from a pool of mRNA derived from multiple tissues and/or 
cell types, as further described in Example 2, infra. 

30 mRNA can be prepared by standard techniques, see 

Ausubel et al. and Maniatis et al . , or purchased 
commercially. The mRNA is then typically reverse- 
transcribed in the presence of labeled nucleotides: the 
index source (that in which expression is desired to be 

35 measured) is reverse transcribed in the presence of 
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nucleotides labeled with a first label, typically a 
fluorophore ( f luorochrome ; Eluor; fluorescent dye) the 
reference source is reverse " transcribed in the presence of 
a second label, typically a . fluorophore , typically 
5 f luorometrically-distinguishable from the first label. As 
further described in Example 2, infra. Cy3 and Cy5 dyes 
prove particularly useful in these methods. After partial 
purification of the index and reference targets, 
hybridization to the probe array is conducted according to 

10 standard techniques, typically under a coverslip. 

After wash, microarrays are conveniently scanned 
using a commercial microarray scanning device, such as a 
Gen3 Scanner (Molecular Dynamics^ Sunnyvale, _CA) . Data on _ 
expression is then passed, with or without interim storage, 

15 to process 500, where the results for each probe are 
related to the original sequence. 

Often, hybridization of target material to the 
genome -derived single exon microarray will identify certain 
of the probes thereon as of particular interest. Thus, it 

20 is often desirable that the user be able readily to obtain 
sufficient quantities of an individual probe, either for 
subsequent arrayed deposition upon an additional support 
substrate, often as part of a microarray having a plurality 
of probes so identified, or alternatively or additionally 

25 as a solitary solid-phase or solution-phase probe, for 
further use. 

Thus, in another aspect, the present invention 
provides compositions and kits for the ready production of 
nucleic acids identical in sequence to, or substantially 
30 identical in sequence to, probes on the genome-derived 
single exon microarrays of the present invention. 

In this aspect, a small quantity of each probe is 
disposed, typically without attachment to substrate, in a 
spatially-addressable ordered set, typically one per well 
35 of a microtiter dish. Although a 96 well microtiter plate 
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can be used, greater efficiency is obtained using higher 
density arrays, such as are provided by microtiter plates 
having 384, 864, 1536, 3456, 6144, or 9600 wells, and 
although microtiter plates having physical depressions 
5 (wells) are conveniently used, any 'device that permits 
addressable withdrawal of reagent from fluidly- 
noncommunicating areas can be used. 

In this aspect of the invention, therefore, a 
fluidly noncommunicating addressable ordered set of 

10 individual probes, corresponding to those on a genome- 
derived single exon microarray, is provided, with each 
probe in sufficient quantity to permit amplification, such 
as by PCR- As earlier mentioned, the ORF-specific 
5' primers used for genomic amplification can have a first 

15 common sequence added thereto, and the ORF-specific 3' 
primers used for genomic amplification can have a second, 
different, common sequence added thereto, thus permitting, 
in this preferred embodiment, the use of a single set of 5' 
and 3" primers to amplify any one of the probes from the 

20 amplifiable ordered set. . 

Each discrete amplifiable probe can also be 
packaged with amplification primers, solutes, buffers, 
etc., and can be provided in dry (e.g., lyophilized) form 
or wet, in the latter case typically with addition of 

25 agents that retard evaporation. 

In another aspect of the present invention, a 
genome -derived single-exon microarray is packaged together 
with such an ordered set of amplifiable probes 
corresponding to the probes, or one or more subsets of 

30 probes, thereon. In alternative embodiments, the ordered 
set of amplifiable probes is packaged separately from the 
genome -derived single exon microarray. 

In some embodiments, the microarray and/or 
ordered probe set are further packaged with recordable 

35 media that provide probe identification and addressing 
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information, and that can additionally contain annotation 
information, such as gene expression data. Such recordable 
media can be packaged with the microarray, with the ordered 
probe set, or with both. 

5 If the microarray is constructed on a substrate 

that incorporates recordable media, such as is described in 
international patent application no. WO 98/12559, then 
separate packaging of the genome-derived single exon 
microarray and the bioinf ormatic information is not 
10 required. 

The amount of amplifiable probe material should 
be sufficient to permit at least one amplification 
sufficient for subsequent hybridization assay. 

Although the use of high density genome-derived 

15 microarrays on solid planar substrates is presently a 
preferred approach for the physical confirmation and 
characterization of the expression of sequences predicted 
to encode protein, other types of microarrays (as herein 
defined) can also be used. 

20 Furthermore, as earlier mentioned, experimental 

verification of the function predicted from genomic 
sequence in process 200 can be bioinf ormatic , rather than, 
or additional to, physical verification. 

For example, where the function desired to be 

25 identified is protein coding, the predicted ORFs can be 

compared bioinf ormatically to sequences known or suspected 
of being expressed. 

Thus, the sequences output from process 300 (or 
process 200), can be used to query expression databases, 

30 such as EST databases, SNP {"single nucleotide 

polymorphism") databases, known cDNA and mRNA sequences, 
SAGE ("serial analysis of gene expression") databases, and 
more generalized sequence databases that allow query for 
expressed sequences. Such query can be done by any 

35 sequence query algorithm, such as BLAST ("basic local 
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alignment search tool") . The results of such query - 
including information on identical sequences and 
information on nonidentical sequences that have diffuse or 
focal regions of sequence homology to the query sequence — 
5 can then be passed directly to process 500, or used to 
inform analyses subsequently undertaken in process 2 00, 
process 300, or process 400. 

Experimental data, whether obtained by physical 
or bioinformatic assay in process 4 00, is passed to process 
10 500 where it is usefully related to the sequence data 

itself, a process colloquially termed "annotation". Such 
annotation can be done using any technique that usefully 
relates the functional information to the sequence, as, for 
example, by incorporating the functional data into the 
15 record itself, by linking records in a hierarchical or 

relational database, by linking to external databases, or 
by a combination thereof. Such database techniques are 
well within the skill in the art. 

The annotated sequence data can be stored 
20 locally, uploaded to genomic sequence database 100, and/or 
displayed 800. 

The methods and apparatus of the present 
invention rapidly produce functional information from 
genomic sequence. Coupled with the escalating pace at 
25 which sequence now accumulates, the rapid pace of sequence 
annotation produces a need for methods of displaying the 
information in meaningful ways. 

FIG- 3 shows visual display 80 presenting a 
single genomic sequence annotated according to the present 
30 invention. Because of its nominal resemblance to artistic 
works of Piet Mondrian, visual display 80 is alternatively 
described herein as a "Mondrian" . 

Each of the visual elements of display 8 0 is 
aligned with respect to the genomic sequence being 
35 annotated (hereinafter, the "annotated sequence"). Given 
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the number of nucleotides typically represented *in an 
annotated sequence, representation of individual 
nucleotides would rarely be readable in hard copy output of 
display 80. Typically, therefore, the annotated sequence 

5 is schematized as rectangle 89, extending from the left 
border of display 80 to its right border. By convention 
herein, the left border of rectangle 89 represents the 
first nucleotide of the sequence and the right border of 
rectangle 89 represents the last nucleotide of the 
10 sequence. 

As further discussed below, however, the Mondrian 
visual display of annotated sequence can serve as a 
convenient graphical user interface for computerized 
representation, analysis, and query of information stored 

15 electronically. For such use, the individual nucleotides 
can conveniently be linked to the X axis coordinate of 
rectangle 89. This permits the annotated sequence at any 
point within rectangle 89, readily to be viewed, either 
automatically - for example, by time-delayed appearance of 

20 a small overlaid window upon movement of a cursor or other 
pointer over rectangle 89 - or through user intervention, 
as by clicking a mouse or other pointing_ device at a point 
in rectangle 89. 

Visual display 80 is generated after user 

25 specification of the genomic sequence to be displayed. 

Such specification can consist of or include an accession 
number for a single clone (e.g., a single BAG accessioned 
into GenBank) , wherein the starting and stopping 
nucleotides are thus absolutely identified, or 

30 alternatively can consist of or include an anchor or 

fulcrum point about which a chosen range of sequence is 
anchored, thus providing relative endpoints for the 
sequence to be displayed. For example, the user can anchor 
such a range about a given chromosomal map location, gene 

35 name, or even a sequence returned by query for similarity 
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or identity to an input query sequence. When visual 
display 80 is used as a graphical user interface to 
computerized data, additional control over the first and 
last displayed nucleotide will typically be dynamically 

5 selectable, as by use of standard zooming and/or selection 
tools . 

Field 81 of visual display 80 is used to present 
the output from process 200, that is, to present the 
bioinf ormatic prediction of those sequences having the 

10 desired function within the genomic sequence. Functional 
sequences are typically indicated by at least one rectangle 
83 (83a, 83b, 83c) , the left and right borders of which 
respectively indicate, by their X-axis coordinates, the 
starting and ending nucleotides of the region predicted to 

15 have function. 

Where a single bioinf ormatic method or approach 
identifies a plurality of regions having the desired 
function, a plurality of rectangles 83 is disposed 
horizontally in field 81. Where multiple methods and/or 

20 approaches are used to identify function, each such method 
and/or approach can be represented by its own series of 
horizontally disposed rectangles 83, each such horizontally 
disposed series of rectangles offset vertically from those 
representing the results of the other methods and 

25 approaches. 

Thus, rectangles 83a in FIG. 3 represent the 
functional predictions of a first method of a first 
approach for predicting function, rectangles 83b represent 
• the functional predictions of a second method and/or second 

30 approach for predicting that function, and rectangles 83c 
represent the predictions of a third method and/or 
approach. >^ 

where the function desired to be identified is 
protein coding, field 81 is used to present the 

35 • bioinf ormatic prediction of sequences encoding protein. 
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For example, rectangles 83a can represent the results from 
GRAIL or GRAIL II, rectangles 83b can represent the results 
from GENEFINDER, and rectangles 83c can represent the 
results from DICTION. 

i 

5 optionally, and preferably, rectangles 83 

collectively representing predictions of a single method 
and/or approach are identically colored and/or textured, 
and are distinguishable from the color and/or texture used 
for a different method and/or approach. 

10 Alternatively, or^'in addition, the color, hue, 

density, or texture of rectangles 83 can be used further to 
report a measure of the bioinf ormatic reliability of the 
prediction. For example, many gene prediction programs 
will report a measure of the reliability of prediction. 

15 Thus, increasing degrees of such reliability can be 

indicated, e.g., by increasing density of shading. Where 
display 80 is used as a graphical user interface, such 
measures of reliability, and indeed all other results 
output by the program, can additionally or alternatively be 

20 made accessible through linkage from individual rectangles 
83, as by time-delayed window ("tool tip" window), or by 
pointer (e.g., mouse ) -activated link. 

As earlier described, increased predictive 
reliability can be achieved by requiring consensus among 

25 methods and/or approaches to determining function. Thus, 
field 81 can include a horizontal series of rectangles 83 
that indicate one or more degrees of consensus in 
predictions of function. 

Although FIG. 3 shows three series of 

30 horizontally disposed rectangles in field 81, display 80 
can include as few as one such series of rectangles and as 
many as can discriminably be displayed, depending upon the 
number of methods and/or approaches used to predict a given 
function. 

35 Furthermore, field 81 can be used to show 
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predictions of a plurality of different functions. 
However, the increased visual complexity occasioned by such 
display makes more useful the ability of the user to select 
a single function for display. When display 80 is used as 
5 a graphical user interface for computer query and analysis, 
such function can usefully be indicated and user- 
selectable, as by a series of graphical buttons or tabs 
{not shown in FIG. 3) . 

Rectangle 89 is shown in FIG. 3 as including 

10 interposed rectangle 84. Rectangle 84 represents the 
portion of annotated sequence for which predicted 
functional information has been assayed physically, with 
the starting and ending nucleotides of the assayed material 
indicated by the X axis coordinates of the left and right 

15 borders of rectangle 84. Rectangle 85, with optional 
inclusive circles 86 (86a, 86b, and 86c) displays the 
results of such physical assay. 

Although a single rectangle 84 is shown in FIG. 
3, physical assay is not limited to just one region of 

20 annotated genomic sequence. It is expected that an 

increasing percentage of regions predicted to have function 
by process 200 will be assayed physically, and that display 
80 will accordingly, for any given genomic sequence, have 
an increasing number of rectangles 84 and 85, representing 

25 an increased density of sequence annotation. 

Where the function desired to be identified is 
protein coding, rectangle 84 identifies the sequence of the 
probe used to measure expression. In embodiments of the 
present invention where expression is measured using 

30 genome -derived single exon microarrays, rectangle 84 
identifies the sequence included within the probe 
immobilized on the support surface of the microarray. As 
noted supra, such probe will often include a small amount 
of additional, synthetic, material incorporated during 

35 amplification and designed to permit reamplif icat ion of the 
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probe, which sequence is typically not shown in display 80. 

Rectangle 87 is used to present the results of 
bioinformat ic assay of the genomic sequence. For example, 
where the function desired to be identified is protein 
5 coding, process 400 can include bioinformat ic query of 
expression databases with the sequences predicted in 
process 200 to encode exons . And as earlier discussed, 
because bioinf ormatic assay presents fewer constraints than 
does physical assay, often the entire output of process 200 

10 can be used for such assay, without further subsetting 

thereof by process ,300. Therefore, rectangle 87 typically 
need not have separate indicators therein of regions 
submitted for bioinf ormatic assay; that is, rectangle 87 
typically need not have regions therein analogous to 

15 rectangles 84 within rectangle 89. 

Rectangle 87 as shown in FIG. 3 includes smaller 
rectangles 880 and 88. Rectangles 880 indicate regions 
that returned a positive result in the bioinf ormatic assay, 
with rectangles 88 representing regions that did not return 

20 such positive results. Where the function desired to be 
predicted and displayed is protein coding, rectangles 880 
indicate regions of the predicted exons that identify 
sequence with significant similarity in expression 
databases, such as EST, SNP, SAGE databases, with 

25 rectangles 88 indicating genes novel over those identified 
in existing expression data bases. 

Rectangles 880 can further indicate, through 
color, shading, texture, or the like, additional 
. information obtained from bioinf ormatic assay. 

30 For example, where the function assayed and 

displayed is protein coding, the degree of shading of 
rectangles 880 can be used to represent the degree of 
sequence similarity found upon query of expression 
databases. The number of levels of discrimination can be 

35 as few as two (identity, and similarity, where similarity 
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has a user-selectable lower threshold) . Alternatively, as 
many different levels of discrimination can be indicated as 
can visually be discriminated. 

Where display 80 is used as a graphical user 

5 interface, rectangles 880 can additionally provide links 
directly to the sequences identified by the query of 
expression databases, and/or statistical summaries thereof. 
As with each of the precedingly-discussed uses of display 
80 as a graphical user interface, it should be understood 
10 that the information accessed via display 80 need not be 
resident on the computer presenting such display, which 
often will be serving as a client, with the linked 
information resident on one or more remotely located 
servers . 

15 Rectangle 85 displays the results of physical 

assay of the sequence delimited by its left and right 
borders . 

Rectangle 85 can consist of a single rectangle, 
thus indicating a single assay, or alternatively, and 

20 increasingly typically, will consist of a series of 

rectangles {85a, 85b, 85c) indicating separate physical 
assays of the same sequence . 

Where the function assayed is gene expression, 
and where gene expression is assayed as herein described 

25 using simultaneous two-color fluorescent detection of 
hybridization to genome -derived single exon microarrays, 
individual rectangles 85 can be colored to indicate the 
degree of expression relative to control. Conveniently, 
shades of green can be used to depict expression in the 

30 sample over control values, and shades of red used to 

depict expression less than control, corresponding to the 
spectra of the Cy3 and Cy5 dyes conventionally used for 
respective labeling thereof. Additional functional 
information can be provided in the form of circles 86 (86a, 

35 86b, 86c) , where the diameter of the circle can be used to 
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indicate expression intensity. As discussed infra, such 
relative expression {expression ratios) and absolute 
expression {signal intensity) can be expressed using 
normalized values. 
5 Where display 80 is used as a graphical user 

interface, rectangle 85 can be used as a link to further 
information about the assay. For example, where the assay 
is one for gene expression, each rectangle 85 can be used 
to link to information about the source of the hybridized 

10 mRNA, the identity of the control, raw or processed data 
from the microarray scan, or the like. 

FIG. 4 is rendition of display 80 representing 
gene prediction and gene expression for a hypothetical BAC, 
showing conventions used in the Examples presented infra. 

15 BAC sequence {"Chip seq.") 89 is presented, with the 
physically assayed region thereof (corresponding to 
rectangle 84 in FIG. 3) shown in white. Algorithmic gene 
predictions are shown in field 81, with predictions by 
GRAIL shown, predictions by GENEFINDER, and predictions by 

20 DICTION shown. Within rectangle 87, regions of sequence 
that, when used to query expression databases, return 
identical or similar sequences {"EST hit") are shown as 
white rectangles (corresponding to rectangles 880 in FIG. 
3) , gray indicates low homology, and black indicates 

25 unknowns (where black and gray would correspond to 
rectangles 88 in FIG. 3) . 

Although FIGS. 3 and 4 show a single stretch of 
sequence, uninterrupted from left to right, longer 
sequences are usefully represented by vertical stacking of 

30 such individual Mondrians, as shown in FIGS. 9 and 10. 

Single Exon Probes Useful For Measuring Gene Expression 

The methods and apparatus of the present 
35 invention rapidly produce functional information from 
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genomic sequence. Where the function to be identified is 
protein coding, the methods and apparatus of the present 
invention rapidly identify and confirm the expression of 
portions of genomic sequence that function to encode 
5 protein. As a direct result, the methods and apparatus of 
the present invention rapidly yield large numbers of 
single-exon nucleic acid probes, the majority from 
previously unknown genes, each of which is useful for 
measuring and/or surveying expression of a specific gene in 

10 one or more tissues or cell types. 

It is, therefore, another aspect of the present 
invention to provide genome-derived single exon nucleic 
acid probes useful for gene expression analysis, and 
particularly for gene expression analysis by microarray. 

15 Using the methods and genome -derived single-exon 

microarrays of the present invention, we have for example 
readily identified a large number of unique ORFs from human 
genomic sequence. Using single exon probes that encompass 
these ORFs, we have demonstrated, through microarray 

20 hybridization analysis, the expression of 12,673 of these 
ORFs in Fetal liver. 

As would immediately be appreciated by one of 
skill in the art, each single exon probe having 
demonstrable expression in Fetal liver is currently 

25 available for use in measuring the level of its ORF ' s 
expression in Fetal liver. 

Diseases of the liver are a significant cause of 
human morbidity and mortality. Increasingly, genetic 
• factors are being found that contribute to predisposition, 

30 onset, and/or aggressiveness of most, if not all, of these 
diseases; although causative mutations in single genes have 
been identified for- some, these disorders are believed-^f or 
the most part to have polygenic etiologies. 

For example, cirrhosis is a major public health, 

35 problem. In the industrialized world, it is among the top 

55 



WO(H/57277 PCT/US01/(M)669 
ten causes of death; among patients aged 45 to 65, it is 
the third leading cause of death. The high prevalence is 
largely the result of alcohol abuse, but other major 
contributors include chronic^ hepatitis , biliary disease and 
5 iron overload. Approximately* 10-15% are cryptogenic. 

Cirrhosis is a broad description encompassing the 
common end stage of many forms of liver injury. Many 
patients with cirrhosis will remain asymptomatic for years, 
while others show generalized weakness, anorexia, malaise, 

10 and weight loss or, occasionally, more severe symptoms. 

The progression from fibrosis, an early 
consequence of liver disease, to cirrhosis, and the 
specific histologic morphology that characterizes cirrhosis 
depend on the extent of injury, the presence of continuing 

15 damage, and the response of the liver to damage. The liver 
may be injured acutely and severely {e.g. necrosis with 
hepatitis), moderately over, months or years {e.g. biliary 
tract obstruction and chronic active hepatitis) , or 
modestly but continuously (e.g. alcohol abuse) . 

20 During the repair process, new vessels 

connecting the hepatic artery and portal vein to the 
hepatic venules form within the fibrous sheath that 
surrounds the surviving nodules of liver cells. These 
vessels restore the intrahepatic circulatory pathway, but 

25 provide relatively low- volume, high-pressure drainage that 
is less efficient than normal and results in increased 
portal vein pressure {portal hypertension) . Thus, 
cirrhosis is not static and its features depend on the 
disease activity and stage. 

30 As cirrhosis is the end stage of many forms of 

liver disease, many genes have been identified that can 
contribute to the development of cirrhosis. These include, 
e.g., the genes responsible for Wilson disease (Online 
Mendelian Inheritance of Man {"OMIM") 277900), type IV 

35 glycogen storage disease (OMIM 232500) , galactosemia {OMIM 
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230400), and a deficiency of alpha- 1-antitrypsin (OMIM 
107400) . There is substantial evidence, however, for as 
yet uncharacterized loci which cause cirrhosis. 

For example, Iber and Maddrey, Prog. Liver Dis. 
5 2; 290-302 (1965), reviewed 13 previously reported families 
and 8 new to this study, each with 2 or more affected 
members. They pointed out that, with a single exception, 
the multiple cases were in the same generation. Within a 
given family, the age of onset, clinical course, and biopsy 
10 findings were very similar, but there were wide differences 
between families. 

Kalra et al . , Hum. Hered. 32:170-175 (1982) 
studied the families of 220 cases of Indian childhood 
cirrhosis and 70 families of age-matched controls. The 
15 hypotheses of autosomal recessive, partial sex- linkage, and 
doubly recessive inheritance were found untenable and the 
authors concluded that multifactorial inheritance was most 
plausible. Lefkowitch et al . , New Eng. J. Med. 307:271-277 
(1982) described 4 white American sibs who died between 
20 ages 4.5 and 6 years of cirrhosis that closely resembled 
that of the childhood cirrhosis of Asiatic Indians. 

Another example of uncharacterized loci which 
cause cirrhosis are those related to the risk of 
alcoholism. 

25 Cloninger, Science 236:410-416 (1987), defined 

two separate types of alcoholism. According to these 
definitions, type 1 alcohol abuse has its usual onset after 
the age of 25 years and is characterized by severe 
psychological dependence and guilt. Type 1 occurs in both 

30 men and women and requires both genetic and environmental 
factors to become manifest. By contrast, type 2 alcohol 
abuse has its onset before the age of 25,- persons with this 
type of alcoholism are characterized by their inability to 
abstain from alcohol and by frequent aggressive and 

35 antisocial behavior. Type 2 alcoholism is rarely found in 
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women and is much more heritable. 

Despite considerable effort to identify genes 
related to the risk of alcoholism, relatively few genes 
have been identified. Some of this work has suggested a 

5 relationship between the metabolism of dopamine and 

alcoholism. Blum et al . , J.A.M.A. 263:2055-2060 (1990) and 
Bolos et al., J. A.M. A. 264:3156-3160 (1990) investigated 
the relationship of the dopamine D2 receptor {DRD2; OMIM 
126450) to alcoholism, but the sample size was small and 

10 their. results were inconclusive. However, Tiihonen et al . , 
Molec. Psychiat. 4, 286-289 (1999), found a markedly higher 
frequency in a population of type 1 alcoholics of the low 
activity allele of the enzyme catechol -0-methyltransf erase 
(COMT, OMIM 116790) , which has a crucial role in the 

15 metabolism of dopamine, suggesting a role for dopamine 
metabolism in increased risk of alcoholism. For a brief 
review of recent progress toward the identification of 
genes related to risk for alcoholism see Buck, Genome 
9:927-928 (1998) . 

20 As another example, multiple genes have been 

shown to predispose to hyperlipoproteinemia or 
hyperlipidemia. Much attention has been focused on these 
disorders because there is a strong association of 
hyperlipidemia, especially hypercholesterolemia, with 

25 development of coronary artery disease. Coronary artery 
disease accounts for at least 25% of all deaths in the 
United States. Coronary .artery disease results when the 
arteries supplying the heart muscle become occluded by 
plaques composed of lipids like cholesterol, blood clotting 

30 components and blood cells. 

The major plasma lipids circulate bound to 
proteins as macromolecular complexes called lipoproteins. 
Although closely interrelated, the major lipoprotein 
classes - chylomicron, very- low-density lipoprotein (VLDL) , 

35 low-density lipoprotein (LDL) , and high-density lipoprotein 
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{HDL) - are usually classified in terms of physicochemical 
properties (e.g., density after centrifugation) . 
Chylomicrons, the largest lipoproteins, carry exogenous 
triglyceride from the intestine via the thoracic duct to 
5 the venous system and into peripheral sites. VLDL carries 
endogenous triglyceride primarily from the liver to the 
same peripheral sites for storage or use. Lipases quickly 
degrade the triglyceride in VLDL to produce intermediate 
density lipoproteins (IDL) and within 2 to 6 h, IDL is 

10 degraded further to generate LDL, which has a plasma half- 
life of 2 to 3 days. While the overall fate of LDL is 
unclear, the liver is responsible for removing 
approximately 70% and active receptor sites have been found 
on the surfaces of hepatocytes. 

15 Several monogenic conditions that lead to 

elevated levels of one or more serum lipoproteins have been 
defined and the responsible gene identified, including, 
e.g., hyperlipoproteinemia type I {OMIM 238600), familial 
hypercholesterolemia (OMIM 143890) , and familial defective 

20 apolipoprotein B (OMIM 107730) . However, in many cases the 
etiology is unknown and there is strong evidence for 
additional uncharacterized loci. 

For example, Zuliani et al . , Arterioscler . 
Thromb. Vase. Biol. 19:802-809 (1999) identified a 

25 Sardinian family with a recessive form of 

hypercholesterolemia with the clinical features of familial 
hypercholesterolemia (OMIM 603813), and found that 
previously identified genes were not responsible for this 
disorder. They proposed that in this new lipid disorder, a 

30 recessive defect causes a selective impairment of the LDL 
receptor function in the liver. Ciccarese et al . , Am. J. 
Hum. Genet. 66:453-460 (2000) recently mapped this novel 
disease locus. 

Another example is designated familial combined 

35 hyperlipidemia (OMIM 144250) which affects approximately 1- 
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2% Of the population in the Western world. This disorder 
can have its basis in mutation in several novel genes, two 
of which have been mapped to chromosome 1 (Pajukanta et 
al., Nature Genet. 18:369-373 (1998)) and chromosome 11 

5 (Aouizerat et al . . Am. J. Hum. Genet. 65, 397-412 (1999)). 
The high frequency of this disorder suggests that most, if 
not all, hyperlipidemias are of multifactorial genetic 
etiology. 

As yet a further example, primary schlerosing 
10 cholangitis (PSC) is a disorder characterized by a patchy 
obliterative inflammatory fibrosis of the large bile ducts. 
Chronic inflammation leads to extensive bile duct 
strictures, cholestasis, and gradual progression to biliary 
cirrhosis. PSC occurs most often in young men and is 
15 commonly associated with inflammatory bowel disease, 
especially ulcerative colitis. The onset is usually 
insidious, with gradual, progressive fatigue, pruritus, and 
jaundice. There is no specific therapy for sclerosing 
cholangitis, and liver transplantation is the only apparent 
20 cure . 

The etiology of PSC is not known, but both 
genetic and immunologic abnormalities have been implicated. 
However, the frequency of HLA-B8 and HLA-DT2, which are 
associated with a number of autoimmune diseases, is higher 
25 in PSC than normal individuals. Prochazka et al . , New Eng. 
J. Med. 322:1842-1844 (1990) found that 100% of 29 patients 
with primary sclerosing cholangitis carried the HLA-DRw52a 
antigen, which is normally present in 35% of the 
population. 

30 As a still further example, sarcoidosis is a ■ 

disease of unknown cause characterized .by non-caseat ing 
granulomas in one or more organ systems. These granul^^mas 
may resolve completely or proceed to fibrosis. The disorder 
is systemic, but the liver is affected in approximately 75% 

35 of cases. Sarcoidosis occurs mainly in persons aged 20 to 
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4 0 yr and is most common in Northern Europeans and American 
blacks. The lifetime risk of developing sarcoidosis is 
particularly high among Swedish men (1.15%), Swedish women 
(1.6%), and African Americans (2,4%). 

5 The much greater frequency in African Americans 

relative to the United States population overall suggests a 
genetic contribution to etiology. Early research studying 
familial aggregation indicated that the disease may have a 
nongenetic basis because the family pattern did not conform 

10 to a simple Mendelian mode of inheritance {Allison, Sth. 
Med. J. 57: 27-32 (1964)). However, Headings et al . , Ann. 
N.Y. Acad. Sci. 278:377-385 (1976) favored multifactorial 
genetic inheritance of susceptibility. Nowack et al . , 
Arch. Intern. Med. 147:481-483 (1987), found an unusually 

15 high frequency of HLA-DR5 in a study of 440 patients with 
sarcoidosis in Marburg, Germany. They also concluded that 
the role of an environmental or infectious agent triggering 
sarcoidosis cannot be envisaged without considering 
genetically linked cof actors. 

20 Other significant diseases of liver are also 

believed to have a genetic, typically polygenic, etiologic 
component. These diseases include, e.g.,. primary biliary 
cirrhosis, Zellweger syndrome, cholestasis- lymphedema 
syndrome, Alstrom syndrome, primary pulmonary 

25 hypertension, Berardinelli -Seip congenital lipodystrophy, 
iron overload in Africa, neonatal cholestatic hepatitis, 
autosomal recessive KID syndrome, familial 
hypotransferrinemia, type I congenital dyserythropoietic 
anemia, porphyria variegata, Finnish lactic acidosis with 

30 hepatic hemosiderosis, Rotor syndrome, essential 
hypertension, ARC syndrome, type II conjugated 
hyperbilirubinemia, Lambert syndrome, ichthyosis congenita 
with biliary atresia, Kabuki make-up syndrome, Meckel 
syndrome, cerebral aneurysm- cirrhosis syndrome, glycogen 

35 storage diseases, polycystic kidney and hepatic disease. 
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isolated Caroli disease, trisomy 18-like syndrome, Osler- 
Rendu-Weber syndrome 3, fatdl intrahepatic cholestasis, 
Coach syndrome, type C Niemann-Pick disease, hereditary 
persistence of fetal hemoglobin and hepatocellular cancer. 
5 Altered responses to a variety of infectious 

agents that target the liver, especially acute viral 
hepatitis, have also been shown or are suspected to have 
genetic bases or contributions. In addition to 
differential susceptibility to primary infectious agents, 

10 these altered responses include predisposition to 

complicating conditions following contact with particular 
infectious agents. These include, e.g., development of 
.hepatocellular carcinoma 2 correlated with Hepatitis B 
infection, and severe hepatic fibrosis following 

15 Schistosoma mansoni infection. 

The central role of the liver in drug metabolism 
results in exposure of this organ to a large variety of 
potentially toxic chemical agents and metabolites. These 
include naturally occurring plant alkaloids and mycotoxins, 

20 industrial chemicals, and, additionally, pharmacologic 
agents used in treating disease. The range of 
manifestations of toxin- and drug- induced liver disease are 
virtually as broad as the range of acute and chronic 
disorders and have also been shown or suspected to have 

25 genetic bases or contributions. 

Such interactions between drugs and genotype have 
been shown in the response, e.g., to the anticonvulsant 
phenytoin, which can cause severe hepatitis-like disease in 
individuals who are impaired in the ability to detoxify a 

30 metabolite of phenytoin in the liver, and in the response 
to the drug sodium valproate, which can produce severe 
hepatotoxicity in certain individuals. The abnormal 
responses to both of these drugs are believed to be 
influenced by underlying genetic factors. 

35 The human genome -derived single exon nucleic acid 
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probes and microarrays of the present invention are useful 
for predicting, diagnosing, grading, staging, monitoring 
and prognosing diseases of human liver, particularly those 
diseases with polygenic etiology. With each of the single 
5 exon probes described herein shown to be expressed at 

detectable levels in human liver, and with about 2/3 of the 
probes identifying novel genes, the single exon microarrays 
of the present invention provide exceptionally high 
informational content for such studies. 

10 For example, diagnosis (including differential 

diagnosis among clinically indistinguishable disorders, 
such as cirrhosis) , staging, and/or grading of a disease 
can be based upon the quantitative relatedness of a patient 
gene expression profile to one or more reference expression 

15 profiles known to be characteristic of a given liver 
disease, or to specific grades or stages thereof. 

In one embodiment, the patient gene expression 
profile is generated by hybridizing nucleic acids obtained 
directly or indirectly- from transcripts expressed in the 

20 patient's liver to the genome-derived single exon 

microarray of the present invention. Reference profiles 
are obtained similarly, using nucleic acids obtained 
directly or indirectly from transcripts expressed by liver 
of individuals with known liver disease. Methods for 

25 quantitatively relating gene expression profiles, without 
regard to the function of the protein encoded by the gene, 
are disclosed in WO 99/58720, incorporated herein by 
reference in its entirety. 

In another approach, the genome-derived single 

30 exon probes and microarrays of the present invention can be 
used to interrogate genomic DNA, rather than pools of 
expressed message; this latter approach permits 
predisposition to and/or prognosis of liver disease to be 
assessed through the massively parallel determination of 

35 altered copy number, deletion, or mutation in the patient's 
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genome of axons known to be expressed in human liver. The 
algorithms set forth in WO 99/58720 can be applied to such 
genomic profiles without regard to the function of the 
protein encoded by the interrogated gene. 
5 The utility is specific to the probe; at 

sufficiently high hybridization stringency, which 
stringencies are well known in the art - see Ausubel et al. 
and Maniatis et al. - each probe reports the level of 
expression of message specifically containing that ORF. 
10 It should be appreciated, however, that the 

probes of the present invention, for which expression in 
the Fetal liver has been demonstrated are useful for both 
measurement in the Fetal liver and for survey of expression 
in other tissues. 
15 Significant among such advantages is the presence 

of probes for novel genes. 

As mentioned above and further detailed in 
Examples 1 and 2, the methods described enable ORFs which 
are not present in existing expression databases to be 
20 identified. And the fewer the number of tissues in which 
the ORF can be shown to be expressed, the more likely the 
ORF will prove to be part of a novel gene: as further 
discussed in Example 2, ORFs whose expression was 
measurable in only a single of the tested tissues were 
25 represented in existing expression databases at a rate of 
only 11%, whereas 36% of ORFs whose expression was 
measurable in 9 tissues were present in existing expression 
databases, and fully 45% of those ORFs expressed in all ten 
tested tissues were present in existing expressed sequence 
30 databases. 

Either as tools for measuring gene expression or 
tools for surveying gene expression, the genome -derived 
single exon probes of the present invention have 
significant advantages over the cDNA or EST-based probes 
35 that are currently available for achieving these utilities. 
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The genome -de rived single exon probes of the 

present invention are useful in constructing genome-derived 

single exon microarrays; the genome -derived single exon 

microarrays, in turn, are useful devices for measuring and 

5 for surveying gene expression in the human. 

Gene expression analysis using microarrays - 
conventionally using microarrays having probes derived from 
expressed message — is well-established as useful in the 
biological research arts (see Lockhart et al . Nature 405, 

10 827-836) . 

Microarrays have been used to determine gene 
expression profiles in cells in response to drug treatment 
(see, for example, Kaminski et al . , "Global Analysis of 
Gene Expression in Pulmonary Fibrosis Reveals Distinct 

15 Programs Regulating Lung Inflammation and Fibrosis," Proc, 
Wati- Acad. Sci. USA 97 (4 ): 1778 -83 (2000); Bartosiewicz et 
al . , "Development of a Toxicological Gene Array and 
Quantitative Assessment of This Technology, " Arch. Biochem. 
Biophys. 376(1) :66-73 (2000)), viral infection (see for 

20 example, Geiss et al . , "Large-scale Monitoring of Host Cell 
Gene Expression During HIV-1 Infection Using cDNA 
Microarrays," l^iroiogy 266 ( 1 ) : 8- 16 (2000)) and during cell 
processes such as differentiation, senescence and apoptosis 
(see, for example, Shelton et al . , "Microarray Analysis of 

25 Replicative Senescence," Curr. Biol. 9{17):939-45 (1999); 
Voehringer et al . , "Gene Microarray Identification of Redox 
and Mitochondrial Elements That Control Resistance or 
Sensitivity to Apoptosis," Proc . Natl, Acad. Sci. USA 
97 (6) :2680-5 (2000) ) . 

30 Microarrays have also been used to determine 

abnormal gene expression in diseased tissues (see, for 
example, Alon et al . , "Broad Patterns of Gene Expression 
Revealed by Clustering Analysis of Tumor and Normal Colon 
Tissues Probed by Oligonucleotide Arrays," Proc. Natl. 

35 Acad. Sci. USA 96 (12) -.6745-50 (1999); Perou et al . , 
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"Distinctive Gene Expression Patterns in Human Mammary 

Epithelial Cells and Breast Cancers, Proc. Natl, Acad. Scl . 
USA 96 (16) : 9212-7 {1999); Wang et al , , "Identification of 
Genes Differentially Over-expressed in Lung Squamous Cell 
5 Carcinoma Using Combination of cDNA Subtraction and 
Microarray Analysis," Oncogene 19 (12 ): 1519-28 (2000); 
Whitney et ai , , "Analysis of Gene Expression in Multiple 
Sclerosis Lesions Using cDNA Microarrays , " Ann. Neurol. 
46(3):425-8 (1999)), in drug discovery screens (see, for 

10 example, Scherf et al . , "A Gene Expression Database for the 
Molecular Pharmacology of Cancer," Nat. Genet. 24(3):236-44 
(2000) ) and in diagnosis to determine appropriate treatment 
strategies (see, for example, Sgroi et al . , "In vivo Gene 
Expression Profile Analysis of Human Breast Cancer 

15 Progression," Cancer Res . 59 (22 ): 5656 -61 (1999)). 

In microarray-based gene expression screens of 
pharmacological drug candidates upon cells, each probe 
provides specific useful data. In particular, it should be 
appreciated that even those probes that show no change in 

20 expression are as informative as those that do change, 
serving, in essence, as negative controls. 

For example, where gene expression analysis is 
used to assess toxicity of chemical agents on cells, the 
failure of the agent to change a gene's expression level is 

25 evidence that the drug likely does not affect the pathway 
of which the gene's expressed protein is a part. 
Analogously, where gene expression analysis is used to 
assess side effects of pharmacological agents - whether in 
lead compound discovery or in subsequent screening of lead 

30 compound derivatives - the inability of the agent to alter 
a gene's expression level is evidence that, the drug does 
not affect the pathway of which the gene's expressed \ 
protein is a part. 

WO 99/58720 provides methods for quantifying the 

35 relatedness of a first and second gene expression profile 
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and for ordering the relatedness of a plurality of gene 
expression profiles. The methods so described permit 
useful information to be extracted from a greater 
percentage of the individual gene expression measurements 
5 from a microarray than methods previously used in the art. 

Other uses of microarrays are described in 
Gerhoid et al., Trends Biochem. Sci. 24 { 5) : 168 - 173 (1999) 
and Zweiger, Trends Biotechnol . 17 { 11) : 429 -436 (199^); 
Schena et al . 

10 The invention particularly provides genome - 

derived single -exon probes known to be expressed in Fetal 
liver. The individual single exon probes can be 

provided in the form of substantially isolated and purified 
nucleic acid, typically, but not necessarily, in a quantity 

15 sufficient to perform a hybridization reaction. 

Such nucleic acid can be in any form directly 
hybridizable to the message that contains the probe's ORF, 
such as double stranded DNA, single- stranded DNA 
complementary to the message, single -stranded RNA 

20 complementary to the message, or chimeric DNA/RNA molecules 
so hybridizable. The nucleic acid can alternatively or 
additionally include either nonnative nucleotides, 
alternative internucleotide linkages, or both, so long as 
complementary binding can be obtained. For example, probes- 

25 can include phosphorothioates , methylphosphonates , 

morpholino analogs, and peptide nucleic acids (PNA) , as are 
described, for example, in U.S. Patent Nos. 5,14 2,047; 
5,235,033; 5,166,315; 5,217,866; 5,184,444; 5,861,250. 

Usefully, however, such probes are provided in a 

30 form and quantity suitable for amplification, where the 
amplified product is thereafter to be used in the 
hybridization reactions that probe gene expression. 
Typically, such probes are provided in a form and quantity 
suitable for amplification by PGR or by other well known 

35 amplification technique. One such technique additional to 
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PGR is rolling circle amplification, as is described, inter 
alia, in U.S. Patent Nos . 5,854,033 and 5,714,320 and 
international patent publicatiems WO 97/19193 and 
WO 00/15779. As is well understood, where the probes are 
5 to be provided in a form suitable for amplification, the 
range of nucleic acid analogues and/or internucleot ide 
linkages will be constrained by the requirements and nature 
of the amplification enzyme. 

Where the probe is to be provided in form 

10 suitable for amplification, the quantity need not be 
sufficient for direct hybridization for gene expression 
analysis, and need be sufficient only to function as an 
amplification template, typically at least about 1, 10 or 
100 pg or more. 

15 Each discrete amplifiable probe can also be 

packaged with amplification primers, either in a single 
composition that comprises probe template and. primers, or 
in a kit that comprises such primers ■ separately packaged 
therefrom. As earlier mentioned, the ORF-specific 

20 5' primers used for genomic amplification can have a first 
common sequence added thereto, and the ORF-specific 3' 
primers used for genomic amplification can have a second, 
different, common sequence added thereto, thus permitting, 
in this embodiment, the use of a single set of 5* and 3' 

25 primers to amplify any one of the probes. The probe 

composition and/or kit can also include buffers, enzyme, 
etc., required to effect amplification. 

As mentioned earlier, when intended for use on a 
genome -derived single exon microarray of the present 

30 invention, the genome -derived single exon probes of the 
present invention will typically average at least about 
100, 200, 300, 400 or 500 bp in length, including (and 
typically, but not necessarily centered about) the ORE. 
Furthermore, when intended for use on a genome -de rived 

35 single exon microarray of the present invention, the 
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genome -derived single exon probes of the present invention 

will typically not contain a detectable label. 

When intended for use in solution phase 

hybridization, however — that is, for use in a 

5 hybridization reaction in which the probe is not first 

bound to a support substrate (although the target may 

indeed be so bound) - length constraints that are imposed 

in microarray-based hybridization approaches will be 

relaxed, and such probes will typically be labeled. 

10 In such case, the only functional constraint that 

dictates the minimum size of such probe is that each such 
probe must be capable of specifically identifying in a 
hybridization reaction the exon from which it is drawn.- In 
theory, a probe of as little as 17 nucleotides is capable 

15 of uniquely identifying its cognate sequence in the human 
genome. For hybridization to expressed message — a subset 
of target sequence that is much reduced in complexity as 
compared to genomic sequence - even fewer nucleotides are 
required for specificity, 

20 Therefore, the probes of the present invention 

can include as few as 20, 25 or 50 bp or ORF, or more. In 
particular embodiments, the ORF sequences are given in SEQ 
ID NOS. 12,674 - 25,129, respectively, for probe SEQ ID 
NOS. 1 - 12,673. The minimum amount of ORF required to be 

25 included in the probe of the present invention in order to 
provide specific signal in either solution phase or 
microarray-based hybridizations can readily be determined 
for each of ORF SEQ ID NOS. 12,674 - 25,129 individually by 
routine experimentation using standard high stringency 

30 conditions. 

Such high stringency conditions are described, 
inter alia, in Ausubel et al . and Maniatis et al . For 
microarray-based hybridization, standard high stringency 
conditions can usefully be 50% formamide, 5X SSC, 0.2 \iq/\il 

35 poly(dA) , 0,2 pg/pl human Cotl DNA, and 0.5 % SDS, in a 
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humid oven at 42°C overnight, followed by successive washes 
of the microarray in IX SSC, 0.2% SDS at 55°C for 5 
minutes, and then 0 , IX SSC, 0.2% SDS, at SS'^C for 20 
minutes. For solution phase hybridization, standard high 
stringency conditions can usefully be aqueous hybridization 
at 65°C in 6X SSC. Lower stringency conditions, suitable 
for cross-hybridization to mRNA encoding structurally- and 
functionally-related proteins, can usefully be the same as 
the high stringency conditions but with reduction in 
temperature for hybridization and washing to room 
temperature (approximately 25°C) . 

When intended for use in solution phase 
hybridization, the maximum size of the single exon probes 
of the present invention is dictated by the proximity of 
other expressed exons in genomic DNA: although each single 
exon probe can include intergenic and/or intronic material 
contiguous to the ORF in the human genome, each probe of 
the present invention will include portions of only one 
expressed exon. 

Thus, each single exon probe will include no more 
than about 25 kb of contiguous genomic sequence, more 
typically no more than about 20 kb of contiguous genomic 
sequence, more usually no more than about 15 kb, even more 
usually no more than about 10 kb. Usually, probes that are 
maximally about 5 kb will' be used, more typically no more 
than about 3 kb. 

It will be appreciated that the Sequence Listing 
appended hereto presents, by convention, only that strand 
of the probe and ORF sequence that can be directly 
translated reading from 5' to 3 ' end. As would be well 
understood by one of skill in the art, single stranded 
probes must be complementary in sequence to the ORF as 
present in an mRNA; it is well within the skill in the art 
to determine such complementary sequence. It will further 
be understood that double stranded probes can be used in 
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both solution-phase hybridization and microarray-based 
hybridization if suitably denatured. 

Thus, it is an aspect of the present invention to 
provide single- stranded nucleic acid probes that have 
sequence complementary to those described herein above and 
below, and double -stranded probes one strand of which has 
sequence complementary to the probes described herein. 

The probes can, but need not, contain intergenic 
and/or intronic material that flanks the ORF, on one or 
both sides, in the same linear relationship to the ORF that 
the intergenic and/or intronic material bears to the ORF in 
genomic DNA. The probes do not, however, contain nucleic 
acid derived from more than one expressed ORF. 

And when intended for use in solution 
hybridization, the probes of the present invention can 
usefully have detectable labels. Nucleic acid labels are 
well known in the art, and include, inter alia, radioactive 
labels, such as 'h, '^P, ^'p, ''s, '"i, '^^I; fluorescent 
labels, such as Cy3 , Cy5 , CyS.S, Cy7, SYBR® 

Green and other labels described in Haugland, 
Handbook of Fluorescent Probes and Research Chemicals, 7th 
ed.. Molecular Probes Inc., Eugene, OR (2000), or 
fluorescence resonance energy transfer tandem conjugates 
thereof; labels suitable for chemiluminescent and/or 
enhanced chemiluminescent detection; labels suitable for 
ESR and NMR detection; and labels that include one member 
of a specific binding pair, such as biotin, digoxigenin, or 
the like. 

The probes, either in quantity sufficient for 
hybridization or sufficient for amplification, can be 
provided in individual vials or containers. 

Alternatively, such probes can usefully be 
packaged as a plurality of such individual genome-derived 
single exon probes. 

When provided as a collection of plural 
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individual probes, the probes are typically made available 
in amplifiable form in a spatially-addressable ordered set, 
typically one per well of a microtiter dish. Although a 96 
well microtiter plate can be used, greater, efficiency is 
5 obtained using higher density arrays. 

If, as earlier mentioned, the ORF-specific 
5' primers used for genomic amplification had a first 
common sequence added thereto, and the ORF-specific 3' 
primers used for genomic amplification had a second, 
10 different, common sequence added thereto, a single set of 
5' and 3' primers can be used to amplify all of the probes 
from the amplifiable ordered set. 

Such collections of genome-derived single exon 
probes can usefully include a plurality of probes chosen 
15 for the common attribute of expression in the human Fetal 
liver . 

In such defined subsets, typically at least 50, 
60, 75, 80, 85, 90 or 95% or more of the probes will be 
chosen by their expression in the defined tissue or cell 
20 type . 

The single exon probes of the present invention, 
as well as fragments of the single exon probes comprising 
selectively hybridizable portions of the probe ORF, can be 
used to obtain the full length cDNA that includes the ORF 

25 by (i) screening of cDNA libraries; (ii) rapid 

amplification of cDNA ends ("RACE"); or (iii) other 
conventional means, as are described, inter alia, in 
Ausubel et al. and Maniatis et ai . 

It is another aspect of the present invention to 

30 provide genome -derived single exon nucleic acid microarrays 
useful for gene expression analysis, where the term 
"microarray" has the meaning given in the definitional \^ 
section of this description, supra. 

The invention particularly provides genome - 

35 derived single-exon nucleic acid microarrays comprising a 
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plurality of probes known to be expressed in human Fetal 
liver. In preferred embodiments, the present invention 
provides human genome -derived single exon microarrays 
comprising a plurality of probes drawn from the group 
5 consisting of SEQ ID NOS.: 1 - 12,673. 

When used for gene expression analysis, the 
genome-derived single exon microarrays provide greater 
physical informational density than do the genome -derived 
single exon microarrays that have lower percentages of 

10 probes known to be expressed commonly in the tested tissue. 
At a fixed probe density, for example, a given microarray 
surface area of the defined subset genome -derived single 
exon microarray can yield a greater number of expression 
measurements. Alternatively, at a given probe density, the 

15 same number of expression measurements can be obtained from 
a smaller substrate surface area. Alternatively, at a 
fixed probe density and fixed surface area, probes can be 
provided redundantly, providing greater reliability in 
signal measurement for any given probe. Furthermore, with 

20 a higher percentage of probes known to be expressed in the 
assayed tissue, the dynamic range of the detection means 
can be adjusted to reveal finer levels discrimination among 
the levels of expression. 

Although particularly described with respect to 

25 their utility as probes of gene expression, particularly as 
probes to be included on a genome -derived single exon 
microarray, each of the nucleic acids having SEQ ID NOS. : 1 
- 12,673 contains an open-reading frame, set forth 
respectively in SEQ ID NOS.: 12,674 - 25,129, that encodes 

30 a protein domain. Thus, each of SEQ ID NOS. 1 - 12,673 can 
be used, or that portion thereof in SEQ ID NOS. 12,674 - 
25,129 used, to express a protein domain by standard in 
vitro recombinant techniques. See Ausubel et al . and 
Maniatis et al . 

35 Additionally, kits are available commercially 
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that readily permit such nucleic acids to be expressed as 
protein in bacterial cells/, insect cells, or mammalian 
cells, as desired (e.g., HAT" Protein Expression & 
Purification System, ClonT^h Laboratories, Palo Alto, CA; 
5 Adeno-X™ Expression System/' ClonTech Laboratories, Palo 

Alto, CA; Protein Fusion & 'Purification (pMAL-) System, New 
England Biolabs, Beverley, 

Furthermore, shorter peptides can be chemically 
synthesized using commercial peptide synthesizing equipment 
10 and well known techniques.: Procedures are described, inter 
aiia, in Chan et ai . (eds.K, Fmoc Solid Phase Peptide 
Synthesis: A Practical Approach (Practical Approach Series, 
(Paper) )., Oxford Univ. Press (March 2000) (ISBNj - 
0199637245); Jones, Amino Acid and Peptide Sy nthesis 
15 (Oxford Chemistry Primers, No 7) , Oxford Univ. Press 

(August 1992) (ISBN: 0198556683); and Bodanszky, Principles 
of Peptide Synthesis (Springer Laboratory) , Springer Verlag 
(December 1993) (ISBN: 0387564314). 

It is, therefore, another aspect of the invention 
20 to provide peptides comprising an amino acid sequence 

translated from SEQ ID NOS . : 12,674 - 25,129. Such amino 
acid sequences are set out in SEQ ID NOS: 25,130 - 37,156. 
Any such recombinantly-expressed or synthesized peptide of 
at least 8, and preferably at least about 15, amino acids, 
25 can be conjugated to a carrier protein and used to generate 
antibody that recognizes the peptide. Thus, it is a 
further aspect of the invention to provide peptides that 
have at least 8, preferably at least 15, consecutive amino 
acids . 



30 

The following examples are offered by way of 
illustration and not by way of limitation. 

EXAMPLE 1 

35 Preparation of Single Exon Microarrays from ORFs Predicted 
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Bioinf ormatics Results 

All human BAG sequences in fewer than 10 pieces 
that had been accessioned in a five month period 
immediately preceding this study were downloaded from 
GenBank. This corresponds to -2200 clones, totaling -350 
MB of sequence, or approximately 10% of the human genome. 

After masking repetitive elements using the 
program CROSS_MATCH, the sequence was analyzed for open 
reading frames using three separate gene finding programs. 
The three programs predict genes using independent 
algorithmic methods developed on independent training sets: 
GRAIL uses a neural network, GENEFINDER uses a hidden 
Markoff model, and DICTION, a program proprietary to 
Genetics Institute, operates according to a different 
heuristic. The results of all three programs were used to 
create a prediction matrix across the segment of genomic 
DNA . 

The three gene finding programs yielded a range 
of results. GRAIL identified the greatest percentage of 
genomic sequence as putative coding region, 2% of the data 
analyzed. GENEFINDER was second, calling 1%, and DICTION 
yielded the least putative coding region, with 0.8% of 
genomic sequence called as coding region. 

The consensus data were as follows. GRAIL and 
GENEFINDER agreed on 0.7% of genomic sequence, GRAIL and 
DICTION agreed on 0.5% of genomic sequence, and the three 
programs together agreed on 0.2 5% of the data analyzed. 
That is, 0.25% of the genomic sequence was identified by 
all three of the programs as containing putative coding 
region . 

ORFs predicted by any two of the three programs 
{"consensus ORFs")' were assorted into "gene bins" using two 
criteria: (1) any 7 consecutive exons within a 25 kb window 
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were placed together in a bin as likely contributing to a 
single gene, and (2) all ORFs within a 25 kb window were 
placed together in a bin as likely contributing to a single 
gene if fewer than 7 exons were found within the 25 kb 
5 window. 

PGR 

The largest ORF from each gene bin that did not 
span repetitive sequence was then chosen for amplification, 

10 as were all consensus ORFs longer than 500 bp. This method 
approximated one exon per gene; however, a number of genes 
were found to be represented by multiple elements. 

Previously, we had determined that DNA fragments 
fewer than 250 bp in length do not bind well to the amino- 

15 modified glass surface of the slides used as support 
substrate for construction of microarrays; therefore, 
amplicons were designed in the present experiments to 
approximate 500 bp in length. 

Accordingly, after selecting the largest ORF per 

20 gene bin, a 500 bp fragment of sequence centered on the ORF 
was passed to the primer picking software, PRIMER3 
(available online for use at 

http://www-genome.wi.mit.edu/cgi-bin/primer/ ). A first 
additional sequence was commonly added to each ORF- unique 

25 5' primer, and a second, different, additional sequence was 
commonly added to each ORF-unique 3' primer, to permit 
subsequent reamplif ication of the amplicon using a single 
set of "universal" 5' and 3' primers, thus immortalizing 
the amplicon. The addition of universal priming sequences 

30 also facilitates sequence verification, and can be used to 
add a cloning site should some ORFs be found to warrant 
further study. 

The ORFs were then PGR amplified from genomic 
DNA, verified on agarose gels, and sequenced using the 

35 universal primers to validate the identity of the amplicon 
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to be spotted in the microarray. 

Primers were supplied by Operon Technologies 
(Alameda, CA) . PGR amplification was performed by standard 
techniques using human genomic DNA (Clontech, Palo Alto, 
■ 5 CA) as template. Each PGR product was verified by SYBR® 
green (Molecular Probes, Inc., Eugene, OR) staining of 
agarose gels, with subsequent imaging by Fluorimager 
(Molecular Dynamics, Inc., Sunnyvale, CA) . PGR 
amplification was classified as successful if a single band 
10 appeared. 

The success rate for amplifying ORFs of interest 
directly from genomic DNA using PGR was approximately 75%. 
FIG. 5 graphs the distribution of predicted ORF (exon) 
length and distribution of amplified PGR products, with ORF 

15 length shown in red and PGR product length shown in blue 

(which may appear black in the figure) . Although the range 
of ORF sizes is readily seen to extend to beyond 900 bp, 
the mean predicted exon size was only 229 bp, with a median 
size of 150 bp {n=9498) . With an average amplicon size of 

20 475 + 25 bp, approximately 50% of the average PGR 

amplification product contained predicted coding region, 
with the remaining 50% of the amplicon containing either 
intron, intergenic sequence, or both. 

Using a strategy predicated on amplifying about 

25 500 bp, it was found that long exons had a higher PGR 

failure rate. To address this, the bioinf ormatics process 
was adjusted to amplify 1000, 1500 or 2000 bp fragments 
from exons larger than 500 bp. This improved the rate of 
successful amplification of exons exceeding 500 bp, 

30 constituting about 9.2% of the exons predicted by the gene 
finding algorithms. 

Approximately 75% of the probes disposed on the 
array (90% of those that successfully PGR amplified) were 
sequence-verified by sequencing in both the forward and 

35 reverse direction using MegaBAGE sequencer (Molecular 
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Dynamics, Inc., Sunnyvale, CA) , universal primers, and 
standard protocols . 

Some genomic clones (BACs) yielded very poor PCR 
and sequencing results. The reasons for this are unclear, 

5 but may be related to the quality of early draft sequence 
or the inclusion of vector and host contamination in some 
submitted sequence data. 

Although the intronic and intergenic material 
flanking coding regions could theoretically interfere with 

10 hybridization during microarray experiments, subsequent 

empirical results demonstrated that differential expression 
ratios were not significantly affected by the presence of 
noncoding sequence. The variation in exon size was 
similarly found not to affect differential expression 

15 ratios significantly; however, variation in exon size was 
observed to affect the absolute signal intensity (data not 
shown) . 

The 350 MB of genomic DNA was, by the above - 
described process, reduced to 9750 discrete probes, which 

20 were spotted in duplicate onto glass slides using 

commercially available instrumentation {MicroArray Genii 
Spotter and/or MicroArray Genlll Spotter, Molecular 
Dynamics, Inc., Sunnyvale, CA) . Each slide additionally 
included either 16 or 32 E. coli genes, the average 

25 hybridization signal of which was used as a measure of 
background biological noise. 

Each of the probe sequences was BLASTed against 
the human EST data set, the NR data set, and SwissProt 
' GenBank (May 7, 1999 release 2.0.9). 

30 One third of the probe sequences (as amplified) 

produced an exact match (BLAST Expect ("E") values less 
than 1 e"^°.°) to either an EST (20% of sequences) or a l^nown 
mRNA (13% of sequences) . A further 22% of the probe 
sequences showed some homology to a known EST or mRNA 

35 (BLAST E values from 1 e'^ to 1 e"^^) . The remaining 45% of 
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the probe sequences showed no significant sequence homology 
to any expressed, or potentially expressed, sequences 
present in public databases. 

All of the probe sequences {as amplified) were 

5 then analyzed for protein similarities with the SwissProt 
database using. BLASTX, Gish et al . , Nature Genet. 3:266 
(1993) . The predicted functional breakdowns of the 2/3 of 
probes identical or homologous to known sequences are 
presented in Table 1. 

10 

Table 1 

Function of Predicted ORFs As Deduced From Comparative 
Sequence Analysis 



Total V6 chip V7 chip Function Predicted from 

Comparative Sequence 
Analysis 



211 


96 


115 


Receptor 


120 


43 


77 


Zinc Finger 


30 


11 


19 


Homeobox 


25 


9 


16 


Transcription Factor 


17 


11 


7 


Transcription 


118 


57 


61 


Structural 


95 


39 


56 


Kinase 


36 


18 


18 


Phosphatase 


83 


31 


52 


Ribosomal 


45 


19 


26 


Transport 


21 


17 


14 


Growth Factor 


17 


12 


5 


Cytochrome 


50 


33 


17 


Channel 



As can be seen, the two most common types of 
genes were transcription factors and receptors, making up 
15 2.2% and 1.8% of the arrayed elements, respectively. 
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EXAMPLE 2 

Gene Expression Measurements From Genome -Derived Single 
5 Exon Microarrays 

The two genome -derived single exon microarrays 
prepared according to Example 1 were hybridized in a series 

10 of simultaneous two-color fluorescence experiments to (1) 
Cy3- labeled cDNA synthesized from message drawn 
individually from each of brain, heart, liver, fetal liver, 
placenta, lung, bone marrow, HeLa, BT 474, or HBL 100 
cells, and (2) Cy5- labeled cDNA prepared from message 

15 pooled from all ten tissues and cell types, as a control in 
each of the measurements. Hybridization and scanning were 
carried out using standard protocols and Molecular Dynamics 
equipment . 

Briefly, mRNA samples were bought from commercial 

20 sources {Clontech, Palo Alto, CA and Amersham Pharmacia 

Biotech (APE)) . Cy3-dCTP and Cy5-dCTP (both from APB) were 
incorporated during separate reverse transcriptions of 1 tag 
of polyA* mRNA performed using 1 pg oligo (dT) 12-18 primer 
and 2 pg random 9mer primers as follows. After heating to 

25 VO^'C, the RNA: primer mixture was snap cooled on ice. After 
snap cooling on ice, added to the RNA to the stated final 
concentration was: IX Superscript II buffer, 0.01 M DTT, 
lOOpM dATP, 100 \iM dOTP, 100 pM dTTP, 50 pM dCTP, 50 pM 
Cy3-dCTP or Cy5-dCTP 50 pM, and 200 U Superscript II 

30 enzyme. The reaction was incubated for 2 hours at 42°C. 

After 2 hours, the first strand cDNA was isolated by adding 
1 U Ribonuclease H, and incubating for 30 minutes at 37°C. 
The reaction was then purified using a Qiagen PGR cleanup 
column, increasing the number of ethanol washes to 5. 

35 Probe was eluted using 10 mM Tris pH 8.5, 
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Using a spectrophotometer, probes were measured 
for dye incorporation. Volumes of both Cy3 and Cy5 cDNA 
corresponding to 50 pmoles of each dye were then dried in a 
Speedvac, resuspended in 30* "pi hybridization solution 
5 containing 50% formamide, 5X SSC, 0.2 pg/pl poly{dA), 0.2 
yg/pl human Cotl DNA, and 0.5' % SDS . 

Hybridizations were carried out under a 
coverslip, with the array placed in a humid oven at 42°C 
overnight. Before scanning, slides were washed in IX SSC, 
10 0.2% SDS at 55°C for 5 minutes, followed by 0 . IX SSC, 0.2% 
SDS, at 55°C for 20 minutes. Slides were briefly dipped in 
water and dried thoroughly under a gentle stream of 
nitrogen . 

Slides were scanned using a Molecular Dynamics 
15 Gen3 scanner, as described. Schena (ed. ) , Microarray 
Biochip: Tools and Technology , Eaton Publishing 
Company/BioTechniques Books Division (2 000) (ISBN: 
1881299376) . 

Although the use of pooled cDNA as a reference 

20 permitted the survey of a large number of tissues, it 
attenuates the measurement of relative gene expression, 
since every highly expressed gene in the tissue/cell type- 
specific fluorescence channel will be present to a level of 
at least 10% in the control channel. Because of this fact, 

25 both signal and expression ratios {the latter hereinafter, 
"expression" or "relative expression") for each probe were 
normalized using the average ratio or average signal, 
respectively, as measured across the whole slide. 

Data were accepted for further analysis only when 

30 signal was at least three times greater than biological 

noise, the latter defined by the average signal produced by 
the E. coli control genes. 

The relative expression signal for these probes 
was then plotted as function of tissue or cell type, and is 

35 presented in FIG. 6. 
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FIG. 6 shows the distribution of expression 
across a panel of ten tissues. The graph shows the number 
of sequence- verified products that were either not 
expressed ("0"), expressed in one or more but not all 
5 tested tissues ("1" - "9"), and expressed in all tissues 
tested ("10") . 

Of 9999 arrayed elements on the two microarrays 
{including positive and negative controls and "failed" 
products) , 2353 (51%) were expressed in at least one tissue 
10 or cell type. Of the gene elements showing significant 
signal — where expression was scored as "significant" if 
the normalized Cy3 signal was greater than 1, representing 
signal 5-fold over biological noise (0.2) - 39% (991) were 
expressed in all 10 tissues. The next most common class 
15 (15%) consisted of gene elements expressed in only a single 
tissue . 

The genes expressed in a single tissue were 
further analyzed, and the results of the analyses are 
compiled in FIG. 7. 

20 FIG. 7A is a matrix presenting the expression of 

all verified sequences that showed expression greater than 
3 in at least one tissue. Each clone is represented by a 
column in the matrix. Each of the 10 tissues assayed is 
represented by a separate row in the matrix, and relative 

25 expression of a clone in that tissue is indicated at the 
respective node by intensity of green shading, with the 
intensity legend shown in panel B. The top row of the 
matrix ("EST Hit") contains "bioinf ormatic" rather than 
"physical" expression data — that is, presents the results 

30 returned by query of EST, NR and SwissProt databases using 
the probe sequence. The legend for "bioinf ormatic 
expression" (i.e., degree of homology returned) is 
presented in panel C. Briefly, white is known, black is 
novel, with gray depicting nonidentical with significant 

35 homology (white: E values < le-100; gray: E values from le- 
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05 to le-99; black: E values > le-05) . 

As FIG. 7 readily shows, heart and brain were 
demonstrated to have the greatest numbers of genes that 
were shown to be uniquely expressed in the respective 

5 tissue. In brain, 200 uniquely expressed genes were 

identified; in heart, 150. The remaining tissues gave the 
following figures for uniquely expressed genes: liver, 100; 
lung, 70; fetal liver, 150; bone marrow, 75; placenta, 100; 
HeLa, 50; HBL, 100; and BT474, 50. 

10 It was further observed that there were many more 

"novel" genes among those that were up-regulated in only 
one tissue, as compared with those that were down- regulated 
in only one tissue. In fact, it was found that ORFs whose 
expression was measurable in only a single of the tested 

15 tissues were represented in sequencing databases at a rate 
of only 11%, whereas 36% of the ORFs whose expression was 
measurable in 9 of the tissues were present in public 
databases. As for those ORFs expressed in all ten tissues, 
fully 45% were present in existing expressed sequence 

20 databases. These results are not unexpected, since genes 
expressed in a greater number of tissues have a higher 
likelihood of being, and thus of having been, discovered by 
EST approaches . 

25 Comparison of Signal from Known and Unknown Genes 

The normalized signal of the genes found to have 
high homology to genes present in the GenBank human EST 
database were compared to the normalized signal of those 
genes not found in the GenBank human EST database. The 

30 data are shown in FIG. 8. 

FIG. 8 shows the normalized Cy3 signal intensity 
for all sequence-verified products with a BLAST Expect 
("E") value of greater than le-30 (designated "unknown") 
upon query of existing EST, NR and SwissProt databases, and 

35 shows in blue the normalized Cy3 signal intensity for all 
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sequence-verified products with a BLAST Expect value of 
less than le-30 ("known"). Note that biological background 
noise has an averaged normalized Cy3 signal intensity of 
0.2. 

5 As expected, the most highly expressed of the 

ORFs were "known" genes. This is not surprising, since 
very high signal intensity correlates with very commonly- 
expressed genes, which have a higher likelihood of being 
found by EST sequence. 

10 However, a significant point is that a large 

number of even the high expressers were "unknown" , Since 
the genomic approach used to identify genes and. to confirm 
their expression does not bias exons toward either the 3 ' 
or 5' end of a gene, many of these high expression genes 

15 will not have been detected in an end-sequenced cDNA 
library . 

The significant point is that presence of the 
gene in an EST database is not a prerequisite for 
incorporation into a genome -derived microarray, and 
20 further, that arraying such "unknown" exons can help to 
assign function to as-yet undiscovered genes. 

Verification of Gene Expression 

To ascertain the validity of the approach 

25 described above to identify genes from raw genomic 

sequence, expression of two of the probes was assayed using 
reverse transcriptase polymerase chain reaction {RT PGR) 
and northern blot analysis. 

Two microarray probes were selected on the basis 

30 of exon size, prior sequencing success, and tissue-specific 
gene expression patterns as measured by the microarray 
experiments. The primers originally used to amplify the 
two respective ORFs from genomic DNA were used in RT PGR 
against a panel of tissue-specific cDNAs (Rapid-Scan gene 

35 expression panel 24 human cDNAs) (OriGene Technologies, 
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Inc. , Rockville, MD) . 

Sequence AL079300_1 was shown by microarray 
hybridization, to be present in cardiac tissue, and sequence 
AL031734_1 was shown by microarray experiment to be present 
5 in placental tissue {data not shown) . RT-PCR on these two 
sequences confirmed the tissue-specific gene expression as 
measured by microarrays, as ascertained by the presence of 
a correctly sized PGR product from the respective tissue 
type cDNAs . 

10 Clearly, all microarray results cannot, and 

indeed should not, be confirmed by independent assay 
methods, or the high throughput, highly parallel advantages 
of microarray hybridization assays will be lost. However, 
in addition to the two RT-PCR results presented above, the 

15 observation that 1/3 of the arrayed genes exist in 

expression databases provides powerful confirmation of the 
power of our methodology — which combines bioinf ormatic 
prediction with expression confirmation using genome- 
derived single exon microarrays - to identify novel genes 

20 from raw genomic data. 

To verify that the approach further provides 
correct characterization of the expression patterns of the 
identified genes, a detailed analysis was performed of the 
microarrayed sequences that showed high signal in brain. 

25 For this latter analysis, sequences that showed 

high (normalized) signal in brain, but which showed very 
low (normalized) signal (less than 0.5, determined to be 
biological noise) in all other tissues, were further 
studied. There were 82 sequences that fit these criteria, 

30 approximately 2% of the arrayed elements. The 10 sequences 
showing the highest signal in brain in microarray 
hybridizations are detailed in Table 2, along with assigned 
function, if known or reasonably predicted. 

35 Table 2 
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Function of the Most Highly 
Expressed Genes Expressed Only in Brain 

Microarray Normal Expressi Homology Gene Function 

Name Signal present GenBank 

in 

GenBank 


AP000217-1 


5,2 


+ 7.7 

- 


High 


S-100 protein, 
b-chain, Ca^* 
binding protein 
expressed in 
central nervous 
system 


AP000047-1 


2 . 3 




High 


Unknown 
Function 


AC006548-9 


1 . 7 




High 


Similar to 
mouse membrane 
glyco -protein 
M6, expressed 
in central 

lliz I. V KJ\Jia oysUcKi 


AC007245-5 


1.5 




High 


Similar to 
amphiphysin, a 
synaptic 
vesicle- 
associated 
protein. Ref 21 


L44140-4 


1.2 


+2.0 


High 


Endothelial 
actin-binding 
protein found 
in nonmuscle 
f ilamin 
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AC004689-9 


1.2 


+ 3.5 


High 


Protein 
Phosphatase 
PP2A, neuronal/ 
downregulates 
activated 
protein kinases 


AL031657-1 


1.2 


+ 3.0 


High 


Unknown 
function/ 
Contains the 
anhyrin motif, 
a common 
protein 
sequence motif 


AC009266-2 


1.1 


+ 3 . 7 


Low 


Low homology to 
the 

Synaptotagmin I 
protein in 
rat/present at 
low levels 
throughout rat 
brain 


AP000086-1 


1.0 


+ 2.7 


Low 


Unknown, very 
poor homology 
to collagen 


AC004689-3 


1.0 




High 


Protein 
Phosphatase 
PP2A, neuronal/ 
downregulates 
activated 
protein kinases 



Of the ten sequences studied by these latter 
confirmatory approaches, eight were previously known. Of 
these eight, six had previously been reported to be 
5 important in the central nervous system or brain. The exon 
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giving the highest signal (AP00217-1) was found to be the 
gene encoding an SIOOB Ca^* binding protein, reported in 
the literature to be highly and uniquely expressed in the 
central nervous system. Heizmann, Neurochem. Res, 9:1097 

5 (1997) . 

A number of the brain-specific probe sequences 
{including AC006548-9, AC009266-2) did not have homology to 
any known human cDNAs in GenBank but did show homology to 
rat and mouse cDNAs . Sequences AC004689-9 and AC004689-3 

iO were both found to be phosphatases present in neurons 
(Millward et al , , Trends Biocbem. Sci . 24 (5) : 186-191 
(1999)). Two microarray sequences, AP000047-1 and 
AP000086-1 have unknown function, with AP000086-1 being 
absent from GenBank. Functionality can now be narrowed 

15 down to a role in the central nervous system for both of 
these genes, showing the power of designing microarrays in 
this fashion. 

Next, the function of the chip sequences with the 
highest (normalized) signal intensity in brain, regardless 

20 of expression in other tissues, was assessed. In this 
latter analysis, we found expression of many more common 
genes, since the sequences were not limited to those 
expressed only in brain. For example, looking at the 2 0 
highest signal intensity spots in brain, 4 were similar to 

25 tubulin (AC00807905 ; AF146191-2; AC007664-4; AF14191-2), 2 
were similar to actin (AL035701-2 ; AL034402-1), and 6 were 
found to be homologous to glyceraldehyde- 3 -phosphate 
dehydrogenase (GAPDH) (AL035604-1; Z86090-1; AC006064-L, 
AC006064-K; AC035604-3; AC006064-L). These genes are often 

30 used as controls or housekeeping genes in microarray 
experiments of all types. 

Other interesting genes highly expressed in brain 
were a ferritin heavy chain protein, which is reported in 
the literature to be found in brain and liver (Joshi et 

35 al., J. Neurol, Sci, 134 (Suppl) : 52 - 56 (1995)), a result 
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duplicated with the array. Other highly expressed chip 
sequences included a translation elongation factor ID 
(AC007564-4) , a DEAD-box homolog {AL023804 -4) , and a Y- 
chromosome RNA-binding motif (Chai et al . , Genomics 
5 49(2) :283-89 (1998) ) {AC007320-3) . A low homology analog 
(AP00123-1/2) to a gene, DSCRl , thought to be involved in 
trisomy 21 (Down's syndrome), showed high expression in 
both brain and heart, in agreement with the literature 
(Fuentes et al . , Mol . Genet. 4 (10) : 1935-44 (1995)). 

10 As a further validation of the approach, we 

selected the BAG AC006064 to be included on the array. 
This BAG was known to contain the GAPDH gene, and thus 
could be used as a control for the ORF selection process. 
The gene finding and exon selection algorithms resulted in 

15 choosing 25 exons from BAG AG006064 for spotting onto the 
array, of which four were drawn from the GAPDH gene. Table 
3 shows the comparison of the average expression ratio for 
the 4 exons from BAC006064 compared with the average 
expression ratio for 5 different dilutions of a 

20 commercially available GAPDH cDNA (Glontech) . 



Table 3 



Comparison of Expression Ratio, for each 
tissue, of GAPDH 




AC006064 (n = 4) 


Control ( n = 5) 


Bone Marrow 


-1.81 ± 0.11 


-1.85 ± 0.08 


Brain 


-1.41 ± 0.11 


-1.17 ± 0.05 


BT474 


1.85 ± 0.09 


1.66 ± 0.12 


Fetal Liver 


-1.62 ± 0.07 


-1.41 ± 0.05 


HBLIOO 


1.32 ± 0.05 


2.64 ± 0.12 


Heart 


1.16 + 0.09 


1.56 + 0.10 


HeLa 


1.11 ±0.06 


1.30 ± 0.15 


Liver 


-1.62 ± 0.22 


-2.07 ± 
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Lung 


-4.95 ± 0.93 


-3.75 ± 0.21 


Placenta 


-3.56 ± 0.25 


-3.52 ± 0.43 



Each tissue shows excellent agreement between the 
experimentally chosen exons and the control, again 
5 demonstrating the validity of the present exon mining 

approach. In addition, the data also show the variability 
of expression of GAPDH within tissues, calling into 
question its classification as a housekeeping gene and 
utility as a housekeeping control in microarray 
10 experiments. 

EXAMPLE 3 

Representation of Sequence and Expression Data as a 
"Mondrian" 

15 

For each genomic clone processed for microarray 
as above -described, a plethora of information was 
accumulated, including full clone sequence, probe sequence 
within the clone, results of each of the three gene finding 

20 programs, EST information associated with- the probe 
sequences, and microarray signal and expression for 
multiple tissues, challenging our ability to display the 
information . 

Accordingly, we devised a new tool for visual 

25 display of the sequence with its attendant annotation 
which, in deference to its visual similarity to the 
paintings of Piet Mondrian, is hereinafter termed a 
"Mondrian". FIGS. 3 and 4 present the key to the 
information presented on a Mondrian. 

30 FIG. 9 presents a Mondrian of BAG AC008172 {bases 

25,000 to 130,000 shown), containing the carbamyl phosphate 
synthetase gene {AF154830 . 1 ) . Purple background within the 
region shown as field 81 in FIG. 3 indicates all 37 known 
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exons for this gene. 

As can be seen, GRAIL II successfully identified 
27 of the known exons (73%) , GENEFINDER successfully 
identified 37 of the known exons (100%), while DICTION 
5 identified 7 of the known exons (19%) . 

Seven of the predicted exons were selected for 
physical assay, of which 5 successfully amplified by PGR 
and were sequenced. These five exons were all found to be 
from the same gene, the carbamyl phosphate synthetase gene 

10 {AF154830. 1) . 

The five exons were arrayed, and gene expression 
measured across 10 tissues. As is readily seen in the 
Mondrian, the five chip sequences on the array show 
identical expression patterns, elegantly demonstrating the 

15 reproducibility of the system. 

FIG. 10 is a Mondrian of BAG AL049839. We 
selected 12 exons from this BAG, of which 10 successfully 
sequenced, which were found to form between 5 and 6 genes. 
Interestingly, 4 of the genes on this BAG are protease 

20 inhibitors. Again, these data elegantly show that exons 
selected from the same gene show the same expression 
patterns, depicted below the red line. From this figure, 
it is clear that our ability to find known genes is very 
good. A novel gene is also found from 86.6 kb to 88.6 kb, 

25 upon which all the exon finding programs agree. We are 
confident we have two exons from a single gene since they 
show the same expression patterns and the exons are 
proximal to each other. Backgrounds in the following 
colors indicate a known gene {top to bottom) : 

30 red = kallistatin protease inhibitor (P29622) ; 

purple = plasma serine protease inhibitor {P05154) ; 
turquoise = al anti-chymotrypsin (POlOll) ; mauve = 40S \ 
ribosomal protein (P08865) . Note that chip sequence 8 and 
12 did not sequence verify. 

35 
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EXAMPLE 4 

Genome-Derived Single Exon Probes Useful For Measuring 
Human Gene Expression 

5 

The protocols set forth in Examples 1 and 2, 
supra, were applied to additional human genomic sequence as 
it became newly available in GenBank to identify unique 
exons in the human genome that could be shown to be 
10 expressed at significant levels in Fetal liver tissue. 

These unique exons are within longer probe 
sequences. Each probe was completely sequenced on both 
strands prior to its use on a genome -derived single exon 
microarray; sequencing confirms the exact chemical 

15 structure of each probe. An added benefit of sequencing is 
that it placed us in possession of a set of single base- 
incremented fragments of the sequenced nucleic acid, 
starting from the sequencing primer 3' OH, (Since the 
single exon probes were first obtained by PGR amplification 

20 from genomic DNA, we were of course additionally in 

possession of an even larger set of single base incremented 
fragments of each of the 12,673 single exon probes, each 
fragment corresponding to an extension product from one of 
the two amplification primers.) 

25 The structures of the 12,673 unique single exon 

probes are clearly presented in the Sequence Listing as SEQ 
ID Nos.: 1 - 12,673. The 16 nt 5 ' primer sequence and 16 
nt 3 ' primer sequence present on the amplicon are not 
included in the sequence listing. The sequences of the 

30 exons present within each of these probes is presented in 
the Sequence Listing as SEQ ID Nos . : 12,674 - 25,129, 
respectively. It will be noted that some amplicons have 
more than one exon, some exons are contained in more than 
one amplicon. 

35 As detailed in Example 2, expression was 
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demonstrated by disposing the amplicons as single exon 
probes on nucleic acid microarrays and then performing two- 
color fluorescent hybridization analysis; significant 
expression is based on a statistical confidence that the 
5 signal is significantly greater than negative biological 
control spots. The negative biological control is formed 
from spotted DNA sequences from a different species. Here, 
32 sequences from E.Coli were spotted in duplicate to give 
a total of 64 spots. 
10 For each hybridisation (each slide, each colour) 

the median value of the signal from all of the spots is 
determined. The normalised signal value is the arithmetic 
mean of the signal from duplicate spots divided by the 
population median, 
15 Control spots are eliminated if there is more 

that a five -fold difference between each one of the 
duplicate spots raw signals. 

The median of the signal from the remaining 
control spots is calculated and all subsequent calculations 
20 are done with normalised signals. 

Control spots having a signal of greater than 
median + 2,4 (the value 2.4 is roughly 12 times the 
observed standard deviation of control spot populations) 
are eliminated. Spots with such high signals are considered 
25 to be "outliers". 

The mean and standard deviation of the modified 
control spot populations are calculated. 

The mean + 3x the standard deviation (mean + 
(3*SD)) is used as the signal threshold qualifier for that 
30 particular hybridisation. Thus, individual thresholds are 
determined for each channel and each hybridisation. 

This means that, assuming that the data is 
distributed normally, there is a 99% confidence that, any 
signal exceeding the threshold is significant. 
35 The probes and their expression data are 
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presented in Table 4, set forth respectively in Example 5. 
Example 5 presents the subset of probes that is 
■significantly expressed in the human Fetal liver and thus 
presents the subset of probes that was recognized to be 

5 useful for measuring expression of their cognate genes in 
human Fetal liver tissue. 

The sequence of each of the exon probes 
identified by SEQ ID NOS . : 12,674 to 25,129 was 
individually used as a BLAST {or, for SWISSPROT, BLASTX) 

10 query to identify the most similar sequence in each of 
dbEST, SwissProt (BLASTX)', and NR divisions of GenBank . 
Because the query sequences are themselves derived from 
genomic sequence in GenBank, only nongenomic hits from NR 
were scored. 

15 The smallest in value of the BLAST {or BLASTX) 

expect {"E") scores for each query sequence across the 
three database divisions was used as a measure of the 
"expression novelty" of the probe's ORF. Table 4 is sorted 
in descending order based on this measure, reported as 

20 "Most Similar {top) Hit BLAST. E Value". Those sequences for 
which no "Hit E Value" is listed are those exons which were 
found to have no similar sequences. 

As sorted. Table 4 thus lists its respective 
probes {by "AMPLICON SEQ ID NO.:" and additionally by the 

25 SEQ ID NO:, of the exon contained within the probe: "EXON 
SEQ ID NO.:") from least similar to sequences known to be 
expressed {i.e., highest BLAST E value), at the beginning 
of the table, to most similar to sequences known to be 
expressed {i.e., lowest BLAST E value), at the bottom of 

30 the table. 

Table 4 further provides, for each listed probe, 
the accession number of the database sequence that yielded 
the "Most Similar {top) Hit BLAST E Value", along with the 
name of the database in which the database sequence is 
35 found ("Top Hit Database Source"). 
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Table 4 further provides SEQ ID NOS . 
corresponding to the predicted amino acid sequences where 
they have been determined for the probe and exon nucleotide 
sequences. These are set out as PEPTIDE SEQ ID NOS.:. The 
5 peptide sequences for a given exon are predicted as 

follows: Since each chip exon is a consensus sequence drawn 
from predictions from various exon finding programs (i.e. 
Grail, GeneFinder and GenScan) , the multiple initial ORFs 
are first determined in a uniform way according to each 

10 prediction. In particular, the reading frame for predicting 
the first amino acid in the peptide sequence always starts 
with the first base of any codon and ends with the last 
base of non- termination codon. Next, for each strand of the 
exon, initial ORFs are merged into one or more final ORFs 

15 in an exhaustive process based on the following criteria: 
1) the merging ORFs must be overlapping, and 2) the merging 
ORFs must be in the same frame. 

The Sequence Listing, which is a superset of all 
of the data presented in Table 4, further includes, for 

20 each probe, the most similar hit, with accession number and 
BLAST E value, from the each of the three queried 
databases . 

Table 4 further lists, for each probe, a portion 
of the descriptor for the top hit ("Top Hit Descriptor") as 

25 provided in the sequence database. For those ORFs that are 
similar in sequence, but nonidentical to known sequences 
(e.g., those with BLAST E values between about le-05 and 
le-100) , the descriptor reveals the likely function of the 
protein encoded by the probe's ORF. 

30 Using BLAST E value cutoffs of le-05 (i.e., 1 x 

10'^) and le-100 (i.e., 1 x 10'^°°} as evidence of similarity 
to sequences known to be expressed is of course arbitrary: 
in Example 2, supra, a BLAST E value of le-3 0 was used as 
the boundary when only two classes were to be defined for 

35 analysis (unknown, >le-30; known <le-30) (see also FIG. 8) . 
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Furthermore, even when the "Most Similar (Top) Hit BLAST E 
Value" is low, e.g., less than about le-100 - which is 
probative evidence that the query sequence has previously 
been shown to be expressed - the top hit is highly unlikely 
5 exactly to match the probe sequence. 

First, such expression entries typically will not 
have the intronic and/or intergenic sequence present within 
the single exon probes listed in the Table. Second, even 
the ORF itself is unlikely in such cases to be present 

10 identically in the databases, since most of the EST and 
mRNA clones in existing -databases include multiple exons, 
without any indication of the location of exon boundaries. 

As noted, the data presented in Table 4 represent 
a proper subset of the data present within the attached 

15 sequence listing. For each amplicon probe (SEQ ID NOs . : 1 
- 12,673) and probe exon (SEQ ID NOs.: 12,674 - 25,129, 
respectively) , the sequence listing further provides, 
through iterated annotation fields <220> and <223>: 

(a) the accession number of the BAG from which 
20 the sequence was derived ("MAP TO")/ thus providing a link 

to the chromosomal map location and other information about 
the genomic milieu of the probe sequence;- 

(b) the most similar sequence provided by BLAST 
query of the EST database, with accession number and BLAST 

25 E value for the "hit"; 

(c) the most similar sequence provided by BLAST 
query of the GenBank NR database, with accession number and 
BLAST E value for the "hit"; and 

(d) the most similar sequence provided by BLASTX 
30 query of the SWISSPROT database, with accession number and 

BLAST E value for the "hit". 



EXAMPLE 5 

35 Genome -Derived Single Exon Probes Useful For Measuring 
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Table 4 {526 pages) presents expression, homology, and 
functional information for the genome -derived single exon 
5 probes that are expressed significantly in human Fetal 
liver . 



\ 
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1. A spatially-addressable set of single exon nucleic acid 
probes for measuring gene expression in a sample derived 
5 from human Fetal liver comprising a plurality single exon 
nucleic probes, said probes comprising any one of the 
nucleotide sequences set out in SEQ ID NOs : 1 - 12,673 or a 
complementary sequence, or a portion of such a sequence. 

10 2. A spatially-addressable set of single exon nucleic acid 
probes as claimed in claim 1 wherein each of said plurality 
of probes is separately and addressably amplifiable. 

3 . A spatially-addressable set of single exon nucleic acid 
15 probes as claimed in claim 1 wherein each of said plurality 

of probes is separately and addressably isolatable from 
said plurality. 

4. A spatially-addressable set of single exon nucleic acid 
20 probes as claimed in any of claims 1 to 3 wherein said 

probes comprise any one of the nucleotide sequences set out 
in SEQ ID NOS.: 12,674 - 25,129. 

5. A spatially-addressable set of single exon nucleic acid 
25 probes as claimed in any of claims 1 to 4, wherein each of 

said plurality of probes is amplifiable using at least one 
common primer. 

6. A spatially-addressable set of single exon nucleic acid 
30 probes as claimed in any of claims 1 to 5 wherein the set 

comprises between 50 - 20,000 single exon nucleic acid 
probes . 

7. A spatially-addressable set of single exon nucleic acid 
35 probes as claimed in any of claims 1 to 6 , wherein the 
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average length of the single exon nucleic acid probes is 
between 200 and 500 bp. 



8. A spatially-addressable s'^t of single exon nucleic acid 

5 probes as claimed in any of claims 1 to 7, wherein at least 

50% of said single exon nucleic acid probes lack 

f 

prokaryotic and bacteriophage vector sequence. 

9. A spatially-addressable set of single exon nucleic 
10 probes as claimed in any of (;;laims 1 to 8 , wherein at 

50% of said single- exon nucleic acid probes lack 
homopolymeric stretches of A or T. 

10. A spatially-addressable set of single exon nucleic acid 
15 probes as claimed in any of claims 1-9 characterised in 

that said set of probes is addressably disposed upon a 
substrate . 

11. A spatially-addressable set of single exon nucleic acid 
20 probes as claimed in claim 10 wherein said substrate is 

selected from glass, amorphous silicon, crystalline silicon 
and plastic. 

12. A microarray comprising a spatially addressable set of 
25 single exon nucleic acid probes as claimed in any of claims 

1-11. 

13 . A single exon nucleic acid probe for measuring human 
gene expression in a sample derived from human Fetal liver 

30 comprising a nucleotide sequence as set out in any of SEQ 
ID NOs . : 1 - 12,673 or a complementary sequence or a 
fragment thereof wherein said probe hybridizes at high 
stringency to a nucleic acid -molecule expressed in the 
human Fetal liver. 



acid 
least 
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14. A single exon nucleic acid probe as claimed in claim 13 
comprising a nucleotide sequence as set out in any of SEQ 
ID NOs . : 12,674 - 25,129 or a complementary sequence or a 
fragment thereof . 

5 

15. A single exon nucleic acid probe for measuring human 
gene expression in a sample derived from human Fetal liver 
which is a nucleic acid molecule having a sequence encoding 
a peptide comprising a peptide sequence as set out in any 

10 of SEQ ID NOs.: 25,130 - 37,156, or a complementary 
sequence or a fragment thereof wherein said probe 
hybridizes at high stringency to a nucleic acid expressed 
in the human Fetal liver. 

15 16. A single exon nucleic acid probe as claimed in any one 
of claims 13 to 15 wherein said single exon nucleic acid 
probe comprises between 15 and 25 contiguous nucleotides of 
said SEQ ID NO. 

20 17. A single exon nucleic acid probe as claimed in any one 
of claims 13 to 15, wherein said probe is between 3 - 25 kb 
in length. 

18. A single exon nucleic acid probe as claimed in any one 
25 of claims 13 - 17, wherein said probe is DNA, RNA or PNA. 

19. A single exon nucleic acid probe as claimed in any one 
of claims 13 - 18, wherein said probe is detectably 
labeled . 

30 

20. A single exon nucleic acid probe as claimed in any one 
of claims 13 - 19, wherein said probe lacks prokaryotic and 
bacteriophage vector sequence . 

35 21 . A single exon nucleic acid probe as claimed in any one 
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of claims 13 - 20, wherein said probe lacks homopolymeric 
stretches of A or T. 



22. A method of measuring gene expression in a sample 
5 derived from human Fetal liver, comprising: 

contacting the microarray of claim 12, with a first 
collection of detectably labeled nucleic acids, 
said first collection of nucleic acids derived 
from mRNA of human Fetal liver; and then 
10 measuring the label detectably bound to each probe of 

said microarray. 



23. A method of identifying exons in a eukaryotic genome, 
comprising ; 

15 algorithmically predicting at least one exon from 

genomic sequence of said eukaryote; and then 
detecting specific hybridization of detectably labeled 
nucleic acids to a single exon probe, 
wherein said detectably labeled nucleic acids are derived 
20 from mRNA from the Fetal liver of said eukaryote, said 

probe is a single exon probe having a fragment identical in 
sequence to., or complementary in sequence to, said 
predicted exon, said probe is included within a microarray 
according to claim 12, and said fragment is selectively 
25 hybridizable at high stringency. 



24. A method of assigning exons to a single gene, 
comprising : 

identifying a plurality of exons from genomic 
30 sequence according to the method of claim 23; and 

then 

measuring the expression of each of said exons in a 
plurality of tissues and/or cell types using 
hybridization to single exon microarrays having a 
35 probe with said exon, 
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wherein a common pattern of expression of said exons in 
said plurality of tissues and/or cell types indicates that 
the exons should be assigned to a single gene. 

5 25. A nucleic acid sequence as set out in any of SEQ ID 
NOs : 1 - 25,129 which encodes a peptide. 

26. A peptide encoded by a sequence as set out in any of 
SEQ ID Nos: 1 - 25, 129. 

10 

27. A peptide comprising a sequence as set out in any of 
SEQ ID Nos: 25,130 - 37,156. 



102 



wo 01/57277 



PCT/USn 1/00669 



£ S s 



E I h- o 

IIP 



O 

(/I y 



J liJ 2 



£ iiJ 2 



9 S 

3 5 



1/S26 



wo (H/57277 



PCT/US« 1/00669 



I 

1.1 

f 




2/526 



wo 01/57277 



PCT/USOl/00669 




3/526 



wo 01/57277 



PCT/USOt/0066'J 



o 



E TS. 




i ^ I 



o 

LU n 
(£ Q 

o 



ill 



4/526 



wo 01/57277 



PCT/USOl/00669 























o 




































o 
a 

1 

li 
^% 
it 

^£ 
















o 

CO 

d 
o 

0 

i 

i 

a 
O) 

e" 

E 

a 

a 
c 

1 

< 
z 
o 

«> 

1 
1 

E 

M 

s 












N; LIPOVITELLIN LV- 














CATALYTIC SUBUNII 






0 gene, partial cds 








in 


if> 
a 


complete cds 










Top Hit Descriptor 


DNA POLYMERASE III, ALPHA CHAIN POLC-TYPE (POLIII) 


LYCOPENE BETA CYCLASE 


HEAT-INDUCIBLE TRANSCRIPTION REPRESSOR HRCA 


PNEUMOLYSIN (THIOL-ACTIVATED CYTOLYSIN) 


Arabtdopsis thaliana DNA chromosome 4, contig fragmerrt No. 67 


Chicken alpha-cardlac actin gene 


Chicken alpha-cardiac acUn gene 


NEL PROTEIN PRECURSOR (NEL-RELATED PROTEIN 2) 


VITELLOGENIN PRECURSOR (VTG) [CONTAINS: LIPOVITELLIN LV-' 
LIPOVlTEaiN LV.2] 


REP1 PROTEIN 


REP1 PROTEIN 


RHODOPSIN 


RHODOPSIN 


Bovine immunodenciency-like virus surface envelope gene, 5' end of cds 


HOMEOBOX PROTEIN CEH-20 


DNA POLYMERASE GAMMA (MITOCHONDRIAL DNA POLYMERASE 


Homo sapiens HERPUD1 gene for stress protein Herp, completo cds 


QV4+IT0691-270400-186-f09 HT0691 Homo sapiens cDNA 


Drosophlla orientacea RIB retrotransposaUe element reverse transcriptas 


HYPOTHETICAL 61.1 KD PROTEIN C23E2.03C IN CHROMOSOME 1 


RHODOPSIN 


COLICIN N IMMUNITY PROTEIN (MICROCIN N IMMUNITY PROTEIN] 


601894910F1 NIH_MGC_19 Homo sapiens cDNAckwro IMAGE:412411' 


8 

S 

lii 
O 

1 

D 

9 

M 
g 

t 

O 

g 
I 

o' 

O 

2 
I 
z 

S 


Cants familiaris skeletal muscle chloride channel CIC-1 (CLCNI) mRNA, < 


Mycobactsrium tuberculosis H37Rv complete genome; segment 103/162 


Human hereditary haamochromatosie regton, histone 2A-1ike protein gene. 
(HLA-H) gene, RoRet gene, and sodium phosphate transporter (NPT3) ge 


lEunlce australis histone H3 (H3) gene, partial ods 


IRC3-GN0042-1 00800-011^:10 GN0042 Homo sapiens cDNA 


|PMO-BT0547-310100-002-b04 BT0547 Homo sapiens cDNA 


— So 
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lOT 




1- 

O 
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Toph 
DatatM 
Sourc 


q: 
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or 
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SWISSPF 


a: 

0. 


Q. 




q: 






Q. 










Q. 


0. 




















SSIMS 


SWISS 


i 

(0 


NT 


SSI MS 


NT 


NT 


NT 


SSI MS 


5SIMS 


JSIMS 


swisi 


SWISS 


NT 


swisi 


iSIMS 


NT 


ESTY 


NT 


SWISi 


SWlSj 


SWIS! 


EST J 


EST ^ 


MT 




IN 


NT 


ESTJ 


EST > 


Top Hit Acession 
No 


P75080 1 


Q55276 1 


P47447 1 


AF1 75425.1 | 


P1 1990 1 


AL161571.2 1 


X02212.1 ] 


X02212.1 1 


•n 
O 


091 062 


P40379 1 


? 
a. 


01 7094 1 


01 7094 1 


L43126.1 1 


P41779 1 


P54098 1 


AB034990.1 1 


! 


s 
1 

Si 
< 


Q10136 1 


016005 1 


§ 


n 

z 

o 

n 
u. 
m 


i 

m 


CM 

in 

i 

< 


Z63660.1 


U91328.1 


AF185255.1 


BF367909.1 


IAW750067.1 


MostSlrrllar 
(Top) Hit 
BLAST E 
Value 


o 
o 
+ 

UJ 

la 
in 


o 
o 

vi 


o 
o 

•f 

tei 


o 
o 

IO 


o 
o 

+ 


5.5E+00I 


5.4E+00I 


5.4E+001 


5.4E+M] 


S.4E+00 


5.4E+00I 


o 
? 

UJ 
•ri 


S.4E+0QI 


5.4E+001 


5.3E+O0I 


o 
UJ 

n 


5.3E+00I 


5.3E+00 


vi 


5.2E+00I 


? 

tri 


S.IE+OOl 


5.1E+00 


5.OE>O0i 


5.0E+00, 


♦ 

UJ 
o 

iri 


5.0E+00 


f 

UJ 

<j> 


8 
+ 


1 4.8E+00 


1 4.8E+Q0 


Expression 
Signal 


1.31^ 


s 


0.73| 


1.54 


3.79| 


2.08} 






CM 

d 


1.58 


0.78| 


0.78| 


1.33 


1.33 


1,52 


0.671 


3.71 


0.62 


1,04| 


d 




0.88| 


0.97 


0.85 


1 0.691 


3.37 


1 13.54 


0.71 


1 12.06 


0.65 


1 


ORF SEQ 
ID NO: 




360591 


317801 




360581 




323261 


32327] 






341931 


341941 


354231 


35424 


29930 














34354 


35217 


31817 




35817 


36723 






1 33547 




Exon 
SEQ ID 
NO: 


200871 


230471 


190021 


23190 


23046] 


24024] 


19507] 


19507 


19895] 


i 20354 


21273 


212731 


i 


22444 


17474 


19211 


20568 


21452 


18282 


22769 


23586 


21430 


22237 


19033 


22591 


22821 


23677 


22626 


16727 


20636 


21017 


Probe 
SEQ ID 

NO: 


I 75701 


1 113491 


1 63991 


10658 


1 113481 


11578 


1 70091 


70O9 


7369 


7811 


8734 


8734 


9949 


9949 


4899 


6614 


8024 


r-^ 


S655 


10274 


11074 


8892 


9739 


6430' 


10096 


10327 


11170 


10131 


4135 


8095 


1 6478 
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Top Hit Descfipta 


BRAIN-SPEOFIC ANGIOGENEStS INHIBITOR 1 PRECURSOR | 


ADHERENCE FACTOR (ADHESION AND AGGREGATION MEDIATING SURFACE ANTIGEN) | 


STRUCTURAL POLYPROTEIN [CONTAINS: MAJOR STRUCTURAL PROTEIN VP2; 
NONSTRUCTURAL PROTHIN VP4; MINOR STRUCTURAL PROTEIN VP3] 


STRUCTURAL POLYPROTEIN [CONTAINS: MAJOR STRUCTURAL PROTEIN VP2: 
NONSTRUCTURAL PROTEIN VP4: MINOR STRUCTURAL PROTEIN VP3] 


60201 7413F1 NCI_CGAP_Brn64 Homo sapiens cON A clone IMAGE:41 53059 5* | 


Buxus harlandil maturase K (matK) gene, partial cds; chloroplast gene for chloroplast product J 


Arabidopsis thaliana DNA chromosome 4, contig fragment No. 52 | 


Mus musculus endomucin (LOCS3423), mRNA | 


601 342758F1 NIH_MGC_53 Homo sapiens cDNA clone IMAGE:3684807 5' | 


Mus musculus endomucin (LOC53423), mRNA | 


< 
Z 

i 

CO 

o 

c 
o 
w 

1 
1 

v> 

3 

» 
3 

E 
n 

1 


Mus musculus per-hexamer repeat gene 3 (PhxrS). mRNA | 


Homo sapiens apoA pdymorphism Kringle IV gene. e)»ns 1 and 2 j 


Ipomoea purpurea chaloone synthase (CHSB) gene Including complete 5'UTR and complete cds j 


Botrytis cinerea strain T4 cDNA library under condition s of nitrogen deprivation | 


THYMOSIN BETA-4 (HUMAN); 


CMO-BT0281-031199-087+04 BT0281 Homo sapiens cDNA 1 


Mus musculus sphingosine Kinase (SPHK1 b) mRNA, complete cds | 


Mus musculus SRY-bac containing gene 1 3 (Soicl 3), mRNA | 


< 

z 
o; 

E 

fO 

i 

w 

to 
o 

1 

g 

"c 

i 

'>- 

U) 

M 
3 
3 
^ 
3 

E 

m 

3 


Mycobacterium fortuitum furA II gene ] 


Mus musculus SH2-containing inositol 5-phosphatase (Ship) gene, exons 16 through 27. and complete cds 


faba bean necrotic yellows virus C2-Eg gene, isolate Egyptian EV V93 | 


faba bean necrotic yetlows wrus C2-Eg gene, isolate Egyptian EV1-93 | 


Arabidopsis thallane DNA chromosome 4. contig fr^ment No. 40 | 


Mus musculus cleavage and polyadenytet'on specificity factor 3 (Cpsf3), mRNA ] 


Hantavirus Z10 segment M G1/G2 glycoprotein (Z10) gene, complete cds j 


Homo sapiens ATP-bindtng cassette, sub-family B (MDRfTAP), member 4 (ABCB4), mRNA ) 


Aspergillus nldulans recQ gene for DNA helicase, exons 1 -4 | 


Aspergilus nidulans recQ gene for ONA helicase, exons 1 -4 | 


Top Hit 
Database 
Source 


SWISSPROT 1 


SWISSPROT 1 


SWISSPROT 


SWISSPROT 


EST HUMAN ] 


NT i 


NT 1 


NT 1 


EST HUMAN ] 


NT ] 


NT ] 


NT 1 


NT ] 


NT ] 


NT 1 


EST HUMAN 


EST HUMAN ] 


NT 1 


NT ] 


NT 1 


NT i 


NT 


NT ^ 


NT 1 


NT 


NT 1 


NT 


NT 1 


NT 


NT 


Top Hit Acesston 
No. 


014514 1 


§ 

<D 
"f 
CL 


P05844 


P05844 


BF344171.1 1 


AF1 86398.1 ] 


AL1 61 552.2 ] 


8393724] 


i 


8393724] 


6679306] 


e67S3oe| 


LI 4005.1 1 


U15947.1 1 


AL1 16458.1 1 


AW088191.1 


BE063527.1 ] 


AF068749.1 ] 


6755601] 


6755601] 


Y1 7062.1 ] 


i 

to 

< 


AJ1 32180.1 ] 


AJ1 32180.1 ] 


AL161540.2 1 


9055193] 


to 

& 

<o 

< 


114192201 


AJ271 844.1 j 


AJ271 844.1 1 


Most Similar 
(Top) Hit 
BUVSTE 
Value 


o 


O 


2.9E+00 


2.9E+00 


2.9E+O0I 


o 
? 




f 

oi 


2.8E-K)0] 


2.8E+00] 


2.7E+00] 


f 


2.7E+00] 


2.7E+00] 


C^i 


2.7E+00 


s 

+ 

UJ 

K 
oj 


2.6E400] 


2.6E+00] 


f 


2,6E+O0j 


f 

LU 
<D 

ci 


2.6E+00] 


2.6E+00I 


f 

UJ 
(D 

ci 


s 

UJ 
(D 

c\i 


8 

♦ 

UJ 
CD 
M 


2.6E+O0I 


2.5E+00 


f 

CNi 


Expression 
Signal 


n 


a> 


<D 

d 


to 
d 


d 






4.88] 


to 
d 


1.68] 


ca 


CO 

a> 




CO 

d 


1.681 


0.63 


? 


4.97] 


1.04] 


1.94] 


2.42] 


5.98 


1.08] 


1.08] 


3.02] 


lO 




2.781 


2.29 




ORFSEQ 
ID NO: 


32645] 


328661 


33260 


33261 


33490] 


26634] 




32740] 




32740] 


25394] 


25395] 


31073] 






33191 




29812] 


31068] 


31069] 






33447] 


33448] 


a 

s 

to 
to 








o 

i 


26641 


Exon 
SEQ ID 

NO; 


197901 


20001 1 


20352 


20352 


205831 


1 


14267 


s 


22031] 


19874] 


12911] 


12911] 


16366] 


20629] 


21436] 


20292 


22888 


17362] 


18362] 


18382 


18612 




20545] 


20545 


22076] 


22752] 




24986] 


141051 


141051 


Probe 
SEQ ID 

NO: 


7262 


7479 


7809 


7809 


8041 


1504 


1 16751 


7348 


9531 


I 10569] 


■n 
tM 


!n 
tM 


1 5740 


8088 


s 
s 


9353 


] 10394] 


i 4781] 


5736 


I S736J 


(M 

s 


7689 


8003 


8003 


9576 


1 10257] 


] 10907] 


] 12390] 


] 15131 


1 1513| 
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Top Hit DoscriptDT 


TEICHOIC ACID BIOSYNTHESIS PROTEIN F | 


TEICHOIC ACID BIOSYNTHESIS PROTEIN F | 


TEICHOIC ACID BIOSYNTHESIS PROTEIN F I 


TEICHOIC ACD BIOSYNTHESIS PROTEIN F | 


Vibrio chderae cbcA gene and cb(B gene for cholera toxins, complete cds | 


QV4-FTO00M 1 0500-205-g07 FT0005 Homo sapiens cDNA | 


Rice DNA for aldolase Cl , complete cds | 


601175779F1 NIH_MGC_1 7 Homo sapiens cDNA clone IMAGE:3531 090 5' | 


Mus musculus EIF4H gene, partial cds; LIMK1 gene, complete cds; and ELN gene, partial cds | 


Chicken alphas collagen type VI mRNA, 3' end 1 


Homo sapiens double C2-llke domains, alpha pOC2A) mRNA | 


VITELLOGENIN 1 PRECURSOR (YOLK PROTEIN 1 ) | 


CD27L RECEPTOR PRECURSOR (T-CELL ACTIVATION ANTIGEN CD27) (T14) | 


CD27L RECEPTOR PRECURSOR (T-CELL ACTIVATION ANTIGEN CD27) (T14) | 


Helicobacter pylori, strain J99 section 47 of 132 of the complete genome | 


1 

! 

r 

0. 

v> 
o 

■<? 

§ 
i 

1 

en 


ENDOCHITINASE B PRECURSOR (CHN^B) | 


SKIN GRANULE PROTEIN PRECURSOR | 


SKIN GRANULE ^kOTEIN PRECURSOR | 


H. sapiens CTGF gene and promoter region 1 


XYLULOSE KINASE (XYLULOKINASE) | 


CO 

n 
n 

n 

ui 

1 
s 

u 
< 
z 

1 

o 

% 

T 

^1 
a. 
< 

8, 

O 

z 

1 


in 

0 

CO 
CO 

« 
uJ 

i 

< 

Z 

o 

1 

3 

T 

i 

a 

o 

o 
z 

i 

JZ 


DEN ITR IFi CATION REGULATORY PROTEIN NIRQ | 


Badllus subtilis chromosomal DNA, reg|on 75 degrees: gIpPFKD operon and downstream 1 


Fragaria x ananassa cytosdic ascorbate peroxidase (ApxSC) gene, ApxSC-c allele, complete cds 


G.domesticus artificial single chain antibody gene (L3) 1 


Bos taurus partial cytb gene for cytochrome b [ 


J7340F Human fetal heart, Lambda ZAP Express Homo sapiens cDNA clone J7340 5' similar to 
PROLYLCARBOXYPEPTIDASE 


Rattus norveglcus ATPase, Ca++ transporting, ubiquitous (Atp2e3), mRNA ) 


MAJOR CENTROMERE AUTOANTIGEN B (CENTROMERE PROTEIN B) (CENP-B) \ 


|M.maz8l dnaK and dnaJ genes homologues coding for DnaK and DnaJ 1 


Top Hit 
Database 
Source 


SWISSPROT I 


SWISSPROT 1 


SWISSPROT 1 


SWISSPROT 1 


Z 


EST_HUMAN 1 


1- 
z 


EST HUMAN | 




1- 
2 


t- 
z 


SWISSPROT 1 


SWISSPROT 1 


SWISSPROT 1 


1- 
z 


EST_HUMAN 1 


SWISSPROT J 


SWISSPROT 1 


SWISSPROT 1 


NT 1 


SWISSPROT 1 


EST HUMAN ' 


EST HUMAN 


SWISSPROT 


1- 
z 


^ 


1- 
z 


1- 

2 


2 

1 

U 
I 

UJ 


t- 
z 


ISWISSPROT 


K 
Z 


Top Hit Acession 
No. 


P13485 1 


P13485 1 


PI 3485 1 


IS 

1 

CL 


030052. 1 1 


Si 

a 
1 


D50307.1 1 


BE297758.1 | 


AF289665.1 1 


M24282.1 1 


45033521 


P02843 I 


P26842 1 


i 


AE001486.1 ] 


AW 8751 26.1 1 


P24091 J 


PI 3673 1 


P13673 1 


X92511.1 1 


§ 

CL 


BE326702.1 1 


BE326702.1 1 


Q51481 


Y14079.1 


AF1 58652.2 


246724.1 


1AJ401081.1 


N86245.1 


1 6978554 


]P07199 


]X60265.1 


E X K § 


2.5E+O0I 


o 
? 

IH 


2.5E+00I 


2.5E+00I 


o 
? 


2.5E+O0I 


2.5E+O0| 


? 


2.5E+O0I 


f 


2.4E+001 


2.4E+001 


2.4E+00I 


o 

? 

liJ 


2.4E+00I 


2.4E+00I 


o 
o 

*• 

UJ 


2.4E+00I 


? 

UJ 
<N 


2.4E+00I 


2.4E+0^ 


§ 

<N 


§ 
*• 

UJ 

N 


? 

UJ 


2.4E+00 


2.4E+00 


2.3E+00 


2.3E+00 


2.3E+00 


f 

UJ 


1 2.3E+00 


1 2.3E+00 


Expression 
Signal 


1.71 


1.71 


1.39| 


1.39| 


0.731 


1.051 


1.751 


0.88| 


! 


o 


6.761 


4.02) 


! 


1.99 


2.33] 


1.62 


8.16] 


2.59 


2.59 








1.63 


1.27 


1.69 


2.27 


11.15 


1.65 


0.91 


2.54 


1 4.61 


1 1.06 


ORF SEQ 
ID NO: 


313341 


313351 


313341 


313351 


322341 


331351 


s 

J 


352471 




28144 


300521 


315361 


33536 


33537] 






34222 


354271 


354281 


355041 




O 
lO 


357111 


359861 


3649^ 


36826 


1 26413 






1 

TO 


\ 


] 20264] 33159 


Exon 
SEQ ID 

NO: 


186011 


186011 


18601 


18601] 


CO 

1 


202441 


21569] 


22264] 


24131] 


15663] 


17607] 


18774] 


20624 


208241 


206941 


211241 


21301 


22446 


224461 


22512 


22636 


22720 


22720 


22977 


23473 


23768 


13890 


16788 


18820 


1 19999 


1 25120 


Probe 
SEQID 
NO: 


5981 


5981 


6586 


1 65861 


6828 


7736 




i 


1 11724] 


3047 


1 5033] 


6161 


8082 


1 80821 


1 81531 


8565 


8762 


1 99511 


tS66 


10017 


1 10141! 


] 102251 


1 10225 


] 10483 


] 10958 


11237 


1296 


4199 


6000 


7477 


7993 


1 7756 
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TOR 
OTEIN 
1 LIGAND- 


TOR 
OTEtN 
1 LIGAND- 


















;6 Mouse 


S 

1 






893965 3' 


Top Hit Descriptor 


Poi^terus ornatlpinnis mitochor^rion, complete genome 


ALPHA-{1.3)-FUCOSYLTRANSFERASE (GALACTOSIDE 3-L-FUCOSYLTRANSFERASE) 
(FUCOSYLTRANSFERASE 4) (FUCT-IV) 


ANNEXINVII(SYNEXIN) 


602069121F1 NIH_MGC_58 Homo sapiens cDNA clone 1 MAGE :40681 73 5' 


6020691 21 F1 N1H_MGC_58 Homo sapiens cDNA clone IMAGE:4068173 5' 


601433673F1 NIH_MGC_72 Homo sapiens cDNA clone IMAGE:3918643 5*. 


Magnaporthe grtsea Qass IV chitin synthase (chs4) gene complete cds 


Rat gene for regucalcin. exonl (non-coding exon) 


Rat gene for regucalcin, exonl (non-coding exon) 


SORTILIN-REUTED RECEPTOR PRECURSOR (SORTING PROTEIN-RELATED RECEP 
CONTAINING LDLR CLASS A REPEATS) (MSORLA) (SORLA-1) (LOW -DENSITY LIPOPR 
RECEPTOR RELATIVE WITH 11 LIGAND-BINDING REPEATS) (LDLR RELATIVE WITH 1 
BINDING REPEATS) (LR11)(> 


SORTILIN-RELATED RECEPTOR PRECURSOR (SORTING PROTEIN-RELATED RECEP 
CONTAINING LDLR CLASS A REPEATS) (MSORLA) {SORLA-1) (LOW -DENSITY LIPOPR 
RECEPTOR RELATIVE WITH 11 LIGAND-BINDING REPEATS) (LDLR RELATIVE WITH 1 
BINDING REPEATS) (LRU )(> 


< 
Z 

o 

o 
» 

t 
% 

X 

OI 

o 
& 

CD 
f 

9 

§ 

o 
o 

I 

01 


< 

z 
o 

o 
» 

i 
1 

X 

OJ 

o 

& 

CD 

9 

§ 

8 

i 

1- 
O 

i5 

01 


1600943401 TI NIH_MGC_17 Homo sapiens cDNA clone tMAGE:2959777 3" 


[MINOR VIRION STRUCTURAL PROTEIN MU-2 


[INSULIN-LIKE GROWTH FACTOR II PRECURSOR (IGF-II) (SOMATOMEDIN A) 


[nl95b02,s1 NCI CGAP_Co10 Homo sapiens cDNA clone IMAGE:1058379 3* 


o 

i 

ui 

i 
i 

o 
< 
z 
o 

i 

t 

X 

n 

u 

1 
1 
§ 

w 
r 
■** 

OI 


in 

i 

ill 
O 

i 

o 
< 
Z 

1 

I 
1, 
00 

Ll_ 

X 

z 

a 

1 

i' 

a) 
o 


bb17h12j<1 NIH_MGC_21 Homo sapiens cDNA clone IMAGE:2963207 3* similar to ob:D4583 


bbl 7h12.)Cl NIH_MGC_21 Homo sapiens cDNA clone IMAGE:2963207 3' similar to gb:D4583 
mRNA for nuclear poro-targeting-compIeK component o( (MOUSE); 


[601594733F1 NIH_MGC_9 Horrw sapiens cDNA clone IMAGE:3948561 5' 


1 TRANSPOSON TY1 PROTEiN A 


qm69b03.x1 Scares .placenta 8tD9vreek3 2NbHP8to9W Homo sapiens cDNA clone IMAGEil 
similar to gb;Y00433 GLUTATHIONE PEROXIDASE (HUMAN); 


Top Hit 
Database 
Source 


NT 1 


SWISSPROT 


SWISSPROT 1 


EST^HUMAN | 


EST_HUMAN | 


EST_HUMAN | 


NT 




NT J 


SWISSPROT 


SWISSPROT 


EST_HUMAN | 


EST HUMAN | 


EST_HUMAN 1 


ISWISSPROT 1 


[SWISSPROT 1 


[EST.HUMAN | 


[EST_HUMAN ] 


[EST HUMAN | 


z 

i 

X 

i5 

UJ 


EST HUMAN 


z 

i 

D 
X 

& 

UJ 


ISWISSPROT 


EST_HUMAN 


1 


58353171 
















































Top Hit Ace 
No. 


Q11127 


Q07076 


BF541 987.1 


BF541 987.1 


i 

m 


AF020528.1 


D67071.1 


D67071.1 


088307 


OB8307 


IBE927220.1 


IBE927220.1 


1 BE 250383.1 


8 


IP51459 


IAA594S74.1 


IAA137027.1 


|AA449012.' 


BE301560.' 


BE301560.' 


IBE741678.- 


IQ04706 


AI290373.1 


Most Similar 
(Top) Hit 
BLAST E 
Value 


2.3E+00I 


o 
o 
+ 

iiJ 

CO 


2.3E+00I 


2.3E+O0I 


2.3E+00| 


o 

? 
^ 

ci 


2.2E+001 


? 


? 

UJ 
OI 

ci 


2.2E+00 


2.2E+00 


2.2E+00I 


1 2.2E+00I 


1 2.2E+00| 


o 
o 
+ 

cvi 


[ 2.2E+00I 


o 
o 
+ 

UJ 
CM 
CM 


o 

UJ 
OJ 

oi 


o 
? 

w 

oi 


o 
? 

oi 


2.2E+00 


o 

9 
W 

rvi 


? 

CM 


2.2E+00 


Expression 
Signal 


0.53| 


1.79 




M 


, . 2.921 


1 7.3li 


1 91.071 


tr> 
■* 




12.27 


1227 


[ 0.951 


1 0.951 


O) 


« 


3,04 


3,58 


a> 
o 


1 25.231 


0.65 


0.65 


12-1^ 


2.57 


1.96 


ORF SEQ 
ID NO: 




34572 


362241 


371191 


371201 


310201 


291431 


294321 


o 


30591 


30592 


1 313731 


1 313741 


1 315931 


o 

CD 

« 


1 32107] 




1 327471 


CM 

i 


33494, 


in 

CO 






35443 


Exon 
SEQ ID 
NO: 


215751 


21633 


232131 


240551 


240551 


242781 


166851 


169881 




18177 


18177 


lA 

§ 


i 186351 


1 188221 


1 190961 


1 193031 


1 180571 


1 168841 


1 201151 


20588 


20588 


1 21791 


1 24793 


22461 


Probe 
SEQID 

NO:- 


§ 


r- 

<» 


1 106811 


1 116121 


1 116121 


1 119501 




1 44031 


1 44031 


5545 


5545 


6016 


1 60161 


6212 


1 64951 


1 67091 


7037 


1 73581 


1 76021 


8046 


8046 


1 9265, 
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UJ ^ i 

§ ' ■ 
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l' 



5 O ■ 



' E 

8 

TO 

T) 

g 



9 

Q. 
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£9 



|C!9 
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CO ^ 
£L g 
0. 



3 o 



3i 

0 ( 

01 i 
t. ! 
z : 

5^ i 
- 5 : 
. lij I 

W 



.1 



III 



H .O *S > 



11 



^ a o 
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u- xt 9 

° ® a 



5^ 



o'l 



1 1 



=,3 



i : 

M I 



f I? 



8- 



LU o 

"1 

cz: 9 
O 
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t 

x: 
o 


|Humanpapillomavinjs type 7 genomic DNA | 


Fugu rubripes neurofibromatosis type 1 (NF1 ). A-kinase anchor protein (AKAP84). BAW protein (BAW). and 
WSB1 protein (WSBI) genes, complete cds 


Fugu rubripos neurofibromatosis type 1 (NF1), A-kinase anchor protein (AKAP84), BAW protein (BAW), and 
WSB1 protein (WSBI) genes, complete cds 


1 Homo sapiens IUIad4 homdog (MAD4) mRNA t 


< 
z 
o 

u 
» 

s 
1 

I 
in 

8 

z 
z 

CD 

2 

o 

i 

8 

z 
z 

i 

o 


1 CMO-NN 1005-1 40300-286.h06 NN1005 Homo sapiens cDNA | 


16021 56687F1 NIH_MGC_83 Homo sapiens cDNA clone IMAGE:4297556 5' | 


i 

o 
in 

1 
1 

a 

6 

t 

o 

1 
I 

O 

z 

"5 

1 


iHomo sapiens mRNA fa KIAA1 157 protdn, partial cds | 


|DNA TOPOISOMERASE III ALPHA | 


]Homo sapiens mRNA for KIAA0905 protein, complete cds 1 


]SYNAPSIN II 1 


jSYNAPSIN II 1 


|Mu3 musculus WW domain binding protein 1 1 (Wbpll-pending). mRNA | 


|CM3-NN0006-300300-132-b12 NN0006 Homo sapiens cDNA j 


Homo sapiens caveolin-1/-2 locus. Contlgl . D7S522. genes CAV2 (exons 1. 2a. and 2b). CAVl (exons 1 and 
2) 


he23tD5.x1 NCI_CGAP_CML1 Homo sapiens cDNA done IMAGE:2919873 3' similar to contains Alu 


GLUCOAMYLASE PRECURSOR (GLUCAN 1.4-ALPHA-GLUCOSIDASE) (1.4-ALPHA-D-GLUCAN 
GLUCOHYDROLASE) 




]yg33f12.r1 Scares infant brain 1 NIB Homo sapiens cDNA clone IMAGE :34345 5' | 


]RCl-BT0313-301299-012-f05 BT0313 Homo sapiens cDNA | 




I602133135F1 NIH_MGC_81 Homo sapiens cDNA clone IMAGE:4288137 5' | 


]IL5-HT0198-29109&O08-C04 HT0198 Homo sapiens cDNA | 


< 
z 

Q 
u 

s 

1 

1 

X 
X 

Q 
O 

o 

1 

-i 


]Pandorina cdemaniae chloroplast rbcL gene for ribulose bisphosphate carboxylase, partial cds | 


[Pandorina cdemaniae chloroplast rbcL gene for ribulose bisphosphate earboxylaso, partial cds | 


2r36e09.n Soares_NhHMPu_Sl Homo sapiens cDNA done IMAGE:665512 5' similar to contains dement 


Top Hit 
Database 
Source 


iNT 1 


NT 


NT 


iNT 1 


iEST_HUMAN | 


|EST_HUMAN | 


|EST HUMAN ! 


lEST.HUMAN | 


NT 1 


[SWISSPROT 1 


iNT 1 


ISWISSPROT 1 


ISWISSPROT j 


iNT j 


|EST HUMAN j 


NT 


EST_HUMAN 


SWISSPROT 


IN 


lEST_HUMAN | 


z 

X 

1- 

(0 
UJ 


INT 1 


|EST_HUMAN | 


]EST_HUMAN | 


[EST HUMAN | 


iNT 1 


NT. 


EST.HUMAN 


Top Hit Acession 
No. 


1X74463.1 1 


CM 
< 


fM 

i 

< 


1 54537331 


IAW900455.1 1 


IAW900455.1 1 


|BFB81 547.1 | 


CD 
o 

1 


1AB032983.1 | 


IQ 13472 1 


]AB020712.1 1 


IQ92777 1 


IQ92777 1 


] 11096333] 


IAW893057.1 ] 


AJ 133269.1 


AW467760.1 


1 

£ 


IAJ271735.1 j 


]R20459.1 ] 


[BE064667.1 | 


|AF1 34844.1 | 


IBF575545.1 | 


|BE1 45374.1 1 


] BE 145374.1 | 




[D63441.1 j 


AA1 95528.1 


Most Similar 
(Top) Hit 
BLAST E 
Value 


1.4E+00| 






! 


1.4E+00I 


f 

lii 


t 1.4E+00| 


1 1.4E+00I 


f 


1.4E+00] 


1 1.4E+00] 


? 

UJ 


1.4E+00] 


I 1.4E+00] 


? 


1.4E+00 


f 


lE+OO 


lE+OO 


o 

o 

UJ 


IE+00 


? 

UJ 


8 
+ 

UJ 


IE+00 


8 

+ 

UJ 


o 
o 
♦ 
UJ 


o 
? 

UJ 


1.4E+00 






















Expression 
Signal 




2.79 


2.79 


1 0.68| 


1 1.14| 


1 1.14j 


i 1.78] 




1 5.04] 


1 2.73] 




] 2.32] 


] 2.32] 


1 0-67] 


1 0.77] 


2.31 




0.68 


] 4.01 ] 


] 2.13] 


] 3.72] 


] 0.56] 


1 0-77] 


] 0.67] 


] 0.671 






2.16 


ORF SEQ 
ID NO: 


27820| 


27922 


27923 




s 


293701 




306571 




1 31811] 




1 319361 


1 31937] 


1 31981] 


[ 32398] 


32720 


32735 






34487] 


[ 34603] 


CO 


35594] 


1 35637] 


35638] 


1 35912] 


1 35913] 


36485 


Exon 
SEQ ID 

NO: 


1 15250| 


15354 


15354 


1 159851 


1 16929] 




1 17267| 


1 18206] 


i 


] 19028) 


] 25116] 


1 19143] 


] 19143] 


1 19181] 


1 19570] 


i 


19870 


20818 


1 21268] 


1 21560] 


I 21660] 


1 


1 22604] 


22645 


] 22645] 


i 22912) 


j 22912] 


23463 


Probe 
SEQID 
NO: 




2802 


2602 


3376 


1 43421 


4342 


4685 


1 5575] 


1 5719] 


1 6425] 


tr> 


? 
•n 

(O 


6S« 


CD 


6911 


7330 


7343 


8277 


8729 


a 
s 


1 9125] 


9158 


10109 


10150 


10150 


10418 


10418 


CO 

1 
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O (D CL 













Emericella nidulans stsrigmatocyslin biosynthellc gene cluster. (stcA). (stcB), (stcC). (stcE). (aflR), (stcF). 
(stcl). (stcJ). (StcK). (StcL), (stcO). (stcQ). (stcS), (stcT). (stcU), (stcV) and (stcW) genes, complete cds 














































Top Hit Descriptor 




R.unicamis complete mitochondrial genome 


African swine fever virus, complete genome 


Drosophila metanogaster D-Titin gene, exons 1-37 


E.faecalis pbp5 gene 


Xylella fastJdiosa, section 15 of 229 of (he complete genome 


[Rattus norwegicus Aquaporin 4 (Aqp4). mRNA 


in 
us 

O 

i 

i 

< 
z 

c 
1 

1 
o 

I 

1 
§ 
§ 


|qd85c03.x1 Soares_testis_NHT Homo sapiens cDNA clone IMAGE:1736260 3' 


Homo sapiens solute carrier family 6 (neurotransmitter transporter), member 14 (SLC6A14), mRNA 


Macgregoria pulchra cytochrome b gene, complete cds; mitochondrial gene la mitochondrial product 


|ye89e03.r1 Soves fetal liver spleen 1NFLS Homo sapiens cDNA clone IMAGE:124924 5* 


|Mus musculus mRNA for ER protein 58 (EP58 gene) 


1 Maize mRNA for enotase (2-phosphD-D-glycBrate hydrolase) 




1 Herpes simplex virus type 1 (strain KOS) UL41 gene 


1 Arabidopsis thaliana DNA chromosome 4, contig fragment No. 84 


Mus musculus silent mating type information regulation 2. (S.cerevisiae, homolog>^ike(Sir2l). mRNA 


I602082582F1 NIH_MGC_61 Homo sapiens cDNA clone IMAGE:4246828 5' 


|lm39h11.x1 NCI_CGAP_Kid11 Homo sapiens cDNA clone IMAGE:21 60549 3' 


lAcetabularia calicutus mitochondrial COXI-like gene 


nt] 


to 
CO 

lij 

1 
1 

< 

i 

i 

X 

K 

0 
0 
2 

X 

z 

i 

s 


|Mus musculus mRNA for stretch responsive muscle (X^hromosome) protein (Srmx gene) 


m 

3 

8 
_j 

t 

3 

< 
z 
a 

a 

5 
1 

£ 

1 

< 


Yersinia pseudotuberculosis psaE, psaF, adhesin (psaA), cheperone (pseB), end usher (psaC) genes, 
complete cds 


Top Hit 
database 
Source 


















HUMAN 1 


_HUMAN 1 






_HUMAN 1 














HUMAN 


■z 

1 

I 






HUMAN 








u 


t- 
z 


NT 


NT 


NT 


NT 


NT 


NT 


NT 




EST 


NT 


NT 


EST 


NT 


NT 


NT 


NT 


NT 


NT 


EST 


EST 


NT 


NT 


EST 


NT 


NT 


NT 


g 




1 


















9739 
















7960 


















Top Hit Aces 
No. 




i! 


U18466.1 


IAJ27 1740.1 


U34740.1 




§ 
§ 


i 


1 

UJ 

m 


1 All 38 582.1 


? 


AF1 97861.1 


|R06037.1 


1AJ404004.1 


1X55981.1 


1272338.1 


1272338.1 


IAL1 61 588.2 


1196 


i 
g 

m 


lAi478339.1 


IAB00308&.1 


S80750.1 


(O 
UJ 

m 


IAJ245772.1 


IY12227.1 


L76301.1 


Most Similar 
(Top) Hit 
BLAST E 
Value 


? 

UJ 


8 

UJ 


00+31 


00+31 


? 

UJ 


1E+00 


00+31 


S 

UJ 


1E+00 


tE+00 


? 

UJ 


lE+00 


00+31 


1E+00 


1E+00 


s 

4- 
LU 


1E+00 


1E+00 


S 
♦ 

UJ 


s 

+ 

UJ 


S 
* 

UJ 


lE+00 


1E+00 


1E+00 


1E+00 


lE+00 


0 
? 

UJ 


























































in 
6 


5 




6 




d 


1.04 


1.52 


CO 

o 
a» 


3.23 


1.25 


0.71 


0-71 


d 


0.78 


2.08 


2.08 


9.72 


0.99 


CO 


0.75 


0.75 


0.73' 


0.68, 


d 


d 




Expressi 
Signal 




















































ORFSEQ 

ID NO: 


291451 




301371 


30138] 


30190 


302131 




1 305551 


1 31139] 


1 311611 


31621 


31807 


1 319311 


1 32225| 


1 327291 


1 328861 


1 328871 


1 32910] 


n 

N- 
CO 


1 33530 


1 33620 


1 34130 


34206 




1 35007 




3 

■n 
« 


Exon 
SEQ ID 
NO: 


166881 


16881 1 


1 17703| 


1 17704| 


17766 


1 17794| 


8 

s 


1 181431 


1 184231 


1 


10851 


19023 


1 191381 


00 
1 


1 19865| 


1 200231 


1 200231 


1 20041] 


24787 


to 

1 


1 207041 


I 212111 


21288 


1 20297 


1 22046 


8 




Probe 
SEQ ID 

NO: 


I 40931 


1 42951 


5131 


5132 


5201 


5230 




1 55101 




1 581 5j 


6242 


6420 


I 6539] 


1 68171 


1 7338! 


7501 


L 7501 


1 7521 


7560 


[ 8074 


1 8163 


1 8672 


8749 


1 9358 


% 


1 9599 


1 
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Top Hit Descriptor 


Xenopus laevis rhodopsin gene, complete cds | 


Agarteus bisporus mRNA for tyrosinase | 


Homo sapiens calcium channel alphel E subunit (CACNA1E) gene, exons 7-49, and partial cds, alternatively 1 
spliced 1 






Pilot whale morbilBvirus phosphoprotein (P) gene, partial cds | 


Oncorhynchus mykiss stl1 mRNA for rhamnose binding lectin STL1 . complete cds | 


Hordeum vulgare gere encoding cysteine proteinase | 


Bos taurus micromdar calcium activated neutral protease 1 (CAPN1) gene, exons 11-20, and partial cds 






z 

\^ 

o 

fXL 
Q. 

a: 

UJ 
CD 
U. 


o 

CO 

i 

g 

< 

Z 
O 

u 

M 

S 
1 

I 

v> 

i3 
3 

to 

Q. 
O 

o 
z 

M 

5 
6 

1 

m 

I 
5 




Homo sapiens cdl cydo protein (PA2G4) gene, axons 2 though 5 | 


SRB-11 PROTEIN | 


V.carteri gene encoding vdvcoct^sin | 


insulin-Bkegrcwvlh factor-binding protein 4 [cattle, pulmonary artery endothelial cells. mRNA. 2028 nt] 


B-CELL RECEPTOR CD22 PRECURSOR (LEU -14) (B-LYMPHOCYTE CELL ADHESION MOLECULE) 
(BL-CAM) 


Homo sapiens endothelin-converting enzyne 2 (ECE2) mRNA, complete cds | 


ac79b08.s1 Stratagene lung (#937210) Homo sapiens cDNA clone IMAGE:668791 3' | 


601443950F1 NIH_MGC_65 Homo sapiens cDNA done IMAGE:3846005 5' | 


601443950F1 NIH_MGC_65 Homo sapiens cDNA done IMAGE:3848005 5' [ 


Rattus norvegicus mRNA for N-acetylglucosamJny(transferase III, complete cds | 


PEROXISOMAL HYDRATASE-DEHYDROGENASE-EPIMERASE (HDE) (MULTIFUNCTIONAL BETA- 
OXIDATION PROTEIN) (MFP) [INCLUDES: 2-ENOYL.COA HYDRATASE ; D-3-HYDR0XYACYL COA 
DEHYDROGENASE ] 


(- « to 


NT 1 


NT 1 


NT 


NT 1 


IN 


IN 


NT 1 


NT 1 


NT 


NT 


NT 1 


SWISSPROT 1 


EST_HUMAN | 




NT 1 


H 
O 

q: 

0. 

to 
tS 

« 


NT 1 


NT 


SWISSPROT 


NT 1 


EST HUMAN 


EST_HUMAN 


EST HUMAN 


NT 


SWISSPROT 


i 

h 

± 


U23808.1 1 


AJ2238ia.1 1 


AF223391.1 


89222451 


D10852.1 1 


§ 


AB039022.1 | 


Z97022.1 1 


AF248054.1 


AF248054.1 


Z97341.2 1 


P04501 1 


AW452782.1 | 


U75902.1 1 


Ul 
< 


g 

in 
0. 


Y11204.1 1 


S52770.1 


P20273 


AF192531.1 


AA775191.1 


m 


i 

s 

so 
UJ 

m 


D1 0852.1 


Q02207 


Most Similar 
(Top) Hit 
BUSTE 
Value 


l.OE+00 1 


o 
o 
+ 
UJ 

o 


o 
? 

UJ 

p 


I.OE+OOl 


I.OE+OO] 


I.OE+OOl 


I.OE+OOl 


1.0E+OOI 


l.OE+00 


l.OE+00 


l.OE+00 1 


o 

? 

UJ 

p 


I.OE+OOl 


1.0E+OOI 


1.0E+00| 


I.OE+OOl 


I.OE+OOl 


s 

♦ 

UJ 

p 


l.OE+00 


l.OE+00 


l.OE+00 


1.0E+00 


1 

UJ 
a 


1 l.OE+00 


f 

UJ 

p 




0.78| 


1.551 




0.64| 


17.21 


to 
o 




2.561 


4.54 


4.54 


1.221 


■ 4.411 


■ 1.561 


1.791 


0.83 1 


in 


1.271 


1.22 




_ 1.36 


7.92 


1.49 


1.49 


1.28 


<M 


ORFSEQ 
ID NO: 




2W13I 


29189 








30394| 


30527| 


31368 


S 
« 

to 


314801 


316411 


316451 


320181 


320621 




324421 


32573 




33093 


33107 


33349 


33350 




33753 


Exon 
SEQ ID 
NO: 


127791 


163451 


16736 






179001 


180101 


181201 


18633 


18633 


187271 


- 16471 1 


18877| 




19258 




19608 


19724 


20016 


202061 


20219 


20444 


20444 


17734 


20832 




3659 


3744 


4144 


4362 


t 51651 


L ^ 


5432 


5486 


6013 


6013 


6111 


6283 


1 62691 


6615 


1 6662| 


1 67421 




7192 


7493 


1 7697 


1. 7710 


1 7902 


1 7902 


1 8084 


8291 
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ETA- 
COA 


JITIN- 1 


JITIN- 


JITIN- 




1 1 

it 










































JLTI FUNCTIONAL B 
)-3-HYDR0XYACYL 


O 
CQ 

LU 

« 


TERASEie) (UBIQI 
16) (UBIQUITIN 


ITERASE 16) (UBIQI 
16) (UBIQUITIN 




Jin (pd) gene, partial < 
i(vpu), Env protein (e 














S 










to 

lii 
O 

i 




CD 

§ 

CM 

lii 
< 


IMAGE:428906 5" 










Top Hit Descriptor 


IXISOMAL HYDRATASE-DEHYDROGENASE-EPIMERASE (HDE) (Ml 
\TION PROTEIN) (MFP) [INCLUDES: 2-ENOYL-COA HYDRATASE ; f 
DROGENASE ] 


JITIN CARBOXYL-TERMINAL HYDROLASE 11 (UBIQUITIN THIOLES 
;IFIC PROCESSING PROTEASE 11) (DEUBIQUITINATING ENZYME ■ 


UBIQUITIN CARBOXYL-TERMINAL HYDROLASE 16 (UBIQUITIN THIOLES 
SPECIFIC PROCESSING PROTEASE 16) (DEUBIQUITINATING ENZYME ' 
PROCESSING PROTEASE UBP-M) 


UBIQUITIN CARBOXYL-TERMINAL HYDROLASE 16 (UBIQUITIN THIOLES 
SPECIFIC PROCESSING PROTEASE 16) (DEUBIQUITINATING ENZYME ' 
PROCESSING PROTEASE UBP-M) 


< 

s 

u 
ta 

.1 
1 

o 

§ 
I 

§ 

o 

t- 
r 

<D 
1 

i 

o 
}- 


1 immunodoflciency virus Gag protein (gag) gene, complete cds; Pol prott 
rtein (vif). Vpr protein (vpr). Tat protein (tat). Rev protein (rev). Vpu proteir 
xMn (nef) genes. > 


n Immunodeficiency virus type 1 (HIV-1), Isolate SF33, 


17581 F1 NIH_MGC_70 Homo sapiens cDNA done IMAGE:3899421 5' 


nusculus chloride channel caldum activated 1 (CIcal). mRNA 


lusculus chloride channel calcium activated 1 (Cloal). mRNA 


J554 GKC Homo sapiens cDNA done GKCCYA1 1 5' 


>U5 laevis zona peltucida C glycoprotein precursor (xlZPC) mRNA, compl< 


JUS laevis zona pdlucicte C glycoprotein precursor (xlZPC) mRNA, compl< 


in Coronavirus gene for membrane protein 


in Coronavirus gene for membrane protein 


1 sapiens MHC binding factor, beta (MHCBFB) mRNA 


> sapiens MHC binding factor, beta (MHCBFB) mRNA 


s 

o 
< 
z 

1 

o 

§ 

I 
u. 

z 

1 

i 

u 

s 
i, 
i 


8823 BM Homo sapiens cDNA done BMFAWC04 5* 


02.r1 Soares_fetal_liver_spleen_lNFLS_S1 Homo sapiens cDNAdono 


i02.r1 SoQres_fetalJiver_spleen_lNFLS_S1 Homo sapiens cDNA done 


in retinoblastoma susceptibility gene exons 1-27, complete cds 


sprdine-rich protein (intron 3} [human. Genomic. 898 nt] 


sum vulgare gene encoding cysteine proteinase 


JMBOMODULIN PRECURSOR (FETOMODULtN)(TM) 






UBIQI 
SPEC 


RCU 


i 

« > 2 


X 


§ 




M 


> 
< 


1 




Huma 


X 


Homo 


X 




> 
< 




■6 


I 


PBR1 


Horde 


|THR< 


Top Hit 
Database 
Source 


SWISSPROT 


SWISSPROT 


SWISSPROT 


SWISSPROT 


|EST HUMAN | 


NT 


iNT 1 


|EST_HUMAN ! 






lEST HUMAN | 


1 iNl 


NT 


iNT 1 


h- 

z 


NT 


NT 


lEST HUMAN 


Z 

i 

I 

1- 

CO 
UI 


lEST.HUMAN 


|EST HUMAN 


1- 
Z 


NT 


iNT 


ISWISSPROT 


i 


















7534291 


7534291 












174562 


174562 


















Top Hit Ac 
No. 


Q02207 


P51784 


IQ9Y5T5 


Q9Y5T5 


[BE147331. 


U42720.2 


IM38427.1 


i 

8 

Ui 
CD 


CD 


CD 


§ 


|U44g52.1 


|U44gS2.1 


1X15498.1 


i 

X 




lO 


IAI077920.- 


IAV758825 


IAA004982 


IAA004982 


|L11910.1 


i 
w 


i 


1 PI 5306 


Most Similar 
(Top) Hit 
BLAST E 
Value 


00+30 


00+30 


00+30 


00+30 


00+30 


00+30 


00+30 


1 

UI 

o 


00+30 


f 00+30 


OE+00 


OE+00 1 


00+30 


00+30 


] 00+30 


00+30 


OE+OO 


00+30 


00+30 


00+30 


00+30' 


00+30 


o 
o 
+ 

LU 
O 


o 
o 
+ 
Ui 
o 


00+30 




















































Expression 
Signal 




«r> 

CD 

6 


0.48 


0.48 


1 2.171 


1.06 


1 1-07| 


1 2.051 


1 ''■34j 


1 1-34| 


CD 
O 

oi 


1 1.33] 


1 1.33] 


m 
d 


]g-o 1 


1 0.62] 


1 0.62 


1 0.73 


4.17 


1 19.78 


1 19.78 


1 0.93 


1.87 


1.57 


1 4.85 


ORF SEQ 

ID NO: 


33754 




33902 


33903 




1 33974 


1 341161 


L 346651 


1 348601 


1 34861 1 


1 34987] 


1 34993] 


1 34994 


? 

CM 

to 
r> 


1 35249] 


1 35510 


35511 


1 35592 


1 35716 


1 35859 


1 35860 


w 


1 36393 


t 30527 




Exon 
SEQ ID 
NO: 


20832 


20953 


20987 


20987 


247911 


\0 

o 

?^ 


21198! 


217121 


219111 


, 219111 


22028] 


] 22034] 


] 22034] 


] 22265] 


] 22265] 


1 22516] 


22516 


[22600] 


[ 22725] 




1 22866 


1 22898 


1 23374 


1 18120 


1 24201 


1 oo 


8291 


8413 


8447 


8447 


8475 


8513 


8659 


9195 


9402 


9402 


CO 
(9 


9534 


9534 


9767 


9767 


10021 


10021 


10105 


10230 


10372 


10372 


10404, 


10853' 


11587 


11837 
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Top Hit Descriptor 


EST368293 MAGE resequertces, MAGN Homo sapiens cDNA | 




Drosophila melanogaster regulator of G-prolein signalling LOCO III mRNA. complete cds | 


Homo sapiens chromosome 21 segment HS21 C1 02 | 


{Apple mosaic vims RNA 2 putative polymerase gene, complete cds | 


SERINEH-HREONINE PROTEIN KINASE MINIBRAIN ] 


^PROBABLE OXIDOREDUCTASE 2K1290.5 IN CHROMOSOME II | 


i Lycopersioon esculentum putative Mil copy 1 nematod&resistance gene | 


|B2 BRADYKININ RECEPTOR (BK-2 RECEPTOR) | 


1 Dank) rerio mRNA for Eph-llke receptor tyrosine kinase rtkS | 


lAMlNO^ACID ACETYLTRANSFERASE (N-ACETYLGLUTAMATE SYNTHASE) (AGS) (NAGS) | 


1 CallithrK jacchus UBE1 gene derived retroposon on the Y chromosome ] 


IXenopus laevis rac GTPase mRNA, complete cds | 


IPROBABLE ENDONUCLEASE IV (ENDODEOXYRIBONUCLEASE IV) t 


I601653583R2 N1H_MGC_55 Homo sapiens cDNA clone IMAGE:3838461 3' j 


I601653583R2 NIH_MGC_.55 Homo sapiens cDNA clone IMAGE:3838461 3' ] 


Enterobacteriaceae sp. JM983 partial gnaES gene for GroES-tike protein and partial groEL gene for GroEL- 
likeprolefn. Isolate JM983 


Enterobacteriaceae sp. JM983 partial groES gene for GroES-like protein and partial groEL gene for GroEL- 
[ike protein, isolate JM983 


io 
at 

i 

o 

8 
< 

Z 
Q 

u 

M 

E 

t 

O 

i 

X 

5 
X 
Z 

fe 

CO 

1 
§ 


I601456337F1 NIH_MGC_66 Homo sapiens cDNA clone IMAGE:3860049 5" J 


IPHOSPHOGLUCOMUTASE (GLUCOSE PHOSPHOMUTASE) (PGM) ] 


[od55d04.s1 NCI_CGAP_GCB1 Homo sapiens cDNA clone IMAGE:1 371847 3' j 


)601110258F1 NIH_MGC_16 Homo sapiens cDNA clone IMAGE:3350750 5* j 


1601110258F1 NIH_MGC_16 Homo sapiens cDNA clone 1MAGE:3350750 5' | 


CO 
CM 

ill 

i 
§ 

o 

i 
t 

o 

g 

X 
< 

8, 

O 

z 

o 
o 
(M 

s 


Homo sapiens X28 region near ALD locus containing dual specificity phosphatase 9 (DUSP9), ribosomal 
protein L18a (RPLISa), Ca2+/Calmodultn-dependent protein kinase 1 (CAMKI), creatine transporter (CRTR). 
CDM protein (CDM). adrenoieukodystrophy protein > 


iDrosophila melanogaster sodium channel protein (para) gene, exons 9,10,1 1,12 and optional segments b, c, dl 
land e. partial cds | 


[Triticum aestivum stripe rust resistance protdn Yr10 (YMO) gene, complete cds | 


[Salmonella typhimurium adenlne-methyttransferase (mod) and restrictjon endonuclease (res) ] 


Top Hit 
Database 
Source 


EST HUMAN | 


NT 1 


NT 1 


INT 1 


iNT 1 


ISWISSPROT 1 


ISWISSPROT 1 


iNT 1 


ISWISSPROT 1 


iNT 1 


ISWISSPROT 1 


|NT 1 


|NT 1 


ISWISSPROT 1 


|EST_HUMAN | 


|EST_HUMAN | 


NT 


NT 


Z 

i 

X 
Ui 


|EST_HUMAN | 


ISWISSPROT 1 


lEST^HUMAN 1 


[EST_HUMAN 


]EST HUMAN 


lEST_HUMAN 


NT 


NT 


iNT 


iNT 


Top Hit Acession 
No. 


AW976184.1 1 


AF245455.1 | 


< 


AL1 63302.2 | 


AFl 74585.1 I 


P49657 1 


Q09632 ] 


U65667.1 1 


Q28642 1 


AJ005O29.1 1 


P22567 1 


AJ003108.1 1 


1 

< 


067551 1 


BE957439.2 | 


BE957439.2 | 


AJ302158.1 


AJ302158.1 


IBF034016.1 1 


IBF034016.1 1 


IP38652 1 


[AA825565.1 i 


IBE258705.1 


]BE258705.1 


1AI680876.1 


U52111.2 


U267161 


JAFI 491 12.1 


1M90544.1 


Most Similar 
(Top) Hit 
BLAST E 
Value 


I.OE+OOl 


9.gE-0l| 


9.9E-01I 


9 


9.9E-01I 


9.9E-01 1 


9 

UJ 
o> 
a 


9.9E-01| 


9 
a> 


9.9E-01I 


9.8E-01I 


9 

UJ 

00 

0> 


9.8E-01 1 


9.8E-01I 


9 

O) 


9 
a> 


9 


S 

00 

a> 


1 9.8E-01I 


9 
o> 


1 9.8E-01I 


1 9.8E-01I 


1 9.8E-01 1 


9 


1 9.8E-01I 


9.8E-01 


9 

UJ 

r- 
a 


9 

Hi 

a 


] B.7E-01 


"sigrtal 


3.0B| 


&» 
b 


0.971 


1.171 


0.94| 


! 14.591 


! 0 83| 


1391 


ci 


1 1.681 


1 1.771 


1 0.89| 


in 
o 
ci 


1 0.95| 


1 0-61 1 


1 0.61 1 


4.86 


4.86 


1-13 


1 1.131 


1 0.771 


0.56 


L_... 4.86 




1.78 


? 


2.51 




1 1.28 


ORF SEQ 
ID NO: 




267421 


267431 


27794| 




n 


31389] 






361421 


256581 






289301 


289331 


28934| 


32634 


32635 


1 330381 


1 33039 


1 341101 




1 364101 


1 364111 


1 371091 




32597 


1 33895 


] 33901 


Exon 
SEQ ID 

NO: 


CM 


14209 


14209 


15222 


16267 


18440 


18648 


217021 


218731 


231281 


13180] 


14907] 


153791 


164671 


1 164701 


16470| 


19779 


19779 


1 201531 


1 20153] 


1 21192] 


1 22828] 


23395 


23395 


) 24040] 


24341 


19743 


1 20980 


[ 20986 














Probe 
SEQ ID 
NO: 


! 121711 


1 16161 


1616 


i 


3665 


5816 


L 60291 


1 91851 


9474 


1 105931 




2336 


2827 


§ 

CO 


3872 


3872 


7250 


7250 


1 7641 1 


1 7641 1 


CP 


o 


1 10874] 


1 10874 1 


] 115971 


12058 


1 7212 
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30/S26 



wo 01/57277 



PCT/USO 1/00669 




31/526 



wo 01/57277 



PCT/US01/()«669 




32/526 



wo 01/57277 



PCTAJSO 1/00669 
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wo in/fiiin 
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(7) 



Top Hit Descriptor 


1 
1 

a 

z 
« 

"a 

1 

i 
1 

3 

it 

q < 

« — • 

n 
II 


Homo sapiens PR01975 mRNA. complete cds | 


Cotumix cotumix japonica sub-species japonica beta^tin mRNA, partial cds | 


Columix cotumw japonica sub-species japonica bet&^actin mRNA. partial cds | 


RAFFINOSE INVERTASE (INVERTASE) | 


RAFFINOSE INVERTASE (INVERTASE) 1 


yf24b02.s1 Soares fetal liver spleen 1 NFLS Homo sapiens cDNA clone IMAGE;1 27755 3' 1 




Archaeoglobus fulgidus, complete genome | 


Or>ctolagus cuniculus immunogtobultn VDJ region gene j 


Orvctdagus cuniculus immunoglobulin VDJ region gene | 


cds 


cds 


MATING-TYPE PROTEIN A-ALPHA Z4 1 


aq14b12.x1 Stanby Frontal NS pod 2 Homo sapiens cDNA clone IMAGE:2030879 ( 


aq14b12.x1 Stanley Frontal NS pool 2 Homo sapiens cDNA clone IMAGE:2030879 | 




Mus musculus neuromedin U precursor (Nmu) gene, partial cds; tPhlP (Tphip) gene, partial cds; CLOCK 
(Qock) gene, complete cds; PFT27 (Pft27) gene, complete cds; and H5AR (H5ar) gene, complete cds 


Mus musculus advillin (Advil-pending). mRNA | 


Mus musculus advillin (Advil-pending), mRNA | 


GLUTAMATE [NMDA] RECEPTOR SUBUNIT EPSILON 3 PRECURSOR (N-METHYL D^PARTATE 
RECEPTOR SUBTYPE 2C) (NR2C) (NMDAR2C) 


GLUTAMATE [NMDA] RECEPTOR SUBUNIT EPSILON 3 PRECURSOR (N-METHYL D-ASPARTATE 
RECEPTOR SUBTYPE 2C) (NR2C) (NMDAR2C) 


Mus musculus cytochrome P450, 2b9, phenobarbitol inducible, type a (Cyp2b9), mRNA | 


MUSCARINIC ACETYLCHOLINE RECEPTOR M2 | 


MUSCARINIC ACETYLCHOLINE RECEPTOR M2 1 


1 H.espersa mRNA for neurofilament NF70 | 


1 H.espersa mRNA for neurofilament NF70 | 


Top Hit 
Database 
Source 


NT 


NT ! 


NT 1 


NT 1 


SWISSPROT 1 


SWISSPROT j 


EST_HUMAN 1 


NT 1 


IN^ 


in! 


NT 1 


NT 


NT 


SWISSPROT 


EST_HUMAN 


EST HUMAN 


NT 1 


NT 


INl 


NT 


SWISSPROT 


SWISSPROT 


NT 


SWISSPROT 


SWISSPROT 




NT 


TopHitAcession 
No. 


8393408 


AF1 18085.1 1 


AF199488.1 | 


AF1M4e8.1 1 


PI 6553 1 


PI 6553 1 


R08600.1 1 


AB021 134.1 1 


114976211 


L27316.1 1 


L27316.1 1 


AF059510.1 


AF059510.1 




AI253399.1 1 


AI253399.1 | 


U72487.1 1 


AF146793.2 


68577521 


68577521 


Q01098 


Q01098 ' 


6753577 


CN 

s 

2 


P30372 




1X86347.1 


Most Similar 
(Top) Hit 
BLAST E 
Value 


7.7E-01 


7.7E-011 


7.7E-01 1 


7.7E.01I 


7.7E-01 1 


7.7E-011 


7.7E-01I 


7.7E-01 1 


7.7E-01J 


7.6E-011 


7.6E-01| 


7.6E-01 


7.6E-01 


7.6E-0lj 


7.6E-0lj 


9 

UJ 
(D 


s 

to 


7.5E-01 


7.6E-01 


7.6E-01 


7.6E-01 


7.6E-01 


7.6E-01 


7.6E-01 


7.6E-01 


1 7.6E-01 


1 7.6E-01 


Expression 
Signal 


0.62 


4.781 


3.171 


3.,7| 


1.331 


1.331 


CO 

6 


0.821 


15.011 


19.731 


19.731 


4.81 


4.81 


1^ 
o 


0.951 


0,95| 


0.98| 


1.34 


1.76 


1.76] 


0.55 


0.55 


0.91 


3.33 


3.33 


1 2.74 


1 


ORF SEQ 
ID NO: 




28734 1 


295211 




310841 


310851 


314791 


352391 




29822 


298231 


31628 


31629 


320391 


304651 


30488] 


322691 


33455 


335221 


335231 


33727 


33728 


34359 


34664 




1 ' ' 36824 


1 36825 


Exon 
SEQ ID 
NO: 


16009 


162621 


170711 


170711 


18375| 


183751 


187261 


222561 


24280] 


173701 


173701 


18857 


18857 


192371 


180431 


180431 


<o 
in 


20551 


20610 


20610 




20808 


21435 


21720 


1 21720 


g 


1 23767 


SEQ ID 
NO: 


3400 


1 36601 


1 


4486 


1 57491 


1 57491 


1 61101 


1 97581 


1 119571 




4790 


6248 


6248 


i 


6935 


6935 


1 71131 


8009 


1 80681 


S 


1 


8287 


1 


1 9203 


1 9203 


1 11236 


1 11236 
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o OD Q- 
^ -S liJ 
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O o Q. 



1.3 



2 



o 
E 
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I 

it 
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o 
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; 5. 



1 15 



^ 1 ' 
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Top Hit Descriptor 


Pdyangium vitellinum (strain PI vtl ) 1 eS rRNA gene 1 


C 

1 

a. 

§ 

.§ 

i 
< 

z 

cc 

E 

M 
O 

ce: 


6021 3931 9F1 NIH_MGC_46 Homo sapiens cDNA clone IMAGE:42981 17 5' | 


V(rl39b12.x1 NCI_CGAP_Ut1 Homo sapiens cDNA clone IMAGE:2427263 3' | 


TRANSCRIPTION-REPAIR COUPLING FACTOR (TRCF) | 


IL2-BT0731-250400-077.G08 BT0731 Homo sapiens cDNA 1 


AV712326 DCA Homo sapiens cDNA clone DCAAUF07 5" | 


yig4a0g.$1 Soares placenta Nb2HP Homo sapiens cDNA clone IMAGE:146872 3' | 


QV4-ST0023-1 60400-1 72-a01 ST0023 Homo sapiers cDNA | 


< 
z 

u 

1 

O 

X 

i 

? 

o 

§ 

S 

I 

a 


Human regereratng protefn (reg) gene, complete cds 1 


65B1 Human retna cDNA Tsp509l-cleaved suWibrary Homo saptens cDNA not directional [ 


Human carboxyl ester lipase (CEL) gene, complete cds | 


601 556863F1 NIH MGC 58 Homo sapiens cDNA clone IMAGE:3826767 5* | 


nac51f10.x1 NCI_CGAP_Bm23 Homo sapians cDNA ckme IMAGE:3406218 3' similar to contains element 


|Homo sapiens postmeiotic segregation increased 2-like 9 (PMS2L9), mRNA 1 


|Homo sapiens postmeiotic segregation increased 2-like 9 (PMS2L9), mRNA | 


Buchnera aphidicola genomic fragment containing (chaperone Hsp60) groEL, DNA biosynthesis initiating 
protein (dnaA), ATP operon (atpCDGAHFEB). and putative chromosome replication protein (gidA) genes. 


Buchnera aphrdicola genomic fragment containing (chaperone Hsp60) groEL, DNA biosynthesis initiating 
protein (dnaA). ATP operon (atpCDGAHFEB). and putative chromosome replication protein (gidA) genes, 
complete cds; and termination factor Rho(rho) gene> 


|Thermotoga marilinia section 97 of 1 38 of the complete genome | 


1 Musmusculusanti-DNA immunoglobulin light chain IgM mRNA, antibody 363p. 136, partial cds | 


iHomo sapiens mRNA for KIAA1184 protein, partiat cds 1 


[Xenopus laevis smooth muscle beta-tropomyosin mRNA, complete cds | 


I601823850R1 NIH_MGC_79 Homo sapiens cDNA clone IMAG&4043485 3' 


Ui 

lii 
O 

i 

o 

8 

"0 
o 

M 

8 

1 

O 

g 
I 

% 

X 

z 

CO 

§ 


GLYCOGEN DEBRANCHING ENZYME (GLYCOGEN OEBRAN CHER) [INCLUDES: 4-ALPHA- 
GLUCANOTRANSFERASE (OUGO-1.4-1.4-GLUCANTRANSFERASE); AMYLO-1,6-GLUCOSIDASE 
(DEXTRIN &-ALPHA-D-GLUCOSIDASE)] 


Top Hit 
Database 
Source 


NT 1 


NT 


EST_HUMAN | 


EST HUMAN | 


SWISSPROT 1 


EST_HUMAN | 


EST HUMAN | 


EST_HUMAN | 


EST HUMAN | 


EST HUMAN | 


INT 1 


lEST HUMAN | 


iNT 1 


lEST.HUMAN | 


EST HUMAN 


iNT 1 


NT 1 


NT 


NT 


NT 


iNT 


INl 


NT 


[EST HUMAN 


|EST_HUMAN 


SWISSPROT 


1 
































N 
(O 
U> 


s 




















Top Hit Ace 
No. 


i 

SJ 
5 


X87885.1 


BF683095.1 


AI85849S.1 




BE091 796.1 


AV71 2326.1 


R80e73.1 


AW806881. 


i 


J05412.1 


iW22302.1 


IM94579.1 


IBF030207.1 


BF439982.1 




1 


AF0Q8210.1 


AF008210.1 


o 

< 


lU 55574.1 


1AB033010.1 


IM92304.1 


IBF107848.1 


IBF317212.1 


P35573 


Most Simitar 
(Top) Hit 

Value 


5.1E-01I 


5.1E-01I 


5.1E-01 


S.IE-OI] 


S.IE-Oll 


5.1E-0lt 


5.1E-01I 


5.1E-01I 


5.lE-0li 


5.1E-01| 


5.iE-qiJ 


5.1E-01I 


9 

UJ 

ui 


1 S.IE-Oll 


5.1E-01 


1 5.0E-01I 


1 S.OE-Oll 


5.0E-01 


9 

UI 

o 

IT) 


1 5.0E-01 


1 S.OE-01 


1 5.0E-01 


1 5.0E-01 


1 5.0E-01 


9 

UJ 

o 

ui 


5.0E-01 


Expreeston 
Signal 


3.981 


0.88 


11.33 




3.031 


0.71 1 


0.791 


1.42 


0.73| 


0.73| 




CO 


v> 
at 
b 


I 2.041 


2.01 






5.46 


5.46 


5.58 


1 0.65 


3.11 


1.78 


t 
b 


n 


1.34 


ORFSEQ 
ID NO: 


25782 






291971 


293001 






323161 


339661 


1 339671 


1 350651 


1 35066] 


1 355551 






I 27326! 


1 273271 


27335 


27336 


28809 


1 26875 


29006 




34057 


33219 


3S004 




13300 




1 


16743 


16852 


S 


19025 


lO 

1 


21046 


o 


22102 


1 221051 


1 225601 


24805 


24385 


14757 


14757 


14765 


14765 


16341 


16410 


16540 


21007 


21143 


20318 


22043 


— 

SE 
















































Probe 
SEQID 
NO: 


fe 

CD 


1692 




5 


4266 


5229 


6422 


6997 


I 85071 


8507 


9602 


n9605 


1 100651 


11874 


12129 


s 


2180 




2189 


3740 


3811 


<N 


3 


8ec4 


9379 
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Top Htt Descriptor 


601306121F1 NIH_MGC_39 Homo sapiens cDNA clone IMAGE:3e40420 5' I 


Homo sapiens hyaluronan synthase 2 (HAS2), mRNA | 


W46t01.s1 Soares infant brain 1NIB Homo sapiens cDNA clone IMAGE;35639 3' I 


< 
Z 
0: 
E 

'0 
a. 

i 
i 

c 

s 
1 
1 

1 

1 
0. 

i 

V) 


602124743F1 NIH_MGC_56 Homo sapiens cDNA clone IMAGE;428161 1 5' | 


602124743F1 NIH_MGC_56 Homo sapiens cDNA clone IMAGE:428161 IS" | 


qieiel 1.x1 NCI_CGAP_Kid3 Homo sapiens cDN A clone IMAGE:1 863980 3' similar to gb;S55700 J 
HYDROXYMETHYLGLUTARYL-COA LYASE PRECURSOR (HUMAN); | 


yb47h08.Bl Stratagene fetal spleen (#837205) Homo sapiens cDNA clone IMAGE:743e7 3' similar to similar 
to gb:M91036_ma2 HEMOGLOBIN GAM^M-A AND GAMMA-G CHAINS (HUMAN) 


601883592F1 NIH_MGC_57 Homo sapiens cDNA clone IMAGE:409581 4 5' 


Homo sapiens KIAA0764 gene product (KIAA0764). mRNA | 


Andis opalinus isolate QS NADH dehvdrogenase subunit 2 (ND2) gene, complete cds; mitochondrial gene 
for mitochondrial product 


Stizostedion vitreum 40S ribosomaf protein SI 1 mRNA, partial cds | 


Homo sapiens membrane-bound amlnopeptldase P (XNPEP2) gene, complete cds 


Homo sapiens transcription factor IGHM enhancer 3, JM1 1 protein. JM4 protein, JM6 protein, T54 protdn, 
JM10 protein. A4 differentiation-dependent protdn. triple LIM domain protdn 6. and synaptophysin genec. 
complete cds; and L-type calcium channel a> 


Mus musculus peptidogl^n recognition protdn-like (Pglyrpl-pending). mRNA | 


Mus musculus protein kinase C, epsilon (Pkce), mRNA | 


Homo eapiens Xq pseudoautosomal region; segment 1/2 | 


CO 

? 

eo 

CM 

lii 

i 
§ 

0 

n 

u 

s 
t 
§ 

X 

< 
0 

Z 


Balaenoptera physalus gene encoding atria! natriuretic peptide 1 


Corvnebacterium sp. ALY-I alyPG gene for pdyguluronale lyase, complete cds I 


PM1-ST0262-261199-001-g01 ST0262 Homo sapiens cDNA | 


Homo sapiens Xq pseudoautosomal region; segment 2/2 1 


Balaenoptera physalus gene encoding atrial natriurelic peptide | 


601 5949eOF1 NIH_MGC_9 Homo sapiens cDNA clone IMAGE;3948734 5' | 


Cantagalo orlhopoxvirus hemagglutinin gene, complete cds | 


RC3-BT03 33-1 80700-111-803 BT0333 Homo sapiens cDNA \ 


RC3-BT0333-180700-l11-a03 BT0333 Homo sapiens cDNA ] 


Mus musculus 12WBvClara cdt 10 kd protdn (mCCiO) gene, completecds | 


Mouse cytokeratin 1 5 gene, complete cds { 


Top Hit 
Detatxase 
Source 


z 

=) 
I 

K 

W 
LU 


NT 1 


EST HUMAN | 


NT 1 


EST_HUMAN | 


EST HUMAN | 


EST HUMAN 


EST HUMAN 


z 

1 

X 
1- 

w 

lU 


NT 


NT 


NT ] 


NT 


NT 


NT 1 


NT ] 


NT ] 


EST HUMAN ] 


NT 


NT 


EST_HUMAN ] 


NT ] 


NT 1 


EST HUMAN ] 


NT ] 


z 
< 

X 
UJ 


EST_HUMAN ] 


NT ] 


INl 


Top Hit Acession 
No. 


BE737392.1 i 


48853901 


R45318.1 1 


66793221 


BF696639.1 | 


BF696639.1 1 


0 

1 

< 


T55325.1 


BF216117.1 1 


7662291 j 


a 

i 

< 


AF304162.1 1 


AF195953.1 ] 


AF1 96779.1 


109466231 


67550831 


AJ271735.1 1 


AW300400.1 1 


AJ008755.1 1 


AB030481.1 ] 


AW 81 7785.1 1 


AJ271736.1 1 


AJ006755.1 1 


BE741629.1 | 


AF229247.1 | 


BE693575.1 | 


BE693575.1 | 


U01247.1 ] 


016313.1 ] 


MostSimflar 
(Top) Hit 
BLAST E 
Value 


3.IE-O1I 


3.IE-OII 


3.IE-OII 


3.IE-OII 


9 

UJ 

w 


3.IE-OII 


3.1E-01 


9 

UI 

<o 


9 

UJ 

« 


3.IE-OI] 


9 

UJ 

n 


3.1E-01] 


3.IE-OI] 


3.1E-01 


9 

LU 

ei 


3.0E-01] 


3.0E-01 ] 


9 

LU 
0 
CO 


3.0E-01] 


9 

Ul 

0 
pi 


3.OE-OII 


3.0E-01] 


3.OE-OII 


3.0E-01 ] 


3.0E-01] 


3.OE-OI] 


3.0E-01] 


3.0E-01] 


3.OE-OI] 


a c 
a 0) 


w 


1.18| 




0.45| 


0.81] 


0.81 1 




0.54 


1.95] 


2.03] 


2.13 


1.9S] 




3,89 


1.62] 


1.37] 


11.51] 


2.05] 


s 




cvi 


1.01] 


1.78] 




0.88) 


4.03] 


4.03] 


4.57] 


3.14] 


ORF SEQ 
ID NO: 


304581 


3307 1| 


34041] 


352951 


354561 


354571 


35525 




i 

« 


s 












25234] 


25419] 


26380] 


26680] 




28997] 


29109] 


29640] 


30601 1 


30726] 


30820] 


30821] 


31062] 


32407] 


Exon 
SEQ ID 
NO; 


24741 1 


20183 


21121 


22314 


22474 


22474] 


22530 




232451 


238621 


24257 


24282] 


24370] 


24624 


25035] 


5382 


2932 


3863 


4147 


5860 


6530 


6643 


7194 


8186 




8321 


8321 


8357 


9578 1 






























Probe 
SEQ ID 
NO: 


1 


1 7671 1 


1 8582| 


9816 


1 


1 9979] 


10035 


10204 


10717] 




1 


1 


! 12108] 


1 

CM 


\ 12535] 




in 

h~ 
tN 


! 1266] 


1 1555] 


32*8 


I 3932 


^ 

9 


1 4611] 


1 5554] 


1 5627] 


m 


1 5695] 


1 5731 1 


] 69191 
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Top Hit Descriptor 


Choristoneura fumiferana diapause associated protein 2 (DAP2) mRNA, complete cds | 


I Vibrio cholerae chromosome II, section 73 of 93 of the complete chromosome | 


endogenous retrovirus) element 




starfish (P.ochraceus) cytoplasmic ecdn gene, complete cds | 


1 T3 receptor-associating cofaetor-1 [human, fetal liver. mRNA. 2930 nl) | 


[Homo sapiens chromosome 21 segment HS21 C007 ] 


Homo sapiens partial sleerin-1 gene | 


Feline caliciviais CFiy68 RNA helicase/cysteine protease/RNA-dependent RNA polymerase pdyprotein 
precursor and capsid protein precursor, genes, complete cds; and unknown gene 


1 Mus musculus SKD1 (Skdl ) gene, complete cds 1 


Homo sapiens chromosome 21 segment HS21C082 1 


CO 

1 
CO 

i 
i 

< 

u 

i 
1 

o 
I 

« 

Q. 

< 

a. 

1- 
O 

1 

u. 

IL 

m 

2 

a 
o 

1 

rt 

s 

IO 

p: 


I601653391R2 NtH_MGC_58 Homo sapiens cDNA clone IMAGE:3826198 3' | 


|601459238F1 NIH.MGC_66 Homo sapiens cDNA done IMAGE:3862809 5' | 


E1B *=ROTEIN. SMALL T-ANTIGEN (E1B 19K) 


|yq84f07,r1 Soares fetal liver spleen 1NFLS Homo sapiens cDNA ck>ne IMAGE:202501 5* | 


[Mouse tfisiis-speciftc protdn (TPX-1 ) gene, exon 1 0 | 


1 Homo sapiens matrix nnetalloprateinase MMP Rasi-1 gene, pronwter region | 


jHomo sapiens matrix metalloproteinase MMP Rasi-1 gene, pronroter regkxi | 


1 Hordeum vulgare receptor-like kinase LRK1 0 gene, partial cds | 


1 Hordeum vulgare receptor-like kinase LRK1 0 gene, partial cds | 


o 

s 

t 

o 

g 

X 

s 

^ 
9 
o 
o 

o 

rt 

i 

rt 
O 
OC 


xg40cl0.xl NCI_CGAP_Ut1 Homo sapiens cDNA clone IMAGE:2e30034 3* similar to contains Alu repetitive 


1 Mouse LIMd LINE DNA | 


1 Arabidopsis thaliana DNA chromosome 4, contig fragment No. 1 7 | 


1 Arabidopsis thaliana DNA chromosome 4, contig fragment No. 1 7 ] 


] Human mRNA for KIAA01 24 gene, partial cds | 


] Zea mays cellulose s vnthasB-4 (CesA-4 ) mRNA, complete cds . | 


] Arabidopsis thaliana DNA chromosome 4. contig fregment No. 41 ( 


]on70d04.s1 Soares_NFL_T_GBC„Sl Homo sapiens cDNA clone IMAGE:1 562023 3' 1 


]602132442F1 NIH_MGC_81 Homo sapiens cONA done IMAGE:4271578 5' | 


Top Hit 
Database 
Source 


iNT 1 


iNT 1 


NT 


|EST_HUMAN 1 


iNT 1 


iNT 1 


iNT 1 


iNT 1 


NT 


iNT 1 


NT 


lEST HUMAN | 


|EST HUMAN | 


|EST_HUMAN | 


ISWISSPROT 1 


|EST_HUMAN | 


iNT 1 


iNT 1 


iNT 1 


iNT 1 


iNT 1 


|EST HUMAN | 


EST HUMAN 


iNT 1 


1 INl 


|NT 1 


]NT 


]NT 


]NT 


]EST HUMAN 


z 
< 

D 
T 

h- 

to 

UJ 


TopHttAcession 
No. 


IAF007768.1 | 


IAE004416.1 1 


AJ230113.1 


BE896785.1 | 


IM26501.1 1 


|S63390.1 1 


|AL1 63207.2 | 


IAJ251 973,1 1 


Ul 3992.1 


IAF134119.1 1 


lALI 63282.2 1 


|BF109040,1 ! 


IBE960712.1 j 


IBF03B595.1 j 


IP04492 j 


IH53238.1 1 


i 


IU89651.2 1 


IU89651.2 1 


IAF085164.1 1 


IAF085164.1 1 


IAW581997.1 1 


AW1 52248,1 


1X58491.1 1 


IAL161 505.2 1 


IAL161505.2 1 


ID50914.1 1 


IAF200528.1 | 


IAL161541.2 1 


1AA936316.1 


IBF576124.1 


Most Simitar 
(Top) Hit 
BLAST E 
Vaiue 


1 2.5E-01I 


1 2.5E-011 


2.5E-01 


1 2.5E-01 1 


t 2.5E-01 1 


3 

IO 
C\i 


1 2.5E-01 1 


I 2.5E-01I 


2.5E-01 


1 2.5E-01I 


[ 2.5E-01I 


[ 2.5E-01 1 


1 2.5E-01I 


1 2.5E-01 1 


1 2.5E-01I 


i 2.5E-01I 


! 2.5E-01 1 


1 2.5E-01I 


1 2.5E-01I 


1 2.5E-01I 


9 
oi 


1 2.5E-01I 


2.5E-01 


1 2.5E^1| 


1 .2.5E^1| 


1 2.5E-0l| 


] 2.5E-01 1 


1 2.5E-01 1 


] 2.5E-01 1 


] 2.4E-0l| 


] 2.4E-01I 


Expression 
Signal 


1 4,691 


1 2.82| 


3.21 


1 0.61 1 


0.65| 


12.861 


s 

d 


0.831 


0.82 


1 1.351 


_____ 4.48 


1 2.311 


r«- 
d 


OI 

oi 


0.72i 


1 3,03| 


1 o.8e| 


1 15.98| 


1 15.981 


a> 
o 
oi 


1 2.091 


CD 
(D 


1.53 


1 1.311 


1 2,03| 


1 2.03| 


rt 


1 5.291 


1 10.131 


1 1.97 


OI 


ORF SEQ 

ID NO: 


1 298961 


1 29927 1 




299761 




305761 




CO 


32778 




! 330251 


[ 33236 1 


i 

rt 
rt 


33625 1 


337981 


1 340291 


1 34271 1 


1 348091 


1 349101 


1 348971 


1 

rt 


354961 


35937 


1 35941 1 


359621 


rt 

§ 

rt 


i 

rt 


1 37153] 




1 25687 1 


1 26014| 


Exon 
SEQ ID 
NO: 


174451 


17471 1 


17499 




129331 


18161 1 


18730] 


CI 

§ 


19914 




201441 


20329] 


203401 




208771 


21110| 


213471 


21961 1 


s 

OI 


219481 


1 

o( 


225051 


22930 


229331 


22953] 


22953] 


23470] 


24122] 


1 25075] 


1 13209] 


1 13495] 


Probe 
SEQ ID 
NO: 


1 


s 


4924 


9 


1 52621 


to 
in 


6114 


CO 

n 


7389 


7413 


1 76321 


I 77861 


1 77971 


8168 


8336 


1 8571! 




1 94351 


ID 

a> 


I 94921 


CM 

a 


[ 100101 


10436 


1 104391 


1 10459] 


1 10459] 


1 10955) 


1 11712] 


1 11740] 




i 
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Top Hit Descriptor 


1 
» 

a 

3 

CO 


C. parasitica espC gene | 


and 11-16 


Q V4-EN0032-1 90500-223-e03 EN0032 Homo sapiens cONA | 


Homo sapiens hypothetical protein ASH1 (ASH1). mRNA | 


HOMEOBOX PROTEIN GLABRA2 (HOMEOBOX^LEUCINE ZIPPER PROTEIN ATHB-10)(HD.ZIP 
PROTEIN ATHB-10) 




Homo sapiens dual OMdase-Bko domains 2 (DU0X2), mRNA i 


F.rubripes DNA encoding for valyl-tRNA synthetase | 


Saccharomyces cerevlsiae HaISp (HAL5) mRNA. complete cds | 


Human hepatocyte growth factor gene, eyon 1 | 


M.auratusmu class glutathione transferase gene | 


PM1-CT0247-141099-001-g06 CT0247 Homo sapiens cDNA 1 




Mus musculus phosphcfructokinase-l 0 isozyme (Pfkc) gene, exons 3 through 7 1 


Homo sapiens mRNA for FLJ00016 protein, partial cds | 


Andes virus strain 0123133 glycoprotein G1 and G2 precursor, gene, partial cds 1 


M. musculus &cp2 gene axon 14 | 


601344648F1 NIH_MGC_8 Homo sapiens cDNA clone IMAGE:3677794 5' | 


Dfctyostelium discoideum random slug cDNA19 protein (rsc19) mRNA. partial cds 1 


Arabtdopsis pyruvate decartxix^ase-2 (Pdc2) gene, complete cds | 


Chlamydia trachomatis section 5 of 87 of the complete genome 1 


DAUGHTERLESS PROTEIN | 


z 

i 

a: 

UJ 
— J 

on 
in 
(- 

X 

0 


Homo sapiens fllamfn 2 (FLN2) mRNA. complete cds | 


Arabidopsis thaDana root gravitropism control protein (PIN2) gene, complete cds | 


Arabidopsis thaliana root gravitropism contrd protein (PIN2) gene, complete cds | 


Homo sapiens cAMP specific phosphodiesterase (PDE4C) gene, exons 2 through 12 | 


Homo sapiens cAMP specific phosphodiesterase <PDE4C) gene, exons 2 through 12 | 


D melanogaster DNA mobile element (hoppel) | 


R.norveglcus mRNA for NTR2 reoeptor ) 


Salvelinus pluvtus mRNA for transferrin, complete cds ] 


Top Hit 
Database 
Source 


NT 1 


in| 


NT 


EST_HUMAN | 


NT 1 


SWISSPROT 


NT 1 


NT 1 


NT 1 


NT 1 


NT 1 


NT i 


EST HUMAN 1 


NT 


NT 


NT 1 


INl 


NT 1 


z 
< 

X 

h- 
(0 
UJ 


INl 


NT 1 


NT 1 


SWISSPROT 1 


SWISSPROT 1 


NT 1 


NT 1 


NT 1 


NT 1 


NT 1 


NT 1 


NT 1 




TopHitAcesslon 

No. 




X83997.1 1 


AF242431.1 


BE826165.1 1 


89220801 


P46607 


X56600.1 1 


114325401 


X91856.1 1 


b 

i 


M75967.1 1 


X61033.1 1 


AW360665.1 | 


U39724.1 1 


AF250371.1 1 


AK02+»27.1 1 


AF028026.1 | 


X91151.1 1 


r*: 

i 

UJ 
CD 


U82511.1 1 


U71122.1 1 


AE001278.1 1 


PI 1420 1 


PI 1420 1 


si 

1 

LL 
< 


AF086907.1 | 


AF086907.1 | 


AF1 57814.1 1 


AF1 57814.1 1 




X97121.1 1 


D89088.1 1 


MostSimtla- 
(Top)Hil 
BLAST E 
Value 


2.0E-01 1 


2.OE-OII 


2.0E-01 


9 
fM 


2.0E-01 1 


9 

<Ni 


2.OE-OII 


2.OE-OII 


9 

UJ 

0 


2,0E-01| 


2.OE-OII 


2.OE-OII 


2.OE-O1I 


2.OE-OII 


2.OE-O1I 


2.OE-O1I 


9 


2.OE-O1I 


2.0E-01I 


9 

Ul 

0 

ci 


2.OE-OII 


9 


2.OE-O1I 


2.OE-OII 


2.OE-OII 


2.0E-01 1 


2.0E-01 1 


2.OE-O1I 


2.OE-OII 


2.0E-01 1 


2.OE-O1I 


9 

N 


Expression 
Signal 


0.781 


0.681 


0.76 


8.43| 


7.091 


0.62 


2.381 


CM 


0.69| 


© 


0.71 1 


^ 


3.63| 


0.68| 


1.18| 


1.531 


6.45] 


2.91] 


0.53] 


1.03] 


0.65] 


4.35] 


0.51 j 


0.51] 


1.98] 


1.791 


1.791 


0.531 


0.53| 


0.72| 


0.88] 


2.77] 


ORFSEQ 
ID NO: 


.290981 


291 52 1 


29552 




30188] 


28621 


30737 1 


312631 


31 361 1 


31 591 1 




319551 


32049] 


8 


32727] 


33181] 




33595] 






34775] 




35132] 


35133] 




35431 1 


354321 


35562] 


35563] 




35789] 


38283] 


Ewjn 
SEQ ID 
NO: 


16626] 


(O 

S5 
S 


17106 


17247 


17757 


16139 


•A 


185381 


186261 


188201 


1 
3 


19158] 


19246] 


19780] 


19863] 


20284] 


20437] 


20883] 


21107] 


217901 


21826] 


21874] 


22160] 


22160] 


22304] 


22449] 


22449] 


22567] 


22587] 


22610] 


22798] 


23248] 


Probe 
SEQID 
NO: 


? 


\ 41021 


4522 


i 46851 


i 5192| 


5226 


5638 


5916 


(D 
0 

s 


6210 


!3 


6560 


6650 


7251 


7336 


7775 


7895 


8142 


: 8658] 


9273 


: 9312 


; 9475 


i 


9681 


CD 

s 


9954 


9954 


10072] 


10072 


10115] 


10304] 


10720] 
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wo 01/57277 PCT/USOl/00669 
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wo 01/57277 



PCT/US01/0(t669 
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wo 01/57277 PCT/USOl/00669 
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